Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,972,154 T→A 6.5% intergenic (+32/+87) PP_1762 → / ← yciK hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,972,1540TA 6.5% 576.2 / 17.6 184intergenic (+32/+87)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (72/95);  new base A (5/7);  total (78/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

CGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGT  >  NC_002947/1972005‑1972302
                                                                                                                                                     |                                                                                                                                                    
ccaggcctaccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCtt                                                                                                                                                      >  2:4604132‑M1/133‑150 (MQ=255)
 caggcctaccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTg                                                                                                                                                     >  2:308174‑M1/132‑150 (MQ=255)
      gcgcAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGc                                                                                                                                                <  2:2431249/150‑1 (MQ=255)
      ctaccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGc                                                                                                                                                >  2:6605891‑M1/127‑150 (MQ=255)
        accccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                              >  2:4122210‑M1/125‑150 (MQ=255)
        accccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                              >  2:1837467‑M1/125‑150 (MQ=255)
        accccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                              >  2:1447669‑M1/125‑150 (MQ=255)
        accccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                              >  2:1253539‑M1/125‑150 (MQ=255)
         ccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaccgcccaataaaccggcatCGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                             >  2:5292721‑M1/124‑147 (MQ=255)
         ccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                             >  1:3886490‑M1/124‑147 (MQ=255)
         ccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtggacctgtacctgacggggccggacagcacaggggggaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCAAGCccgcgac                                                                                                                                             >  1:1003249‑M1/124‑147 (MQ=255)
         cAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGTGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCAAGCCCGcgaa                                                                                                                                             <  1:2244945/150‑4 (MQ=255)
         cAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCAAGCCCGcgaa                                                                                                                                             <  2:3681068/150‑4 (MQ=255)
          aGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCGTGGCCGCcccc                                                                                                                                            >  1:5972797/1‑150 (MQ=255)
           gggCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                           <  2:370867/150‑1 (MQ=255)
           ccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaatatacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                           <  1:556590‑M1/29‑1 (MQ=255)
            ggCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                          <  1:4052689/150‑1 (MQ=255)
            ggCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                          <  2:3246372/150‑1 (MQ=255)
            ggCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                          <  1:4630716/150‑1 (MQ=255)
            cggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCAagcccgcgaagag                                                                                                                                          >  2:1377333‑M1/121‑144 (MQ=255)
             gCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                         <  2:2954247/150‑1 (MQ=255)
             gCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                         <  2:6415230/150‑1 (MQ=255)
              gcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAg                                                                                                                                        >  2:5569719‑M1/119‑150 (MQ=255)
              gcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAg                                                                                                                                        >  2:1232506‑M1/119‑150 (MQ=255)
                   aacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                             <  2:3431605‑M1/27‑4 (MQ=255)
                     cAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCAcc                                                                                                                                 >  1:705810/1‑150 (MQ=255)
                     cAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCAcc                                                                                                                                 <  1:4858972/150‑1 (MQ=255)
                      ccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCg                                                                                                                                <  2:3593520‑M1/40‑1 (MQ=255)
                             tcGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCt                                                                                                                         <  2:4585200/149‑1 (MQ=255)
                               gACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCa                                                                                                                       >  2:6783063/1‑150 (MQ=255)
                                aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                      <  1:5807391/150‑1 (MQ=255)
                                aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                      <  2:2354199/150‑1 (MQ=255)
                                aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                      <  1:6806901/150‑1 (MQ=255)
                                aCTGCATGTTCCTGGTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                      <  1:1346804/150‑1 (MQ=255)
                                  tGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGGGTTCCCCCTCCCACACGATAACCGCAAACCAAAAgg                                                                                                                    >  1:444843/1‑150 (MQ=255)
                                    cATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGAt                                                                                                                  <  1:2534184/150‑1 (MQ=255)
                                      tgttcctgttccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTg                                                                                                                <  1:540137/150‑1 (MQ=255)
                                       gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTgc                                                                                                               <  2:5452100/150‑1 (MQ=255)
                                          cctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAATCCTTGAGGGTTgcggt                                                                                                            <  2:3366790/150‑3 (MQ=255)
                                           gaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                              <  1:976439‑M1/26‑1 (MQ=255)
                                            aagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGCCGCAACCCTCAAGGATTGCGCcct                                                                                                          >  1:743564‑M1/89‑150 (MQ=255)
                                              ttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCGTGCCCGCTCCCACACGCTCACCGCAAACCCCAAGGATTTCGCccggc                                                                                                        >  1:6237079/1‑150 (MQ=255)
                                                cctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgcct                                                                                                      <  2:88379/150‑1 (MQ=255)
                                                  tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                    >  1:614395‑M1/83‑150 (MQ=255)
                                                  tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                    >  1:4991206‑M1/83‑150 (MQ=255)
                                                   ggAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                   <  1:4557454/150‑1 (MQ=255)
                                                   ggAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                   <  1:4267518/150‑1 (MQ=255)
                                                   ggAGCGCTGCATGCAGCTGGAAGGGCGCCAACAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                   <  2:1032052/150‑1 (MQ=255)
                                                   catgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                   >  1:5547137‑M1/82‑150 (MQ=255)
                                                    gAGCGCTGCATGCAGCTGGAAGGGCGCCATCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                  <  1:5507682/150‑1 (MQ=255)
                                                    gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                  <  2:1600509/150‑1 (MQ=255)
                                                    gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                  <  2:2536545/150‑1 (MQ=255)
                                                    atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaatgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                  >  1:2855050‑M1/81‑150 (MQ=255)
                                                     aGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTggg                                                                                                 <  1:2300466/150‑1 (MQ=255)
                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                <  1:6017149/150‑1 (MQ=255)
                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                <  2:4530478/150‑1 (MQ=255)
                                                      gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                >  1:1052230‑M1/79‑150 (MQ=255)
                                                      gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                >  2:5869644‑M1/79‑150 (MQ=255)
                                                      gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                >  1:5229403‑M1/79‑150 (MQ=255)
                                                      gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                >  2:2634884‑M1/79‑150 (MQ=255)
                                                       ccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAg                                                                                               >  2:3018551‑M1/78‑150 (MQ=255)
                                                        gcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                              <  2:3742970/150‑1 (MQ=255)
                                                        cggtgtacctgtacctggtggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTTCGCCCTGCCTGTGGGAGc                                                                                              >  1:3453500‑M1/77‑150 (MQ=255)
                                                         cTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCg                                                                                             <  1:351825/150‑1 (MQ=255)
                                                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCgg                                                                                            <  1:5129382/150‑1 (MQ=17)
                                                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGcc                                                                                                            <  1:3274936/134‑1 (MQ=255)
                                                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGGGGGAGCgg                                                                                            <  2:258893/150‑1 (MQ=255)
                                                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTAGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                            <  2:6260008/150‑1 (MQ=255)
                                                          gtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                            >  1:1857941‑M1/75‑150 (MQ=255)
                                                          gtgtacctggacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                            >  2:3199229‑M1/75‑150 (MQ=255)
                                                           gcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCggg                                                                                           <  2:2439454/150‑1 (MQ=17)
                                                            catgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGc                                                                                          <  2:1322914/150‑1 (MQ=255)
                                                              tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCaa                                                                                        <  2:6201260/150‑1 (MQ=255)
                                                               gcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCAGCCTGGGGGAGCGGGCAAg                                                                                       >  2:5028213/1‑150 (MQ=255)
                                                                ctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCCAGCccgcgaa                                                                                                                                             <  1:305839‑M1/27‑4 (MQ=255)
                                                                ctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                             >  2:305839‑M1/69‑92 (MQ=255)
                                                                caGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGc                                                                                      <  1:1509727/150‑1 (MQ=255)
                                                                 aGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCGTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGCTTGCGCCCTGCCTGGGGGAGCGGGCCAGcc                                                                                     >  2:4694314/1‑150 (MQ=255)
                                                                   cTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGGCCCCACCCACACGGGCAACGCAAACCCCAAGGATTGGGCCCTGCCCGGGGGGGCGGGCacccccc                                                                                   >  2:2272451/1‑144 (MQ=255)
                                                                   cTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                        >  2:6038992/1‑129 (MQ=255)
                                                                    tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGTGCGGGCAAGCCcgc                                                                                  <  2:4831363/150‑1 (MQ=14)
                                                                    tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGTTTGCGCCCTGCCTGTGGGAGCGGGCAAGCCcgc                                                                                  <  2:541871/150‑1 (MQ=255)
                                                                           gcgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCGTCGCGGGGGTGGCCGCCCCCCCCCCAACACCCCCACCCCCAACGGTTGGGGCCCGCCGGGGGGGGCGGGGGAGCCCCCCCAGAgg                                                                           >  1:3772473/1‑150 (MQ=255)
                                                                            cgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGc                                                                                                                                                    >  2:705809/1‑76 (MQ=255)
                                                                                   aaGGCCTTATGCGGCCGGCGTTGGCTGCGGGCTGATCGATTTCCTGTACCGGCCTCTGCGCGGGCTGGCCGGCTCCCACAGGATCACCGCAACCCTCACGGATTGCGCCCTGCCTGTGGGCGCGGGCAAGCCCGCGAAGAGGCCGATGCa                                                                   <  2:3305760/150‑1 (MQ=38)
                                                                                          tATGCGGCGGGCTTTGGCTGGGGGCTTAGCGCTTTCGTGTCCCGCCCTCTTCGGGGGCTTGCCCGCACCCACTGGTCCACCCCAACCCAAATGGTTTGCGCCCTCCCTGTGGGCGGGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTAt                                                            <  2:213067/150‑1 (MQ=255)
                                                                                          tATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGTCCGCGAAGAGGCCGATGCAGGTTTAt                                                            >  2:3830821/1‑150 (MQ=255)
                                                                                          caggcggcgGGCGTGGGCTGCGGGCTGATCGCTGGCCTGTACCGGCCTCTCCGCGGGCTTGCCCGCCCCCACCGGATCACCGCACCCCTCACGGATTGCGCCCGGCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTAt                                                            <  2:6585890/147‑1 (MQ=37)
                                                                                                 gcgggGTGGGCGGGGGGTTTATCGCTTTCGTTCCCCGGCCTCGGCGGGGGCGGGCCCGCCCCCCCCGTATCCCCCAAACCCTCTGGGTTTGCGCCCGGCGGGTGGGCGGGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGt                                                     <  2:2852241/149‑1 (MQ=255)
                                                                                                                 gCTGATCGCTTGCCTGTCCGGGCCTCTGGGGGGGCTTGCCCGCCCCCACTGGATCACCGCAACCCCCACGCTTTGCGCCCGGCCTTTGGGTGGGGGAAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGtt                                     <  2:200031/150‑1 (MQ=255)
                                                                                                                     ccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                    <  1:656004‑M1/68‑1 (MQ=255)
                                                                                                                     ccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                    >  2:656004‑M1/16‑83 (MQ=255)
                                                                                                                      tCGCTTTCCTGTACCGGCCTCTTCGCGGGCGTGCCCCCGCCCACACGCTCACCCCAACCCACCAAGCTTTCCCCCCGCCCGTGGGGGCCGGCCAGCCCCCCCAAAGGCCCAGGCCGGTGTATTTGGCGGTGTACCCCAGCCCGGGCAccc                                >  2:4050535/1‑150 (MQ=255)
                                                                                                                                    CGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCAGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTTCGGc                  <  2:5672430‑M1/150‑1 (MQ=255)
                                                                                                                                       ccTCTTCGCGGGCTTGCCCCCACCCACAGGATCACCGCCACCCGCAAAGATTTCGCCCTGCCGGTGGGTGCGGGCCAGCCCGCGCCGAGGCCGGTGCAGGGTTGTTGGGCGGTTAACGCCGGCCCCTGCACCCCCGGGCGGTCCGGcccc               >  2:948353/1‑150 (MQ=255)
                                                                                                                                           TTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATca           >  2:407218‑M1/1‑150 (MQ=255)
                                                                                                                                                 GGCTGGCCCGCTCCCCCCGGATCACCGCAACCCTCACGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtaca     >  2:5142152‑M1/1‑150 (MQ=255)
                                                                                                                                                  gcgtGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag    <  2:5417943/147‑1 (MQ=255)
                                                                                                                                                  GCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag    <  2:2471619‑M1/150‑1 (MQ=255)
                                                                                                                                                   CTTGCACGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacagg   <  2:5090124‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGT  >  NC_002947/1972005‑1972302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: