Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	32.3%	48.1 / 12.5	31	intergenic (+31/‑72)	yfbL/yfbM	putative peptidase YfbL/DUF1877 domain‑containing protein YfbM
Reads supporting (aligned to +/- strand):  ref base G (3/18);  new base A (10/0);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 6.45e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTATTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTG  >  NC_000913/2386637‑2387105                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            tttATTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAAATTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCGGGAAAAAACGACATgaccttttttttttctcc                                                                                                                                                                                                                                     >  2:564591/1‑234 (MQ=255)    attccggtttttgttttCTTTGCCCACCTTAATTATGAGAAAAACGATTTCCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCAt                                                                                                                                                                                                                                 <  1:457800/236‑1 (MQ=255)           gTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAAAAAAAAAGACAAGCACAATTTTTTTTTccccctcttttttt                                                                                                                                                                                                                          >  1:210422/1‑240 (MQ=255)                                                            cATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGAc                                                                                                                                                                                                                                                         <  2:330940/171‑1 (MQ=255)                                                             aTAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTgg                                                                                                                                                                  <  1:294121/251‑1 (MQ=255)                                                                  gATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTAtt                                                                                                                                                             <  1:538421/251‑1 (MQ=255)                                                                         gACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGa                                                                                                                                                      <  1:123443/251‑1 (MQ=255)                                                                              cGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTg                                                                                                                                                 <  1:322943/251‑1 (MQ=255)                                                                                         gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACATGCACA‑TTTTTTTTTTCCTCATCTTTTTTTAGAAGCGGAACAAGAAA‑ATTAACGGGGG‑GAAAATAAAAT‑GGGA‑TGA‑TTTGCTATTTTTTAGAAATTTACTCAGaaaa                                                                                                                                       >  1:38376/1‑250 (MQ=255)                                                                                                  aaTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTCCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGaaan                                                                                                                                                   <  2:272309/229‑2 (MQ=255)                                                                                                  aaTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTCCTCGAGGCGGGG‑CCCGAAAAAAAGACAAGCACA‑TTTTTTTTTTTCTCCTCTGTTTTTAGATGCGAAACAAGAAA‑ATCAACGGAGG‑AAAAAAAAA‑A‑GGGAATGA‑TTGGCTAttttttagaaaa                                                                                                                                                   >  1:272309/1‑221 (MQ=255)                                                                                                     aCCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACt                                                                                                                          <  2:142878/250‑1 (MQ=255)                                                                                                            gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGTTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGt                                                                                                                  <  1:120875/251‑1 (MQ=255)                                                                                                                 gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGACAAGCACA‑TTTTTTTTTTCCTCATCTGTTTTTAAAAGCGAAACATGAAAATTAAACG‑AGG‑AAAAAAAAA‑A‑GGGGAAGA‑ATGGGTATTTTTTAGAAATTGAATCAGAAAAAATTAA‑ACAACTGCTTGTGTCTACt                                                                                                             >  1:386695/1‑251 (MQ=255)                                                                                                                                gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCa                                                                                                                                                        <  1:260037/194‑1 (MQ=255)                                                                                                                                            aGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAAAGCCCA‑ATTTATTTTTCCCCATTTGCTTTAAGATGAGAAAAATGAAA‑ATTAACGGAGG‑AAAAATAAA‑AAGGGAA‑GA‑ATGGGCATTTTGCAGAAATTGAAACAGAAAAGATTAA‑CCAACTTCTTTAGGCTACT‑TAAAAAACTTTGATGGaaaaaaataa                                                                                  >  1:542995/1‑241 (MQ=255)                                                                                                                                                   ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGAGGGATAATATACATGACACGc                                                                          <  1:410828/251‑1 (MQ=255)                                                                                                                                                        aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAACGAAATGCACA‑TTTTTTTTTTCCTCATCTGTTTTTAAAAGCGGAAAATGAAA‑ATTAACGGGGG‑GAAAAAAAA‑A‑GGGGATGA‑TTGGCTATTTTTTAGAAAATTACTCCGAAAAAAATAA‑ACAAATGCTTTAGTTTACT‑GAAAAAACTTTTATGGAGTATATAAAAGACACGCtttt                                                                      >  1:152591/1‑249 (MQ=255)                                                                                                                                                               ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTAc                                                              <  2:210422/251‑1 (MQ=255)                                                                                                                                                                   ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGcc                                                           <  1:577676/250‑1 (MQ=255)                                                                                                                                                                        tataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑TCAAAAAAAACGAAAAGCACA‑ATTTATTTTTCCCCCTCTTTTTTTAAATTCGAAAAAAGAAA‑ATTAACAGGGG‑AAAAAAAAA‑T‑GGGAATGA‑TTGGGTATTTTTTAGAAAATTATCCAAAAAAAAATAA‑CAAACTGCTTTATTTTACT‑GAAAAACCTTTTATTGATAATATACATTAAACACTTTTATGATTACGCCgcccgg                                                     >  2:61960/1‑246 (MQ=255)                                                                                                                                                                                tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAACGACATGCAAA‑TTTTTTTTTTCCTCATCTGTTTTTAGATGCGGAACAAGAAA‑TTTAACGGGGG‑AAAAAAAAA‑T‑GGGGATGA‑TTTGGTATTTTTTAAAAAATGACTCAGAAAAAAATAA‑ACAACTGCTTGAGTCTACT‑TAAAAAACTTTTATGGATAAAATAAATGACACGCTTTCAGGATTAAGCCGACTGGGTATTGa                                              >  1:501979/1‑250 (MQ=255)                                                                                                                                                                                       caacaGTAAATAATAAAAATTTCTCGAGGCGAGT‑CCAGAAAAAGCGACATGCATA‑TTTTATTTTCCCTCATCTGTTTTAAGATTCGAGACATGAAA‑TTCAACGGAGC‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGATTTCAGGATTACGCCGACTGGATATTTATAAACGCt                                      <  2:344891/251‑1 (MQ=255)                                                                                                                                                                                              aaaTAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAACGACAAGCAAA‑TTTTATTTTTCCCCCTCTGTTTTTAAAAGCGAAACCTGAAA‑TTTTACGGAGGGAAAAAAAAA‑T‑GGGA‑TGAATTGG‑TTTTTTGTAGAAATTGACTCAGAAAAGATTAAACCACCT‑CTTTAGGCTTCTTAAAAAACCTTT‑ATGGATTATAAAAATGACACGCCTTCTAGATTTTGCCGACTGGATATTTTTAAACGCTGGGAttt                               >  1:406440/1‑251 (MQ=255)                                                                                                                                                                                                   ataaAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTac                          <  2:38376/251‑1 (MQ=255)                                                                                                                                                                                                      aaaaTTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACt                       <  1:564591/251‑1 (MQ=255)                                                                                                                                                                                                         aTTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTg                    <  1:61960/251‑1 (MQ=255)                                                                                                                                                                                                                gAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTTTTTTCCCCCCTCTGTTTTTAAAAGCGAAACATGAAA‑TTCAACGGAGG‑GAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAAATGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑TAAAAAACTTTGTTGGATAATATAAATGACACACTTTCAGGATTACGCCGACTGGATATTTGTAAAAGCTGGGATTTTTTACACTTTTGGCTGa              >  2:295745/1‑250 (MQ=255)                                                                                                                                                                                                                      ggT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTTAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCAc       <  1:157779/251‑1 (MQ=255)                                                                                                                                                                                                                           naGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGA‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACctct   >  2:274988/2‑251 (MQ=255)                                                                                                                                                                                                                            aGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTg  >  1:564532/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            TTTATTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCACA‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAA‑TTCAACGGAGG‑AAAAATAAA‑T‑GGGAATGA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAA‑CCAACTGCTTGAGTCTACT‑GAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTG  >  NC_000913/2386637‑2387105

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: ‑