Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,358,236	A→G	G300G (GGA→GGG)	thrA →	homoserine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,358,236	0	A	G	97.3%	132.2 / ‑2.0	37	coding (900/1281 nt)	thrA	homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/17);  minor base . (1/0);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

GAATATATGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGAGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAGATAAAACACGGGAAATGATGGATTCAG  >  CP000730/1358091‑1358372                                                                                                                                                  |                                                                                                                                         gAATATATGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGgc                                                                                                                                        >  2:42856/1‑148 (MQ=255)  aaTATATGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGgca                                                                                                                                       >  1:324966/1‑148 (MQ=255)       atGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGgcaggcag                                                                                                                                  <  2:295994/148‑1 (MQ=255)       atGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGgcaggcag                                                                                                                                  <  1:325762/148‑1 (MQ=255)       atGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGgcaggcag                                                                                                                                  <  2:406396/148‑1 (MQ=255)       atGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGgcaggcag                                                                                                                                  <  1:96315/148‑1 (MQ=255)         gAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGtt                                                                                                                                <  1:330505/148‑1 (MQ=255)          aaaaTGGCAAAGTTAATTCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGttt                                                                                                                               >  1:381785/1‑148 (MQ=255)          aaaaTGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGttt                                                                                                                               >  1:412189/1‑148 (MQ=255)          aaaaTGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGttt                                                                                                                               >  1:205607/1‑148 (MQ=255)                       tAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGt                                                                                                                  >  2:237340/1‑148 (MQ=255)                           gCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCg                                                                                                              <  1:54882/148‑1 (MQ=255)                               cGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGt                                                                                                          <  1:112301/148‑1 (MQ=255)                                    gAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTg                                                                                                     >  2:171053/1‑148 (MQ=255)                                      aCCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGAt                                                                                                   <  1:25058/148‑1 (MQ=255)                                      aCCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGAt                                                                                                   <  1:367674/148‑1 (MQ=255)                                              tAATTGATAAAAAGCATCAATTAGGAGGTGTAGAGGATGAATATAACGCGATTTCTGTTATTGGTGATGCCGTTGGTGACACGCTGTTTTATGGAAAGGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATATATTGaa                                                                                           >  2:374242/1‑148 (MQ=255)                                                   gATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAg                                                                                      >  1:377129/1‑148 (MQ=255)                                                   gATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAg                                                                                      >  1:201081/1‑148 (MQ=255)                                                           gCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTAGGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTc                                                                              <  2:270107/148‑1 (MQ=255)                                                                  ttAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAAt                                                                       >  1:233272/1‑148 (MQ=255)                                                                  ttAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAAt                                                                       <  2:296315/148‑1 (MQ=255)                                                                         cTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGAttt                                                                >  1:10254/1‑148 (MQ=255)                                                                          tGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTa                                                               <  1:76538/148‑1 (MQ=255)                                                                                    gAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGc                                                     >  2:271328/1‑148 (MQ=255)                                                                                       tataACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGccac                                                  <  2:205607/148‑1 (MQ=255)                                                                                       tataACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGccac                                                  <  2:412189/148‑1 (MQ=255)                                                                                       tataACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGccac                                                  <  2:381785/148‑1 (MQ=255)                                                                                         taACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGccacca                                                >  1:322453/1‑148 (MQ=255)                                                                                           aCGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAGGG‑GCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACAt                                             >  2:354525/1‑148 (MQ=255)                                                                                           aCGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACcaca                                              >  1:70256/1‑148 (MQ=255)                                                                                                            ggTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAg                             >  1:410762/1‑148 (MQ=255)                                                                                                            ggTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAg                             >  1:289522/1‑148 (MQ=255)                                                                                                            ggTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAg                             >  2:409965/1‑148 (MQ=255)                                                                                                                 tGCCGTTGGTGACACGATGTTTTATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAGATaaa                        <  1:84969/148‑1 (MQ=255)                                                                                                                                     tttATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAGATAAAACACGGGAAATGATGGATTc    <  1:371405/148‑1 (MQ=255)                                                                                                                                       tATGGAAAAGGGGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAGATAAAACACGGGAAATGATGGATTCAg  >  1:29742/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                         GAATATATGAAAATGGCAAAGTTAATGCATCGGTAGAACCAACGTTAATTGATAAAAAGCATCAATTAGCAGCTGTAGAGGATGAATATAACGCGATTTATGTTATTGGTGATGCCGTTGGTGACACGATGTTTTATGGAAAAGGAGCAGGCAGTTTAGCAACAGGTAGTGCCGTTGTCAGTGATTTATTGAATGTAGCATTATTCTTTGAATCAGATTTACACACATTGCCACCACATTTTGAATTAAAGACAGATAAAACACGGGAAATGATGGATTCAG  >  CP000730/1358091‑1358372

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