Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB143	5,394,361	(G)6→7	coding (188/1731 nt)	CJ019_05073 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,394,355	1	.	G	92.3%	123.9 / 4.6	40	coding (182/1731 nt)	CJ019_05073	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (16/20);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GCATCAGAAAGTTGCGCGAGCTGATTCTGGAACTGGCGGTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGC‑GGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATACCCAATCGCT  >  ACB143/5394226‑5394485                                                                                                                                   |                                                                                                                                   gCATCAGAAAGTTGCGCGAGCTGATTCTGGAACTGGCGGTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGC‑ggggg                                                                                                                               <  1:288770/135‑1 (MQ=255)   aTCAGAAAGTTGCGCGAGCTGATTCTGGAACTGGCGGTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGC‑ggggggg                                                                                                                             >  4:365696/1‑134 (MQ=255)         aaGTTGCGCGAGCTGATTCTGGAACTGGCGGTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCa                                                                                                                       <  1:189370/136‑1 (MQ=255)                                    gctgTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCg                                                                                              <  5:326746/131‑1 (MQ=255)                                    gCGGTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCg                                                                                              <  4:128005/134‑1 (MQ=255)                                    gCGGTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCg                                                                                              <  5:290231/134‑1 (MQ=255)                                       gTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCgtag                                                                                           >  7:146583/1‑134 (MQ=255)                                        tGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGCGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGtagtag                                                                                        <  8:244701/136‑1 (MQ=255)                                        tGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGtagtag                                                                                        <  4:216781/136‑1 (MQ=255)                                        tGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGtagtag                                                                                        <  7:283554/136‑1 (MQ=255)                                                  gTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCt                                                                              <  2:212138/136‑1 (MQ=255)                                                       tACAGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACa                                                                         <  6:21983/136‑1 (MQ=255)                                                          cGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAgcg                                                                      >  2:95536/1‑136 (MQ=255)                                                            cAAGAGCGTAATGATGAGCCGGCCAGTGAGGTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCt                                                                    <  7:210871/136‑1 (MQ=255)                                                            cAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTATCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCt                                                                    <  6:224996/136‑1 (MQ=255)                                                            cAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCt                                                                    <  6:341116/136‑1 (MQ=255)                                                            cAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCt                                                                    <  1:289287/136‑1 (MQ=255)                                                            cAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCt                                                                    <  3:261528/136‑1 (MQ=255)                                                            cAAGAGCCTAATGATCAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGAGTAGTAGCTAAAGGTCTAGACAGCGCt                                                                    <  3:20831/136‑1 (MQ=255)                                                              agagCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGC‑GGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGaa                                                                 >  2:137255/1‑136 (MQ=255)                                                                 gCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACt                                                               <  3:137649/136‑1 (MQ=255)                                                                     aaTGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTAcc                                                            >  1:149809/1‑135 (MQ=255)                                                                       tgatgaGCCGGCCAGTGAGTTTCTAGCGCGAATTGCATTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGAAAAGAAAGGAAATGTTAAAGCTGGAGTAGTAGCGAAAGGTCTAGACAGAGATTAACTACCca                                                          >  8:307233/1‑135 (MQ=255)                                                                              ccGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTg                                                   >  7:33653/1‑135 (MQ=255)                                                                              ccGGCCAGTGAGTGGCTAGCGCGAATTGCCTTAGAGAGTGCACTGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTgg                                                  <  2:102013/136‑1 (MQ=255)                                                                               cGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGa                                                 >  5:341898/1‑136 (MQ=255)                                                                               cGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGa                                                 >  1:333135/1‑136 (MQ=255)                                                                                  ccAGTGAGTGGATAGAGCGAATTGCCTTTGAGACTTCAAAGATAGTGGCTGGGGTGTAATGCAAAGAAAGGAAATGTTAAATCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGatat                                              <  4:274463/136‑1 (MQ=255)                                                                                    aGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATg                                             >  1:279054/1‑135 (MQ=255)                                                                                     gTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGC‑GGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGatat                                              <  2:81409/132‑1 (MQ=255)                                                                                     gTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGC‑GGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGatat                                              >  1:81409/1‑132 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTATTAGCGCAAGGTCTAGACAGAGCTGAAATACCCACAACTTGTATATGTACATCGTTCGATGATCTAATAAATCCAGAATAccc         <  7:262605/136‑1 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTTGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATAccc         <  3:69699/136‑1 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTTGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATAcc          >  4:69699/1‑135 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATAccc         <  1:166042/136‑1 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATAcc          >  5:324645/1‑135 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATAcc          >  7:330269/1‑135 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATAcc          >  2:350742/1‑135 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATAcc          >  2:14120/1‑135 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCAAGAATaac          >  2:198615/1‑133 (MQ=255)                                                                                                                       acaGCTAGTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGAGGATCTAATAAATCCAGAATAcc          >  6:100576/1‑135 (MQ=255)                                                                                                                              gTGGCGGGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATACCCAATCGCt  >  8:337922/1‑136 (MQ=255)                                                                                                                                   |                                                                                                                                   GCATCAGAAAGTTGCGCGAGCTGATTCTGGAACTGGCGGTGCGCGGTAAGTTGGTACCGCAAGAGCCTAATGATGAGCCGGCCAGTGAGTTGCTAGCGCGAATTGCCTTAGAGAGTGCACAGCTAGTGGC‑GGGGGGTAAGGCAAAGAAAGGAAATGTTAAAGCTGGCGTAGTAGCGCAAGGTCTAGACAGCGCTGAACTACCCACAACTTGGATATGGACATCGTTCGATGATCTAATAAATCCAGAATACCCAATCGCT  >  ACB143/5394226‑5394485

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑