Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	41.2%	32.4 / 13.4	36	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (14/6);  new base A (0/14);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 4.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAG  >  NC_000913/1222448‑1222741                                                                                                                                                      |                                                                                                                                                   cAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGa                                                                                                                                                   >  1:89797/1‑150 (MQ=255)         gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaa‑ca                                                                                                                                             <  2:204522/147‑1 (MQ=255)         gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAaa‑caaaca                                                                                                                                         <  1:313192/151‑1 (MQ=255)          cattacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAg                                                                                                                                        <  1:609930/151‑1 (MQ=255)                      ccGCAACATTTTCTGCAACCTGTGCGCTGTCAGTTTTTCCCTTTTCGTTGGTTTCTGCGAAAGTAATATGGGTTTCTGACGGAAAACCT‑TTTTTAATATTGTAGCCCTGCCCCCAGGTCATTGGTGCGAAAA‑AAAACAGGCCCGCCAGTa                                                                                                                             <  2:266490/150‑1 (MQ=255)                      ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTa                                                                                                                             >  2:619754/1‑150 (MQ=255)                       cGCAACATTTTCTGCACCCTGTGCTCTGTCATTTTTTGCCATTTCGTTGGTTTCTGCGCAAGTAATAGAGGTTTCTGACGGAAAACCG‑TTTTTAATTTTGTAGTCCTGCGCCCAGGTCATTGGGGCAAAAA‑AAAACAGGCCCGCCAGTaaa                                                                                                                           <  2:378459/151‑1 (MQ=255)                        gCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTaa                                                                                                                            >  1:583530/1‑149 (MQ=255)                             attttCTGCACCCTGTGCGCTGTCATTTTTTGCCATTTCGTTGGCTTTGGCGCAAGTAATATAGTTTTCTGACGGCAACCCT‑TTTTAAATATTGTAGTCCGGCGCCCGGTTCATTGGTGCAAAAA‑AAAACAGGCCCGCCAGTAAAGc                                                                                                                         <  2:353028/145‑1 (MQ=255)                               tttCTGCACCCTGTGCGTTGTCATTTTTTGCCTTTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGAAACCCT‑TTTTTAATATTGTAGTCCTGCGCCCAGTTCATTGGTGAAAAAA‑AAAACAGGCCCGCCAGTAAAGCAAAtt                                                                                                                    <  1:485278/150‑1 (MQ=255)                                      aaCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                            >  1:253165/1‑151 (MQ=255)                                           gtgCGCTGTCATTTTTTGCCTTTTCTTTGGTTTCTGCGCAATTAATATAGTTTTCTGACGGAAAACCT‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑AAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcatc                                                                                                         <  1:640970/149‑1 (MQ=255)                                           gtgCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                       >  2:552112/1‑151 (MQ=255)                                             gcgtTGTCTTTTTTGGCCTTTTCTTTGGTTTCTGCGCAAGTAATATAGGTTTCTGACGGCAACCCG‑TTTTTAATATTGTAGCCCTGCCCCCAGTTTATTGGTCCAAAAA‑CAACAAGGCCCGCCAGTAAAGCAAATTTTTtcatcatca                                                                                                        <  2:455594/148‑1 (MQ=255)                                               gcTGTCATTTTTTGCCTTTTCTTTTGTTTCTGCTAAAGTAATATGGTTTTCTGACGGCAAACCT‑TTTTAAATATTGTAGCCCTGCCCCCAGGTATTTGGTGAAAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcatca                                                                                                        <  2:630516/146‑1 (MQ=255)                                                    ctttttttGCCTTTTCGTTGGCTTCTGCGCAATTAATATAGTTTTCTGACGGAAAACCT‑TTTTTAATATTGTAGTCCTGCGCCCAGTTCTTTGCTGCAAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                   <  2:56120/143‑1 (MQ=255)                                                       tttttGCCTTTTCGTTGTTTTCTGCGCAAGTAATATAGTTTTCTGACGGAAACCCT‑TTTTAAATATTGTAGTCCTGCGCCCAGGTCTTTGGTGCAAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAtt                                                                                               <  1:576442/147‑1 (MQ=255)                                                                tttCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                    >  1:572928/1‑149 (MQ=255)                                                                    gTTGGCTTCTGCGCAAGTAATATAGGTTTCTGCCGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                              >  2:140540/1‑151 (MQ=255)                                                                          ttCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacc                                                                         >  1:148184/1‑150 (MQ=255)                                                                                     tAATATAGTTTTCTGCCGCCAACCTT‑TTTTTATTATTGTAGCCCCGCCCCCAGTTCATTGGGGCGAAAA‑AAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGc                                                              <  1:95597/150‑1 (MQ=255)                                                                                         ataGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATc                                                         >  1:82380/1‑151 (MQ=255)                                                                                                   gACGGCACACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGACATTGGTGGGAAAA‑CCAACAGACCCGCCAGTAAAGCATATTTTTTCATCATCATCCCTTATTTCTTTTTACCCAGACTTGCACCTCCCGTACCGCCCATCACTGCCCCtt                                               >  2:388304/1‑151 (MQ=255)                                                                                                      cgCAAACCCTTTTTTAATTTTTTAGTCCTCCCCCCAGTTCTTTGGTGAAAAAA‑CAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                              <  2:67857/149‑1 (MQ=255)                                                                                                                  tttAATATTGTAGCCCGGCCCCCAGTGCTTTGTCGAAAAAA‑CAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                  <  1:401638/150‑1 (MQ=255)                                                                                                                         ttgtgttCCTGCCCCAAGGTATTGGGTGCAAAAA‑AAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                 <  1:635048/138‑1 (MQ=255)                                                                                                                         ttGTATCCCTGCCCCAAGGTATTGGGTGAAAAAA‑AAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCa                           <  2:89780/150‑1 (MQ=255)                                                                                                                         ttGTAGCCCTGCCCCAAGTTCATTGGTGCAAAAA‑AAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGa                                <  1:672875/145‑1 (MQ=255)                                                                                                                         ttGTAGCCCGGCCCCCGGGTCTTTGGTGCAAAAA‑AAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                            <  1:71366/149‑1 (MQ=255)                                                                                                                          tgtatgcccggcccccGGTCCTTGGTGCCAAAACCAACCGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCa                           <  2:395159/138‑1 (MQ=255)                                                                                                                                 cTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                  >  1:644009/1‑151 (MQ=255)                                                                                                                                      cccAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc             >  2:149635/1‑151 (MQ=255)                                                                                                                                           gTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc          >  2:411515/1‑149 (MQ=255)                                                                                                                                              aTTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg      >  1:527680/1‑150 (MQ=255)                                                                                                                                               ttGGTGCAAAAA‑AAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGt       <  1:113110/148‑1 (MQ=255)                                                                                                                                                aggggagAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCATATTTCATTTTACCCAGAATTGCACCACCCTTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAGCGCCAg        <  2:318022/140‑1 (MQ=255)                                                                                                                                                   tGCTAAAA‑CAAACAGGCCAGCCAGTAAAGCAATTTTTTTCATCAACATTACTTATTACATTTTACCCAGAATTGCACCACACGTACCGCCAATCATGGCACCTTTAATCGCTCCATAGATTCTGTTGCCGGTCCGAACCCCAGTGACAg  >  2:360372/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAA‑CAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAG  >  NC_000913/1222448‑1222741

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