Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,748	0	C	A	44.2%	37.1 / 29.9	43	E55* (GAG→TAG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (16/8);  new base A (0/19);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 3.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTT  >  NC_000913/3177602‑3177895                                                                                                                                                    |                                                                                                                                                    ttCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCaaa                                                                                                                                                  <  2:21598/150‑1 (MQ=255)    aCCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt                                                                                                                                              >  2:390020/1‑151 (MQ=255)        tctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTccc                                                                                                                                          >  1:42172/1‑151 (MQ=255)         ctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGTTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc                                                                                                                                           <  1:28300/149‑1 (MQ=255)           tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccg                                                                                                                                         <  2:629576/149‑1 (MQ=255)            cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg                                                                                                                                      >  2:209623/1‑151 (MQ=255)                        cAACCATCCCCAGTAGCCAGGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGG‑TACCAAGACTGGGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGa                                                                                                                            <  1:715599/149‑1 (MQ=255)                                 ccAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGACAACCCGAGTCTGCTCAATAAGCACAACTTCATCGCGCCCTGGGTCAAAGGG‑AAGCAAGACTGCGGGGTGACCGCGCTCAAAATTTTCCCGCCGCACCTCATGACTCAttt                                                                                                                     <  2:181748/147‑1 (MQ=255)                                              gtttccgcgGGTTTCGTCCGCGAAAACCCAAACCTCCTAAACCACCCCACCTTCACCACCAACGGGGTAAAAGGG‑AACCAAGCCTGCGGTGTGACCGCGAAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa                                                                                                    <  1:160581/146‑1 (MQ=255)                                                 ttCGTTGGTGTCGTACGCGGCAATCCGAATCTGCCCAATCAGCACACCTTCACCACGCACTGGGTCAAAGGG‑TAGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                  <  2:323469/150‑1 (MQ=255)                                                 tgcGCTGGTGTCGTCCGCGGCAACCGGAATCTGCCCAACCACCACAACTCCATCACGCACTGGCCAAAAGGG‑TACCAAGACTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                  <  1:81383/148‑1 (MQ=255)                                                   cGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCACCACGCCCTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTCCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                                 <  1:499277/149‑1 (MQ=255)                                                    gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                               >  2:356987/1‑150 (MQ=255)                                                      tGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa                                                                                             >  2:579867/1‑150 (MQ=255)                                                          gtCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                         >  2:670733/1‑150 (MQ=255)                                                             gtccgcgCCAATCCAAACCTGCTCAACCACCACAACTTCATCCCCCACTGGGAAAAAGGT‑TACCAAGACGGGGGCGTCCCCGCTCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt                                                                                           <  1:532992/142‑1 (MQ=255)                                                             gTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                      >  1:251066/1‑150 (MQ=255)                                                             gTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACGGGGTAAAAGGG‑TAGCAAGACTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat                                                                                       <  1:416586/149‑1 (MQ=255)                                                                   ggCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGT‑TAGCAAGACTGCGGGTTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                                <  2:648753/150‑1 (MQ=255)                                                                    gCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt                                                                               >  2:357628/1‑150 (MQ=255)                                                                     cAATCCGAATCTGCTCAATCACCACACCTCCATCACGCACTGGGTAAAAGGG‑TAGAAAGACTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa                                                                              <  1:389974/150‑1 (MQ=255)                                                                      aaTCCGAATCTGCTCAAACAGCACAACTCCACCACGCACTGGG‑CAAAGGGTTAGCAAGACTGCGGCGTCACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt                                                                            <  1:81752/151‑1 (MQ=255)                                                                       aTCCGAATCTGCTCAATCAGCACAACTTCACCCCCCCCTGGGTCAAAGGG‑TAGCAAGACTGCGGTGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                           <  2:519421/151‑1 (MQ=255)                                                                         ccaaaTCTGCTCAATCAGCACAACTTCATCCCCCACTGGGTAAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaa                                                                         <  2:794465/148‑1 (MQ=255)                                                                          cGAATCTGCTCAATCAGCACTACTTCATCACGCACTGGGTTAAAGGG‑TAGCATGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGCTGAATAGACGATGATGAATTCTGTAAAGATCTAATg                                                                           >  1:644689/1‑148 (MQ=255)                                                                                    cAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                              >  1:764236/1‑151 (MQ=255)                                                                                     aaTCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                              <  1:771121/150‑1 (MQ=255)                                                                                     aaTCAGCACAACTTAATCCCCCACTGGGTAAAGGGG‑TAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                              <  1:39710/150‑1 (MQ=255)                                                                                            acaACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtg                                                      >  2:536007/1‑151 (MQ=255)                                                                                             caccTTCATCACCCACTGTGTCAAGGGG‑TAACAAGACTGCGTCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtg                                                      <  1:622653/147‑1 (MQ=255)                                                                                             caACTTCATCCCGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                     <  2:248500/151‑1 (MQ=255)                                                                                                      cgcGCTCTGGGTCAAGGGG‑TAGAAAGACTGCGGCGGACCCGGGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                   <  2:629568/142‑1 (MQ=255)                                                                                                           aCTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGCCCGCGCAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                         <  2:759252/149‑1 (MQ=255)                                                                                                             tGGTTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAtt                                     >  2:682747/1‑151 (MQ=255)                                                                                                              gggTCAAGGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAttt                                    <  2:785838/151‑1 (MQ=255)                                                                                                              gggTCAAAGGG‑TAGAAAGACTGCGGGGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaata                                        <  1:737589/147‑1 (MQ=255)                                                                                                                           aGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                     >  2:541331/1‑122 (MQ=255)                                                                                                                             cAAGACTGCGGGGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttttt                       <  2:225495/150‑1 (MQ=255)                                                                                                                              aaaacggcggggTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                     <  1:541282/113‑1 (MQ=255)                                                                                                                                      ggcgTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGt              <  1:98250/150‑1 (MQ=255)                                                                                                                                       gcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa             >  1:549229/1‑150 (MQ=255)                                                                                                                                                cgcTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGt   >  2:671281/1‑151 (MQ=255)                                                                                                                                                 gcTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGtt  >  2:797050/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTT  >  NC_000913/3177602‑3177895

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