Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,139,881	C→T	V578M (GTG→ATG)	asmA ←	putative assembly protein AsmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,139,881	0	C	T	96.8%	103.6 / NA	31	V578M (GTG→ATG)	asmA	putative assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/13);  minor base G (0/1);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGATTGACCTTGCATGTCGTTTAACGTCACGACGCCATCCTTCAACGTCAAATCGGT  >  NC_000913/2139643‑2140121                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 aaGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGTCTGTa                                                                                                                                                                                                                                      <  2:362177/251‑1 (MQ=255) aaGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGt                                                                                                                                                                                                                                       >  1:367880/1‑250 (MQ=255)     gAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTg                                                                                                                                                                                                                                  <  2:315232/251‑1 (MQ=255)      aaCGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGa                                                                                                                                                                                                                                 <  1:281681/251‑1 (MQ=255)                          tACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCGTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATTTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATGCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCATTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGc                                                                                                                                                                                                                  <  1:75312/246‑1 (MQ=255)                                          aaTGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTg                                                                                                                                                                                             >  2:99328/1‑251 (MQ=255)                                                          aCGCCGGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGaa                                                                                                                                                                             <  2:357113/251‑1 (MQ=255)                                                              cTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAaatc                                                                                                                                                                          >  1:169714/1‑250 (MQ=255)                                                                             acccgTCTTATCAGGCCTACAATACTCGTCACCTGGTAACCTACATCTTCTCCCGCAACTTCTTTACATATTTGCCATTACGGGACTCTTTATTCCGCTCGGCCCAGTCATTCAGGCGCCGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTCCTGGCGTTTCTATCAGGAAATCAATCAGTTTGCTTTcc                                                                                                                                                          >  2:173627/4‑251 (MQ=255)                                                                                cgTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCGTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCg                                                                                                                                                       <  2:292208/251‑1 (MQ=255)                                                                                cgTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCAGTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGGTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCg                                                                                                                                                       <  1:6911/251‑1 (MQ=255)                                                                                 gTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGt                                                                                                                                                      >  1:131824/1‑251 (MQ=255)                                                                                      aTCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTcc                                                                                                                                                   >  2:234748/1‑249 (MQ=255)                                                                                                  ggACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGAcc                                                                                                                                     <  1:270400/251‑1 (MQ=255)                                                                                                         tCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAatat                                                                                                                              >  1:352867/1‑251 (MQ=255)                                                                                                                     cTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGgtgt                                                                                                                  >  1:317281/1‑251 (MQ=255)                                                                                                                       aCATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGgtgt                                                                                                                  <  1:217390/249‑1 (MQ=255)                                                                                                                           cttctCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAgg                                                                                                            >  1:374474/1‑251 (MQ=255)                                                                                                                                      aCTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGcc                                                                                                 <  1:143438/251‑1 (MQ=255)                                                                                                                                         tcttcACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGcc                                                                                                 >  1:53383/1‑248 (MQ=255)                                                                                                                                             cacaTCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTc                                                                                           >  1:96662/1‑250 (MQ=255)                                                                                                                                               caTCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGCTTCAAc                                                                                        >  1:79994/1‑251 (MQ=255)                                                                                                                                                                         ttCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGcac                                                              >  1:97828/1‑251 (MQ=255)                                                                                                                                                                                 ggCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGAtt                                                      <  2:367880/251‑1 (MQ=255)                                                                                                                                                                                     cAGTCATTCAGGCGACGTTTCGCTTCGTCCTGGAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGATTGAcc                                                  >  1:210436/1‑251 (MQ=255)                                                                                                                                                                                      aGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGATTGACCt                                                 >  2:16082/1‑251 (MQ=255)                                                                                                                                                                                                               atacTGTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTATCCCCGTTCCAGCCACCCACGTCCCGAATATCTAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAATATGCCTTCCCCTGTCAGCGCCAGCACAGACGATTGACCTTGCATGTCGTTTAACGTCACGACGc                        >  2:297110/4‑251 (MQ=255)                                                                                                                                                                                                                    gTAGATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGATTGACCTTGCATGTCGTTTAACGTCACGACGCCATcc                   >  2:325160/1‑251 (MQ=255)                                                                                                                                                                                                                    gTAGATGTTGGCGCAGTAACTGAGCCAGTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCAGAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGATTGACCTTGCATTTCGTTTAACGTCACGACGccttc                    <  2:14759/250‑4 (MQ=255)                                                                                                                                                                                                                       gATGTTTGCGCAGTAACTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCATGTTTGATCTGCCAGATTCATCATGCTTTCCCCTGTCAGCGCCAGCACAGGCTATTGACCTTGCATGTCGTTTAACGTCACGACGCCATCCTTc                <  1:372135/251‑1 (MQ=255)                                                                                                                                                                                                                                     aaCTGATCCATTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGATTGACCTTGCATGTCGTTTAACGTCACGACGCCATCCTTCAACGTCAAATCggt  <  1:206607/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 AAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTGGGTGTCGCAGGTTTGATCTGCCAGATTCAACATGCCTTCCCCTGTCAGCGCCAGCACAGGCGATTGACCTTGCATGTCGTTTAACGTCACGACGCCATCCTTCAACGTCAAATCGGT  >  NC_000913/2139643‑2140121

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