Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 825,248 G→C 20.2% G287G (GGC→GGG)  ybhR ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913825,2480GC20.2% 111.6 / 3.4 62G287G (GGC→GGG) ybhRputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base C (11/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 2.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACC  >  NC_000913/825014‑825482
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
aCGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGGGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCACCTTCAAACCAAA‑AACGccca                                                                                                                                                                                                                                           >  2:442325/1‑237 (MQ=255)
       aaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACGATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCCTAAAGAAAAA‑CAACGCGATaaaag                                                                                                                                                                                                                                     >  2:1170016/1‑236 (MQ=255)
       aaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAGAAAAA‑AACGCCGATAAACg                                                                                                                                                                                                                                     >  2:932817/1‑238 (MQ=255)
       aaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑AACGcccgaaaaag                                                                                                                                                                                                                                     >  2:616357/1‑230 (MQ=255)
       aaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGGTAATCCAAGTTCGG‑TTTTGCAGCCCTACCGGCCGGGTTTCCACCGGAGAA‑ACGTAACCGGAAAGGGGAAAGGCGGGCATCCTAAAGAAAAA‑AACGCCCAGaaaaa                                                                                                                                                                                                                                     >  1:919682/1‑236 (MQ=255)
        aCATCGGGTACCCCTCTGCCCCTGTCGTGGCCGTTACCACCAGTAGCGGCCACAACCTATCCCAAACAATACCCAGACTCGCATCCTTCAAATAAAT‑CTGTTTG‑GTAATGTCCGTAAAGG‑GGCAAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAAAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAA‑CACGCCGATAAACGcc                                                                                                                                                                                                                                   <  2:586059/239‑1 (MQ=255)
            cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCGTCATAAAGAAAAA‑AAACCCGATAAAAGCct                                                                                                                                                                                                                                  >  2:260705/1‑236 (MQ=255)
            cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCCTAAAGACCAA‑AACACCCGTAAACGCccgct                                                                                                                                                                                                                               >  2:170512/1‑239 (MQ=255)
            cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACCCCCATCAACGCct                                                                                                                                                                                                                                  >  2:148078/1‑236 (MQ=255)
            cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGAATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAGAGCAA‑AACCCCGATAAAAGccc                                                                                                                                                                                                                                  >  2:937052/1‑235 (MQ=255)
             gcgTACGCCGCTGACCCTTCCTTGGCCTTAACCCCAAGTAGCGGCCAAAAACTATCCCACAAAATATCCAGCCTCGCTTCCTTCAAATAAAC‑CTGTTGG‑GTAATGTCCGTAAAGT‑GGCAAATAGGGTTAACCCACTTAAGG‑TTTGGCAGCCATACCGGCATGTTTTCCACCGGAAAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCctgctg                                                                                                                                                                                                                              <  2:620549/239‑1 (MQ=255)
              cgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACCCCGGTAAACGCCTGCcgt                                                                                                                                                                                                                             >  2:1166454/1‑239 (MQ=255)
              cgTACGCCGCTGACCCTGTCGTGGCAGTTATCACCAGTAGCGGCCACAACCTATTCCACATACTATCCAGACTCGCATGCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGACTAGGGTTAATCCACGTCTGG‑TTTTCCAGCCATACCGGCTTGTTTTCCACAGGAGGC‑ACACAACCCGAAAGGTGAAGTGCGGGCCACCTACAGACACA‑CACAccctacaaacccccatcg                                                                                                                                                                                                                             >  1:988821/1‑223 (MQ=255)
               gTACCCCGCTGCCCCTGTGGTGGCCGTTATCACGAGTAGGGGCAAAAAATTTTTCCCCACAATACCAAGACTCCCATCCTTAAAAAAAAT‑CTGTTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAACCCACGTAAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAAAA‑ACGTAACCGGAAGGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCtgtt                                                                                                                                                                                                                            <  2:389140/239‑1 (MQ=255)
                           cccTTTGGTGCCCGTTTCCCCCAGTAGCGCCCACAAACTTTCCCCAACAATACCAGGCCCCGAATCTTTAAAAAAATT‑CTGCTTG‑GTAAGGCACGAAAAGG‑GGCAAATGGGGTTAACCCAGGTCAGG‑TTTTGCAGCCATCCCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACaga                                                                                                                                                                                                               <  1:1166454/239‑1 (MQ=255)
                                ccgtgccttttacccccATTACCGGCCACAAAATGTTCCACACAATACCCAGACCGGCACACTAAAAATAAAC‑CTGCTGG‑TTAATGTCCGAAAAGG‑GGCGAATAGGGTAAATCCCCGAAGGG‑TTTTGCAGCCATCCGGGCATGTTTTCCACCGGAGAA‑ACGTAACGGAAAGGGAGAAGGGGGGGCATCAAAAAAACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGtga‑t                                                                                                                                                                                                         <  1:670728/224‑1 (MQ=255)
                                    ggCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAACAAA‑CACGCCGATAAACGCCCGCTGGTGGG‑GTGAACaagagggttgaa                                                                                                                                                                                                      >  2:1247951/1‑228 (MQ=255)
                                     gCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCTTAAAGACAAA‑CACGCCGCTAAAAGCCCGCCGTTGTG‑TTTAACAGAGGGA‑TGAAAt                                                                                                                                                                                                    >  2:762030/1‑239 (MQ=255)
                                      ccGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGA‑GGCGAATAGGGTTAATCCACTTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGTGGTCATCATAAAGACAAA‑CACCTCGATTACAGCCTCGCGTTGTGTTTGAAAagagagtaaaaac                                                                                                                                                                                                     >  2:694743/1‑230 (MQ=255)
                                        tttATCCCCATTGCCGCCCACAATTTGTTCCACACAATACCAAGACT‑GCATCCCTAAAAAAAATCCTGC‑TG‑GGAATGTTCCTAAAATGGGGGAATAGGGTAAACCCCCGTCAGT‑TTTGGAAGCCATACCGGCATTTTTCCCACGGGAAAA‑ACTTAACGGAAAAGGAGATTGGGGGGAAACAAAAAGAAAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCaa                                                                                                                                                                                                 <  1:585749/238‑1 (MQ=255)
                                                        gCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGGAGCCTAAAAACAAA‑AACCCCGCAAAAACCCCGCTGTTGTG‑TTGAACAACATGA‑TGAAAACAAAAAAACGAAAaccacc                                                                                                                                                                                 >  1:1014227/1‑239 (MQ=255)
                                                        gCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑AACGCCGATAAACGCCCGCTGTTGTG‑TTTAACAGCGGGG‑TGAAAACCACAGACCGAATCccac                                                                                                                                                                                  >  2:459423/1‑238 (MQ=255)
                                                          cacaAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATCGGGTTAATCCACGTCAGG‑TTTTGCAGCCAAACCGGCATGGTTTCCACCGGAGAC‑AACTAACCGGAAAGGAGAATAGCGGGGATCGCGAAAACAAA‑CAACCCGAAAAAAACCACCTGTTGTG‑TGGCAAAGAGAGA‑TGAAAACACAAGACCGAAACCCCCCAg                                                                                                                                                                               >  1:833023/1‑239 (MQ=255)
                                                                         acaATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGGATCCTAAAGACAAA‑CACACCGATAAACGCCTGCTGTTGGT‑TTGAAAAAATGGA‑TGAAATCAAAAAACCGAATCCCCCCAGCGATAAACCATaaa                                                                                                                                                                 >  2:114765/1‑238 (MQ=255)
                                                                             tatCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACCCCGATAAACGCCGGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAAAAGACCGAATCCCACCACCGGTAAACCATAAATCAcc                                                                                                                                                            >  2:749689/1‑239 (MQ=255)
                                                                             tatCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATAACGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGAATCTTAAAGACAAA‑ACAGCCCAGAAACGCCACCTTGTGGG‑TTGTACAGAGGGT‑TGAAAACAAACAACCGAACCCCACCACCCGTAAAcacaaaaaaaca                                                                                                                                                            >  1:984217/1‑228 (MQ=255)
                                                                                 cAGACTCGCATCCTTCAAATAAAT‑CTGTTTG‑GTAATGTCCGTAAAGG‑GGCAAATGGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATAATAAAGAAAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCg                                                                                                                                                        <  2:856289/239‑1 (MQ=255)
                                                                                     gtccccatgcttaaataaat‑CTGTTTGTGTAGTGTCCGAAAGGT‑GGGGTAGGAGTTTAGTCCCGGTCAGT‑TTTGGCAGCCATCCCGGGA‑GTTTTCCCCCAGAAAACCCGTCACCGAAAAGGAGATTGGCGGGCATCAAAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTa                                                                                                                                                      <  1:1156007/226‑1 (MQ=255)
                                                                                         cATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGCA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCCTAAAAACAAA‑CAACACGCTAAAAGACAGCTGGTGTG‑TTTAACAAAAGGA‑GTAAAACCACAGAACGACTACCACCACAGAGAAAcacaaaataacaaacagaaaaaaa                                                                                                                                                >  1:1016939/1‑216 (MQ=255)
                                                                                         cATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGGACAA‑CAACCCGATAAACGCCCGCTGTTTGG‑TTGAACAGAGGGA‑ATAAATAAACAAACCCAAACCACCCCACGATAAAAAAAAAAACACCAACCTaaaaaaa                                                                                                                                                >  1:1188282/1‑235 (MQ=255)
                                                                                         cATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAAAACGAATCCCACCACCGATAAACCATAAATCACCATCTTaaaagag                                                                                                                                                >  2:812194/1‑235 (MQ=255)
                                                                                         cATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCCGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCCCCACCGAGAAACCAAAAATCACCCCCGTAAAgaaa                                                                                                                                                >  2:367521/1‑236 (MQ=255)
                                                                                                          tcttttttgttgtggcctaaaaGT‑GGCGAATAGGTTTAACCCACGGTAGT‑TTTTCACACCATCCCGGAATGTTTTCCACCGGAGAA‑CCGTAACCGGAAAGGAAAATGGCGGGAAACAAAAAGAAAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGAt                                                                                                                                 <  1:932817/220‑1 (MQ=255)
                                                                                                             cTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGGTAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAAACCGAAACCCACCAGCGCTAAACCCTAAATAACCATCGTAAAGTAGAAAAGAGCCAGCGactcg                                                                                                                              >  2:520897/1‑234 (MQ=255)
                                                                                                             cTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTTAACAGAGTGA‑TGAAATCAACAGACCCAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCgg                                                                                                                             >  2:126388/1‑239 (MQ=255)
                                                                                                                      tGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATTATAAAGACCAA‑CAAGCCGGTAAAACCACGCTGTTGGG‑GTTGACAAAGGGAGTGAAAT‑AACAAAACGAATTCCACCACCAGTAAACCAAAAAACACACCCGTCAACAAGAAAAACGCCGCCCAACCGGGGAAgggg                                                                                                                     >  1:754882/1‑239 (MQ=255)
                                                                                                                               atggggCGAA‑AGGGTTAAACCCCATTCGGTTTTTGCAGCCATACCGGCATTTTTTCCACCGGAAAA‑ACGTAACGGGAAAGGAGAATGGCGGGCATAAAAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATa                                                                                                             <  2:984217/235‑1 (MQ=255)
                                                                                                                                    gcgatgaggttaaacccccgagaGG‑TTTGGCAGCAATGCGGGAATGTTTTCCACGGGAAAA‑ACGTAACGGGAAAGGAAAAGGGCGGGCATAAAAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGcc                                                                                                        <  2:746306/216‑1 (MQ=255)
                                                                                                                                           gggTTATCCCACTTCAGA‑TTTGGAGCCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATAc                                                                                                 <  2:20950/238‑1 (MQ=255)
                                                                                                                                             gTTAACCCGCGTCAGG‑TTTTGCAGCCATCCCGGAATGTTTTCCACCGGAAAA‑ACGTAACCGGAAAGGAGAATGGCGGGAATCAAAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCa                                                                                               <  2:1014227/238‑1 (MQ=255)
                                                                                                                                                 ctccgccgaaGG‑TTTTGCAGCCATCCCGGCATGTTTCCCACGGGAAAA‑TCGTAACCGGAAAGGAGAATGGCGGTCATCAAAAAGAAAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGc                                                                                          <  2:1230393/230‑1 (MQ=255)
                                                                                                                                                 acccccGTCAGG‑TTTGCCACCCATCCCGGCTTTTTTCCCCCGGGAAAA‑CCGTAACCGGAAAGGATAATGGCGGGCATAAAAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGc                                                                                          <  2:1074978/234‑1 (MQ=255)
                                                                                                                                                 aTCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCCGCTGTTGTG‑GTGAACAGAGAGA‑TGAAATCAACAGACCGAATCCCACCAACGGTCAACCATAAATACCAATCGTAAAGAAGAACACCGACAAGGATCCGGGGAAGGGGAGTTggaacaacaaaaaacaaacacc                                                                                           >  1:407605/1‑218 (MQ=255)
                                                                                                                                                    ccggTCAGG‑TTTTGCAGCCATCCCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAg                                                                                       <  2:112457/236‑1 (MQ=255)
                                                                                                                                                      ggttaGG‑TTTGGCGGCCATCCCGGCATGTTTTCCACCGAAAAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATAAAAAAGCAAAC‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTAGATA‑CGCCCAGATACCAATCGCCAGca                                                                                     <  2:405762/235‑1 (MQ=255)
                                                                                                                                                      cGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGGTAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAGGGGGATTTGAAA‑GCCCCCGAAACCACTccccccca                                                                                     >  2:898513/1‑235 (MQ=255)
                                                                                                                                                      cGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAAACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAAAAGCGCCAACGAACCGGGGAAGGGGAATTGAAACCGCCCA‑AAACCAATCGCCAGcc                                                                                     >  2:434764/1‑238 (MQ=255)
                                                                                                                                                      cGTCAGG‑ATTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGACTGCTGTTGTG‑TTGAACAGAAGGA‑GGTAATCACAAGACCGAACCCCACCCGCGATTAACCCAAAATCACCACATAAAAGAAGAAAAGACACACCGATCCCGCCCAGGGGGCTTTATA‑AGCCCacaaaaaaaaacccaac                                                                                      >  1:20950/1‑222 (MQ=255)
                                                                                                                                                                                        cACGGGAAAA‑ACGAACCCGAAAAGGAGATTGGGGGGAAAAAAAAAGAAAAA‑CACGCCGATAAACGCCTGCTGTTGGG‑TTGAACAGAGTGA‑TGAAAGCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACg                                                    <  1:434764/239‑1 (MQ=255)
                                                                                                                                                                                        cACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGGGGCCTGGAAGGGGGCGACAATTaaa                                                     >  2:1094314/1‑237 (MQ=255)
                                                                                                                                                                                                        cccGAAAAGGAGAATGGCGGCAAACATAAAGACAAC‑CTCGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGc                                     <  2:964620/238‑1 (MQ=255)
                                                                                                                                                                                                        aCCGAAAAGGAGAATGGCGGGCATAATAAAGAAAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGc                                     <  2:527361/239‑1 (MQ=255)
                                                                                                                                                                                                          cgaagaggggtaTGGGGGGCATAATAAAGAAAAA‑CACGCCGATAAACGCCGTCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCCCCAGCGATAAACCATAAATTCCCATCGTAAAGTAGAACAGCGCCAGCTATCCGGGGAAGGGGATTTGATT‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCg                                   <  1:204276/228‑1 (MQ=255)
                                                                                                                                                                                                             aagggAGAGTGGCGGTAATAATAAGGCACAA‑CACCCCGTAAAACGACCGCTGTGGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGGGATAAACCATAAATCACCAGCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGAt                                 <  1:613985/235‑1 (MQ=255)
                                                                                                                                                                                                               aGGAGAATGGGGGGCAT‑AAAAAAACAAACCACCCCGATAAACGCCTGTTGCTGTG‑TTGAACAGGGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGaa                              <  1:749689/239‑1 (MQ=255)
                                                                                                                                                                                                               aGGAGAATGGCGGATATAAAAAAGAAAAA‑CACGCCGATAAACGCCTGGTGTTGTG‑TTGAACAGAGTGA‑TGAAAGCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCCCAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGaa                              <  2:545955/239‑1 (MQ=255)
                                                                                                                                                                                                                                catAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGTGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGGGGCTGCAATTAACGCCGGTACGGCTTTGCCGATGAAGGTCTGCCAGGGGGGGa             >  1:1030022/1‑239 (MQ=255)
                                                                                                                                                                                                                                    aaagcaaaccacGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCCAGGTGGTGAgcg          <  2:754882/230‑1 (MQ=255)
                                                                                                                                                                                                                                          aa‑CACGCCGATAAACGCCTGCTGTTGGG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CCCCAAGATCCCAACCCCCAGCACCATGGGCGCCTGGAAGGGCGCGCGATTTACACCCGGAACGGGCTTCCCGTGGAATATCGGCAAGGGGGGGAGCGGGGaaaa   >  1:1033805/1‑238 (MQ=255)
                                                                                                                                                                                                                                           accacGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCCAGGTGGTGAGCGGCGAAAc   <  2:1279645/237‑1 (MQ=255)
                                                                                                                                                                                                                                           accacGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTACCGATGAAGATCTGCCAGGTGGTGAGCGGCGAAAc   <  2:749693/237‑1 (MQ=255)
                                                                                                                                                                                                                                           a‑cacGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCCCAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCCAGGTGGTGAGCGGGGAAAcc  >  2:844732/1‑239 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
ACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTG‑GTAATGTCCGTAAAGT‑GGCGAATAGGGTTAATCCACGTCAGG‑TTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAA‑ACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAA‑CACGCCGATAAACGCCTGCTGTTGTG‑TTGAACAGAGTGA‑TGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATA‑CGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACC  >  NC_000913/825014‑825482

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