Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,626	0	G	A	41.5%	47.0 / 15.8	41	F180F (TTC→TTT)	ycfT	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (17/7);  new base A (0/17);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 2.53e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACC  >  NC_000913/1174479‑1174775                                                                                                                                                     |                                                                                                                                                      accGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaaa                                                                                                                                                    <  1:667699/149‑1 (MQ=255)       aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAcc                                                                                                                                              >  2:818813/1‑151 (MQ=255)           aaaGAAGCCATCAGCAGATGGCGGCGCAAGGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑GGTAATAATGCAAATTGCGGATCACACGGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCca                                                                                                                                          <  1:521958/151‑1 (MQ=255)           aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGcc                                                                                                                                           >  1:870515/1‑150 (MQ=255)                    aacagcagATGGCGGCCCACCGGCCCCTTTTTAACACAGGTCATTATTTTTGCCCCAAACCAGGCGCCAAGGC‑TGTAAAAAGGAAAATGGCGGATCACCCTTTTCTTCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTcagca                                                                                                                                 <  1:756255/149‑1 (MQ=255)                        gcagATGGCGGCGCAACGGCCCCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTTTTCATTCCCCACCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTac                                                                                                                             <  2:764197/151‑1 (MQ=255)                           gATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaa                                                                                                                           >  2:418832/1‑150 (MQ=255)                                        cggcccttttttAACACAGGTCATTATTTTTGCCCCAAACAATG‑GCCCAGGCTTGTAAGAGGAAAAATTGCGCATCACCCTGTTCATCCCCCCCCACGGTGGGGGAACGAAATTAACCGCCACACTCAGCAGTCCAAACAAGGCGAAt                                                                                                                <  1:811616/139‑1 (MQ=255)                                          gCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGCACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTgg                                                                                                           >  1:699194/1‑151 (MQ=255)                                                ctTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgcgc                                                                                                     >  1:201350/1‑151 (MQ=255)                                                     aagcaGGTCATTTTTTTGGCCAAAAACATTGCCCCAAGCG‑TGAAATAAGGCAAATGGCGGAACCCACGGTCCATCCCCCACCAAGGGGGGGAAAGAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                                         <  1:392926/139‑1 (MQ=255)                                                      cacaGGTATTTATTGTCGCGCCAACCCAGGCGCCAGGGC‑TGAAATAAGGAAAATGGCGGATCACACTGTTCATTCCCCACCACGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                                               <  1:307836/151‑1 (MQ=255)                                                          gttCATTTTTGTTGCGCAAACCCATGCGCAAAGGT‑TGTAATAAGGCAAATTGCGGATCACACTTTTCTTTCCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                                               <  2:899563/145‑1 (MQ=255)                                                           gTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTa                                                                                           >  2:592244/1‑150 (MQ=255)                                                            tCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                                          >  2:138015/1‑150 (MQ=255)                                                              attattGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                                        >  1:330186/1‑150 (MQ=255)                                                                tattGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                                     >  2:223316/1‑151 (MQ=255)                                                                  ttgttgCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                                                   >  2:726030/1‑151 (MQ=255)                                                                  ttgttgCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                                                   >  1:263556/1‑151 (MQ=255)                                                                    gttgCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAc                                                                                 >  2:548476/1‑151 (MQ=255)                                                                     tggcgcAAACCCAGGCGCAAGGGC‑TGAAAAAAGGAAAATTGCGGATCACCCTGTTCATCCCCCCCCAGGGGGTGGAAACGAAATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAc                                                                                 <  1:650498/148‑1 (MQ=255)                                                                       gcgcCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                                               >  2:216035/1‑150 (MQ=255)                                                                        cgcAAAACCATGCGCAAAGGT‑TGTAATAAGGAAAATTGCGGATCACCCTTTTCTTCCCCCCCCAGGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAc                                                                             <  2:870746/151‑1 (MQ=255)                                                                                                 aataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACAATCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                                      >  1:444652/1‑150 (MQ=255)                                                                                                 aaaaaGGCAAATTGCGGATCACCCTGTTCATCCCCCCCCACGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACg                                                         <  2:409138/144‑1 (MQ=255)                                                                                                    aaGGAAAATGGCGGATCACCTTGTTCTTCCCCCCCCACGGCGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACa                                                   <  1:444993/150‑1 (MQ=255)                                                                                                      ggaaaaTTGGGAACCACCCTTTTAATCCCCCACCAGGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGa                                                 <  1:630645/147‑1 (MQ=255)                                                                                                       gaaaaTTGCGGATCCCACTGTTCATTCCCCACCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                                                <  1:463546/148‑1 (MQ=255)                                                                                                         aaaTTGCGGATCACACTTTTCTTCCCCCCCCCCGGGGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                                               <  1:230218/149‑1 (MQ=255)                                                                                                          aaTTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                                             >  1:50474/1‑150 (MQ=255)                                                                                                             tGCGGATCACACTTTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                                          <  2:539217/150‑1 (MQ=255)                                                                                                                           ttCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCg                            <  2:330285/150‑1 (MQ=255)                                                                                                                               ttCCCCCCCGGGGCGTGGAAACAAAATGACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                          <  2:59631/148‑1 (MQ=255)                                                                                                                                 ccccccCGCGGCGGGGGAACAAAATTAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg                               <  2:558434/141‑1 (MQ=255)                                                                                                                                 ccccccCAGGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                          <  1:762280/146‑1 (MQ=255)                                                                                                                                     cccagggggggggAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATccc                 <  1:167823/141‑1 (MQ=255)                                                                                                                                         cggggggggAACAAAATGACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATccc                 <  1:477736/141‑1 (MQ=255)                                                                                                                                         cGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGt               <  2:155299/149‑1 (MQ=255)                                                                                                                                               ggaaaCGAAATAATCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAAc       <  1:518617/148‑1 (MQ=255)                                                                                                                                                  aaCGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTc     >  2:254794/1‑150 (MQ=255)                                                                                                                                                    cGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCAcc  >  1:254406/1‑151 (MQ=255)                                                                                                                                                     gAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCAcc  >  2:545643/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGC‑TGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACC  >  NC_000913/1174479‑1174775

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