Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA36.7% 28.2 / 21.7 30intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (1/18);  new base A (11/0);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 2.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACT  >  NC_000913/2386645‑2387084
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
cGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGGCAGAAAAAAAGACAAGT‑ACaatatattttt                                                                                                                                                                                                                >  2:115051/1‑238 (MQ=255)
             tctGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGCCAGAAAAAACGACAAGC‑ACACTTTATTTTTTCTCATCTG‑tttt                                                                                                                                                                                                  >  1:209052/1‑239 (MQ=255)
                           ttatgGCCAACATAATTTTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑Cg                                                                                                                                                                                         <  1:306084/233‑1 (MQ=255)
                                ggCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGCTACACAAAGAGAAAAGCAACAATTAATTTTcctccctct                                                                                                                                                                                                        >  2:1011039/1‑209 (MQ=255)
                                                                                                 ttGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACattttattt                                                                                                                                                                                                                  <  1:703968/140‑1 (MQ=255)
                                                                                                       ggTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACa                                                                                                                                                                                    <  1:410021/162‑1 (MQ=255)
                                                                                                           aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGGCAGAAAAAACGAAAAAC‑ACAATTTTTTTTTTCTCATCt                                                                                                                                                                                                        >  2:703968/1‑140 (MQ=255)
                                                                                                            tACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTcc                                                                                                                                                                                                               <  1:242913/132‑1 (MQ=255)
                                                                                                                    acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACg                                                                                                                                                                        <  1:806909/161‑1 (MQ=255)
                                                                                                                              ttATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAAAAAAGGC‑AATTTTTATTTTTCC‑CCTCTGTTTTTAAAATG‑‑GAAACAAGAAAATTTACAGAGGGAAAATAAAAAGGAAT‑GATTTGGGTATTTT‑TTGAAAATTAATCAGCAAAGAATAAACCACCACTTTAGGCTACTTCAAAAAcca                                                                              >  2:778829/1‑238 (MQ=255)
                                                                                                                                 tCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGt                                                                                             <  1:539483/221‑1 (MQ=255)
                                                                                                                                           gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGa                                                                                                                                                                                                                                   <  2:988849/82‑1 (MQ=255)
                                                                                                                                           gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAACGa                                                                                                                                                                                                                                   >  1:988849/1‑82 (MQ=255)
                                                                                                                                             tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGataat                                                                  <  1:1003525/235‑1 (MQ=255)
                                                                                                                                                 gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACattttat                                                                                                                                                                                                                    <  1:170459/90‑1 (MQ=255)
                                                                                                                                                 gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACattttat                                                                                                                                                                                                                    <  1:188635/90‑1 (MQ=255)
                                                                                                                                                 gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGata                                                                    <  1:79233/229‑1 (MQ=255)
                                                                                                                                                  tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                                                                                                       <  1:135385/194‑1 (MQ=255)
                                                                                                                                                       gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTTCCCAAAAAAAGACAAGC‑ACATTTTTTTTTTTCTCATCTT‑TTTTTAGATG‑TGAAAAATGAAAATTAAAAGAGGGAAAAAAAAAGGGAATGGA‑TTGG‑TATTTTTTAGAAATTGAATCAGAAAAAAATAAACAAAAGCTTTAGGTTAC‑TTGAAAAACATTCTTGGAAAAtataaaataac                                                        >  2:1003525/1‑228 (MQ=255)
                                                                                                                                                           gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCCGAAAAAACGACAAAC‑ACATTTTATTTTTCCTCAt                                                                                                                                                                                                          >  2:170459/1‑90 (MQ=255)
                                                                                                                                                           gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGTCCGAAAAAAAGAAATGC‑ACATTTTATTTTTCCTCCTCTG‑TTTTTAGATG‑CGAAACAAAAAAATTAAAGGGGGAGAAAAAAAATGGAAT‑GT‑TTTGGTCTTTTTTTGAAAATTACTCAGAAAAAATTAAACAACAGTCTTAGTTTACTTGACAAAACTTACTTGGGAAATATACAttaaccga                                                     >  2:79233/1‑228 (MQ=255)
                                                                                                                                                            tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGata                                                                    <  1:265248/218‑1 (MQ=255)
                                                                                                                                                                      ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTTCCGAAAAAAAGACAAGC‑ACATTTTTTTTTTTCCCATCTG‑TTTTTAGATGAGGAAACA‑GAAATTTAAATGAGGGAAAATAAATGTGGAT‑GG‑TTTGGTATTTTTTAGAAAAATAATTAGAAAAAAATAACCAAATGGTTGGGTGTTA‑TGGAAAAACATTTTTGGTTAATATAAAAGAAACACGTTTAGGGTTAAGcc                                      >  2:33792/1‑238 (MQ=255)
                                                                                                                                                                      ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGCCAGAAAAAAAGAAAAGA‑ACATTTTATTTTTTCTCATCTG‑TTTTTAGATG‑AGAAAAATGAAAATAAACGGAGGAAAAAAAAATGGGGAT‑GA‑TTGGGTATTTTTCAGAAAATGAAACAGAAAAAAATAAATAAATGCTTTGGTCTAC‑TGGAAAAACTTTGATGGATAATATAAATGa                                                          >  2:265248/1‑218 (MQ=255)
                                                                                                                                                                            taCAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGa                               <  1:33792/239‑1 (MQ=255)
                                                                                                                                                                               aacaGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTg                                                                           <  1:319016/192‑1 (MQ=255)
                                                                                                                                                                                    tAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGata                                                                    <  1:786007/194‑1 (MQ=255)
                                                                                                                                                                                      aataatAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATaaa                     <  1:825275/239‑1 (MQ=255)
                                                                                                                                                                                              aaaTTTCTCGAGGC‑GGGTCCCGAAAAAAAGAAATGC‑ACATTTTATTTTTCCTCCTCTG‑TTTTTAAATG‑CGAAAAATGAAATTTAACGGAGGAAAAAAAAATGGGAAT‑GA‑TTGGGTATTTTTTAGAAAATGAATCAGAAAAAAATAAACAACTTCTTGAGTCTTC‑TGCAAAAACTATTATGGATAATATACATGa                                                          >  2:786007/1‑194 (MQ=255)
                                                                                                                                                                                                            ctgttccgaaaaaaGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAAGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACt  <  2:209052/228‑1 (MQ=255)
                                                                                                                                                                                                                          gCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGata                                                                    <  1:802352/157‑1 (MQ=255)
                                                                                                                                                                                                                          gCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGata                                                                    >  2:802352/1‑157 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
CGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGC‑ACATTTTATTTTCCCTCATCTG‑TTTTAAGATG‑CGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAT‑GA‑TTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACT  >  NC_000913/2386645‑2387084

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