Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913345,3280CA37.1% 65.4 / 12.9 35intergenic (+337/‑76)yahL/yahMuncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (21/1);  new base A (0/13);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 9.48e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGAAAAAAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCA  >  NC_000913/345119‑345562
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
tgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAaga                                                                                                                                                                                                                >  1:959647/1‑238 (MQ=255)
               cGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                                                                                                                >  1:857363/1‑239 (MQ=255)
                aataaaaaacGCCCCAGTAGTGTTTTCGCAGCCAGGCGGTATAAAAAAAATGATCCCAGGGCTATGTTTGGTGATTGTGTTGGTTGACTTTAATATAAAAAAAATCATAATCCCATCTTGTTTGTGTGGGTTTTTTTCTTTATCTTCCCCACCCCCCTCAAGCCTAAAAA‑AAAA‑TATGCTGAAGGGCTCAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGa                                                                                                                                                                                               <  2:575806/235‑1 (MQ=255)
                                                      gttaaaaaaataTGATCTGATGCAGATGTTTGGTGATTGTGTGGTTTGAAATTCATATGAAAAAAATCAAAACCCCATCATGTTGGTGTGGGTTCTTTTCTGTATCTAACGCACCGCACCCAAGCTTGACAA‑AAAA‑TATGCTGTAGGGCTCATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCgg                                                                                                                                                          <  2:770511/232‑1 (MQ=255)
                                                                     tctcATGCAGATGTTTTGTGAATGTGTTGGTTGACCTTCATATGAAAAAAATCATAATTCCCTCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGGATGGTa                                                                                                                                                                                  >  1:379322/1‑199 (MQ=255)
                                                                          tGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttc                                                                                                                                     >  2:786727/1‑239 (MQ=255)
                                                                               aTGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCt                                                                                                                                                                                                                           >  1:940965/1‑148 (MQ=255)
                                                                                 tttttGTGATGTTGTTGGTGACCTTTCATATAAAAAATATCATATCCCCAAGTTTTTTGGGTGGGTTCTTTTCTGTTTCTTCCACACCAC‑CCCAATCCTAAAAGAAAAA‑TTGGTTGTAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                                                                                                              <  2:276016/238‑1 (MQ=255)
                                                                                             tgtTGGTTGACATTCATATAAAAAAAATCAAAATCCCATCATTTTTGTGGGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAA‑AAAA‑TATGCTGAAAGGCTCATATAAAAAACCGCCATTAGc                                                                                                                                                                                                                            <  1:6173/133‑1 (MQ=255)
                                                                                                tggttaTCATTCAAATGAAAACATTCATATTTCCATCATGTGTGTGTGGAGTTTTGTCTATATCTGCCACACCGCACTCAACCCTAACAA‑AAAA‑TTTGCTTAAAGCCAAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                                                                                                                <  2:959647/238‑1 (MQ=255)
                                                                                                   ttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcat                                                                                                             >  1:327610/1‑238 (MQ=255)
                                                                                                       cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGaa                                                                                                                                                                                                                  >  2:6173/1‑133 (MQ=255)
                                                                                                       cATTCATAAGAAAAAAATCATAATCCCATCATGTTTGTGTGGGTTTTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAA‑AAAA‑TATGCTGTAAGGCTCATAAAAAAAAc                                                                                                                                                                                                                                      <  2:751967/113‑1 (MQ=255)
                                                                                                       aaTTCAAAAGAAAAAATTCATAATTCAATCATGTTTGTGTGGGGTTTTTTCTGTCTCTTACCCATCGCACTCAACACAGAAAA‑AAAA‑TATGCTGTAGGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg                                                                                                        <  2:857363/238‑1 (MQ=255)
                                                                                                            atatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg                                                                                                   >  1:400892/1‑239 (MQ=255)
                                                                                                                 aaaaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGc                                                                                                                                                                                                                            >  1:751967/1‑113 (MQ=255)
                                                                                                                                    catTTTTGTGGGGGGTTTTTTCTGTATCTTACGCACCGCACTAAACCCTGAAAA‑AAAA‑TAGGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAagag                                                                                                                                                                                                               <  1:205299/107‑1 (MQ=255)
                                                                                                                                      tttttGGGTGGGTTTTTTCCTGAATCTCACGAACCGAAC‑CAAACCTGAAAA‑AAAATTATGGTGAAAGCCAAAAAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTAtt                                                                         <  2:327610/237‑1 (MQ=255)
                                                                                                                                      tGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAgg                                                                                                                                                                                                                    >  1:114391/1‑100 (MQ=255)
                                                                                                                                       tgttGTGTGTGTTCTTTTTTGTATCTTGCGCACCGCACTAACACCAAAAAG‑AAAA‑TGAGCTGAAGGTCTCAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTAtt                                                                         <  2:459044/236‑1 (MQ=255)
                                                                                                                                              tgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACaa                                                                                                                                                                                                     >  2:205299/1‑107 (MQ=255)
                                                                                                                                                 ggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGc                                                              >  1:112322/1‑239 (MQ=255)
                                                                                                                                                  gTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATgggg                                                               >  1:459044/1‑237 (MQ=255)
                                                                                                                                                            tATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGaa                                                   >  1:8717/1‑239 (MQ=255)
                                                                                                                                                                        cGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGc                                                                                                                                                                                                                            >  2:240352/1‑58 (MQ=255)
                                                                                                                                                                        cGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTAATATCAAAAACCGCCATTAGc                                                                                                                                                                                                                            <  1:240352/58‑1 (MQ=255)
                                                                                                                                                                                 cccTGAAAA‑AAAA‑TATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacg                                                                                                                                                                                                         <  1:834112/67‑1 (MQ=255)
                                                                                                                                                                                         aaaaaa‑TTTGCTGAAAGGCTCATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa                                                                                                      <  2:642922/158‑1 (MQ=255)
                                                                                                                                                                                            aaa‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGa                                                                                                                                                                                               >  2:834112/1‑68 (MQ=255)
                                                                                                                                                                                                    cTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACa                                                                                            >  1:642922/1‑160 (MQ=255)
                                                                                                                                                                                                      gTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGAAAAAAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCt           >  1:182016/1‑239 (MQ=255)
                                                                                                                                                                                                               atatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                                                                                                                     >  1:395483/1‑60 (MQ=255)
                                                                                                                                                                                                               atatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGa                                                           >  1:176826/1‑182 (MQ=255)
                                                                                                                                                                                                               atatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGAAAAAAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCa  >  1:806969/1‑239 (MQ=255)
                                                                                                                                                                                                               atatAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGa                                                           <  2:176826/182‑1 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
TGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAG‑AAAA‑TATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGAAAAAAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCA  >  NC_000913/345119‑345562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑