Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,065,580	0	T	A	33.3%	46.2 / 10.6	34	T63S (ACC→TCC)	yihO	putative sulfoquinovose transporter
Reads supporting (aligned to +/- strand):  ref base T (16/6);  new base A (1/10);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 7.90e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGC  >  NC_000913/4065358‑4065814                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              aCGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAAt                                                                                                                                                                                                                  <  1:100303/250‑1 (MQ=255)      cAGTGCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑Gcc                                                                                                                                                                                                           <  2:341965/251‑1 (MQ=255)              caGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTa                                                                                                                                                                                                     >  2:131826/1‑249 (MQ=255)                        aaaaGTCCAAACATAAAAAACAGCGGGGGGGCAAGCGTCGTTTTAACCGGCAAACAAAAGGTGGGGGCGATAAACGGGAGCGTTGAAATTAACGCCGCGGAAACCGCTGCATATAAAAAAAAGGGACAAAATTTGCCTTGAGGCCAAA‑TTTTTTCCCG‑CGAGGCGAGTAAAAATCGGGG‑GAGCATATGGGAAAACGCGTAAAAAAA‑CTCCG‑CCACCGGAAAGATAATCCC‑GCCATAGTAAGCAGGCATccc                                                                                                                                                                                         <  1:221964/250‑1 (MQ=255)                                                         aTCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑cca                                                                                                                                                                                                                               >  2:206697/1‑181 (MQ=255)                                                         aTCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑cca                                                                                                                                                                                                                               <  1:206697/181‑1 (MQ=255)                                                                ttttAACCGGAAAACAAAGGGTGGTGGAGAAAAACGGGGGGTTTGAAATTAACGCCCCCGAACCCGCTAAATAAAAAAAAAAGGGACAAAATTTCCCTTAGGGCCAAA‑TATTTTTCCG‑CAGGTCGAAAAAAAACCCGGT‑GAGCATATCGGTAAACGGGGAAAAAAC‑CTCCC‑CCACAAAAAAAATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc                                                                                                                                                  <  2:137363/250‑1 (MQ=255)                                                                                     ggggggcaaaaatgtGGAGCGTTGCAAGTGACGCCGCGGGAACCGGTGAATAAAAAAAAAAGGACCAAAATTTCCCTTT‑GGCCCAAATATTTTCCCG‑CGAGTCGA‑AAAAAAACCGGTTGAGCATAGCGGGAAACGCGG‑AAAAAA‑CCTCGGCGCCCAAAAAAATAATCCC‑GCCATAGTAAGCAGGTATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGaa                                                                                                                            <  1:379722/237‑1 (MQ=255)                                                                                              tAAACTGAAGAGTTGCAATTAACGCCGCCGCCACCGCTGCGTATAAACTAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGCGTCGAGTAAAAATCCGCT‑GAGCATATCGGTAAACGCGGTCAAAAA‑CTTCG‑CCAACAACAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACCGAGTACCGATGCAGAGCATTAAATTAGAACCGACgacg                                                                                                                   >  2:66505/1‑248 (MQ=255)                                                                                               aaaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGc                                                                                                                  >  2:533042/1‑251 (MQ=255)                                                                                                        ggTTGAAATTAACGCCGCCGAACCCGCTGCAAATAAAAAAAAGGGACGAAATTTGCTTTTAGCCCCAT‑TTTTTTTCCG‑CGAGTCGAGAAAAAACCCGGT‑GAGCATATCGGAAAACGCGGAAAAAAC‑CTCCG‑CCACAAGAAAAATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccataccat                                                                                                         <  1:11150/250‑1 (MQ=255)                                                                                                         gTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccataccat                                                                                                         >  2:122661/1‑250 (MQ=255)                                                                                                           tGCAATTAACGCCGCCGGACCCGCTGCATCTAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGGGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAgg                                                                                                      >  1:203708/1‑251 (MQ=255)                                                                                                               aTTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCCTCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGAt                                                                                                  >  2:173541/1‑251 (MQ=255)                                                                                                                  aaCGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGCCGTCGCCCATACCATAGGCGAtttt                                                                                               >  1:536393/1‑251 (MQ=255)                                                                                                                        gccgGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCg                                                                                         >  2:152013/1‑251 (MQ=255)                                                                                                                                             ataaaGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAt                                                                                                                                                                      >  2:122496/1‑153 (MQ=255)                                                                                                                                             aaaaaGGGACAAAATTTGCCTTTAGGCCAAA‑TATTTTCCCG‑CGAGTGGAGTAAAAACCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAt                                                                                                                                                                      <  1:122496/151‑1 (MQ=255)                                                                                                                                               aaaGGGACGAAATTTGCCTTTAGGCCCAC‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAACGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTg                                                                  >  2:574770/1‑251 (MQ=255)                                                                                                                                                   ggACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGAtt                                                              >  2:189934/1‑251 (MQ=255)                                                                                                                                                   ggACAAAATTTGCCTTTATGCCAAA‑TATTGTCCCG‑GGAGTCAAGTAAAAATCGGGG‑AAGAATATCTTTAACCGCGGTAAAAAC‑CTTCC‑CCCCCTGAAA‑AAAATCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCTATGCAGAGCATTAAATTCGAACCGACGGCGCCCATACCATAGGCGATTTGTTCCCGCAGATGCAGTTTTAAGGTCAGTGGAtt                                                              <  2:569388/251‑1 (MQ=255)                                                                                                                                                                   tagcccaaa‑TATTTTCCCG‑CGAGTCGAATAAAAACCCGGT‑GAGCATATCGGAAAACGCGGAAAAAAC‑CTCCG‑CCACCAGAAAAATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAAt                                               <  1:268828/243‑1 (MQ=255)                                                                                                                                                                            atttttgcccG‑CGAGGCGATTAAAAATCCGGT‑TAGCATATCGTTAAACG‑GGAAAAAAATCTGCG‑CAACAAAAAAGAAAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTGTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACGCCTta                                      <  2:282869/243‑1 (MQ=255)                                                                                                                                                                               tttttCCGGCGAGTAGATGAAAAATCCGGG‑GAGCATATGGGAAAACGAGGAAAAAAA‑CTCCC‑CCACAAAAAAGAAAAACCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCTGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTtat                                     <  2:439743/250‑1 (MQ=255)                                                                                                                                                                               tttttCCG‑CGAGTCGAGAAAAAATCCGGG‑AAGCATATCGGAAAACGCGGAAAAAAA‑CTTCG‑CCCCAAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTtata                                    <  1:131826/250‑1 (MQ=255)                                                                                                                                                                                    ccG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGTTAACCGCGGAAAAAAC‑CTTCG‑CCACCAGAAAGATAACCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGccc                              <  1:122661/251‑1 (MQ=255)                                                                                                                                                                                    ccG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGACCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGTTTATGGTCAGACATAATCACTCCTTATAACGccc                              >  2:561453/1‑251 (MQ=255)                                                                                                                                                                                             cgagaAAAAATCCGGT‑GAGCATATCGGTAACCGCGGAAAAAAA‑CTCCG‑CCACAAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGc                      <  1:127656/251‑1 (MQ=255)                                                                                                                                                                                             cgagaAAAAATCCGGG‑GAGCATATCGGAAAACCCGGTAAAAAA‑CTCCC‑CCACAAGAAAAATAACCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAggg                       <  1:390585/250‑1 (MQ=255)                                                                                                                                                                                                      atgcGGT‑GAGCATATAGATAAACGCGGTAAAACA‑CTTCG‑CCACCAGAAAGATAATACC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAACATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGAGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat             >  2:35943/4‑251 (MQ=255)                                                                                                                                                                                                          ggT‑GAGCATATCGGTAAACGCGGAAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCAATAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGACGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAAAGCCCCGCAGCGCGACGTttattatt         >  1:74515/1‑251 (MQ=255)                                                                                                                                                                                                                     cccgatcccgcgGAAAAAAA‑CTCCC‑CCACCAGAAAAATAACCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGc  <  1:480754/240‑1 (MQ=255)                                                                                                                                                                                                                             cgcgGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGc  >  2:480754/1‑240 (MQ=255)                                                                                                                                                                                                                              gcgGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTAt                                                            >  1:166128/1‑181 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              ACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAA‑TATTTTTCCG‑CGAGTCGAGTAAAAATCCGGT‑GAGCATATCGGTAAACGCGGTAAAAAA‑CTTCG‑CCACCAGAAAGATAATCCC‑GCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGC  >  NC_000913/4065358‑4065814

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