Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,448,521 T→G 11.2% F763C (TTT→TGT)  PP_3062 → tail length determination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,448,5210TG11.2% 320.9 / 7.4 161F763C (TTT→TGT) PP_3062tail length determination protein
Reads supporting (aligned to +/- strand):  ref base T (92/45);  new base G (21/2);  total (113/48)
Fisher's exact test for biased strand distribution p-value = 2.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

GAGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTC  >  NC_002947/3448373‑3448669
                                                                                                                                                    |                                                                                                                                                    
   ggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggg                                                                                                                                                  <  2:1265296/150‑1 (MQ=255)
   ggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggg                                                                                                                                                  <  1:6392204/150‑1 (MQ=255)
    gCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgggg                                                                                                                                                 <  2:2236722/150‑1 (MQ=255)
    ctgccggcagCGCCGCCGTAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGgggggg                                                                                                                                                >  1:3155468/5‑150 (MQ=255)
    ctgccggcagCGCCGCCGTAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTggggg                                                                                                                                                >  1:1176366/5‑150 (MQ=255)
     cTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggggg                                                                                                                                                <  2:5382822/150‑1 (MQ=255)
     cTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTggggg                                                                                                                                                >  2:4661380/1‑150 (MQ=255)
       aCGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGt                                                                                                                                              <  2:290395/150‑1 (MQ=255)
          gcggcgCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGc                                                                                                                                           >  2:3042913/1‑150 (MQ=255)
           cggcgCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCt                                                                                                                                          <  2:1998250/150‑1 (MQ=255)
             gcgcTGCCGGCCGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGGGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTgg                                                                                                                                        <  2:5756389/150‑1 (MQ=255)
             gcgcTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCGCCGTGGTGCCGGTGATCGGGACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTgg                                                                                                                                        <  1:5134163/150‑1 (MQ=255)
             gcgcTGCCGGCAGCCGGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATGGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTgg                                                                                                                                        <  1:5134162/150‑1 (MQ=255)
               gcTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGc                                                                                                                                      <  2:2315031/150‑1 (MQ=255)
               gcTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGc                                                                                                                                      <  2:5713818/150‑1 (MQ=255)
                    cGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCg                                                                                                                                 >  1:2686299/1‑150 (MQ=255)
                     ggCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGTGTCGTGTGGGGGGTGGCTGGGCAAGCGa                                                                                                                                >  2:4195596/1‑150 (MQ=255)
                     ggCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCGGGGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGGGGCTGGGCAAGCGa                                                                                                                                <  1:1814915/150‑1 (MQ=255)
                     ggCAGCCTGGCGGGCACCCGTGCCGGCGCTGCAGCGGGGGCGGCAATTCGCTCCGTGGTGCCGGTGAACGGAACCGCGGTGGGTGGCGCGGTCGGCGCTGTCCTCGGCGGCATGGGGGGCGAGTCGTTCGGGGGGTGGCTGGGCAAGCGa                                                                                                                                <  1:6173259/150‑1 (MQ=255)
                      gCAGCCTGGCGGGCCCCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGAt                                                                                                                               <  2:6856181/150‑1 (MQ=255)
                      gCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGAt                                                                                                                               >  1:6019426/1‑150 (MQ=255)
                       cAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATGGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATg                                                                                                                              <  2:4463117/150‑1 (MQ=255)
                       cAGCCGGGCGGGCACCCGTGCCGGTGCTGCAGCGGGGGCGGCGAGTGCCCCCGGGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGGGGCATGGGGGGCGCGTCGTTTGGGGGGTGGCTGGGCAAGCGATg                                                                                                                              <  2:2221692/150‑1 (MQ=255)
                          ccTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGtt                                                                                                                           >  1:2748709/1‑150 (MQ=255)
                           cgggCGGGCCCCCCGGCCGGGGCGGCCGCGGGGGCGGCCCGTGGCACCGCGGGGCCGGTGACCGGCCCCGCGGTGGGCGGCGCGGCCGGCGCCGTCCCGGGGGGCATGGGGGGCGCGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGttt                                                                                                                          <  2:2748709/148‑1 (MQ=255)
                           cTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATGGGCGCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTGGGGGGCATGGGGGGCGAGCCGTTTGGGGGGTGGCTGGGCAAGCGATGGttt                                                                                                                          <  2:295944/150‑1 (MQ=255)
                           cTGGCGGGCACCCGGGCCGGTGCTGCCGCGGGGGCGGCTATGGGCTCCGTGGGGCCGGCGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGATTGGGGGGTGGCTGGGCAAGCGATGGttt                                                                                                                          <  2:5816872/150‑1 (MQ=255)
                            tGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGTTTg                                                                                                                         >  1:2153496/1‑150 (MQ=255)
                            tGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTg                                                                                                                         >  1:4261495/1‑150 (MQ=255)
                            tGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTg                                                                                                                         >  2:201120/1‑150 (MQ=255)
                            tGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTg                                                                                                                         >  2:2069024/1‑150 (MQ=255)
                            gggCGGGCACCCTGGCCGGGGCTGCAGCGGGGGCGGCGATTGGCGCCGGGGGGCCGGGGATCGGTACCGCTGTGGGGGGCGCGGTCGGCGCTGGCCTCGGGGCCATGGGGGGCGAGACGCTTGGGGGGTGGCTGGGCAAGCGATGGTTTg                                                                                                                         <  2:5440321/149‑1 (MQ=255)
                                gggCACCCTCGCCGGTGCTGCAGCGGGGGCGGCTCTTGGCGCCGGGGTGCCGGTGCTCGGGCCCGCGGGGGGTGGCGCGGTCGGCGCGGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGa                                                                                                                     <  2:1032517/150‑1 (MQ=255)
                                gggCACCCGCGCCGGCGCGGCCGCGGGGGCCGCTCTGGGCGCCGTGGGGCCGGCGCCCGGGACCGCCGGGGGCGGCGCGCGCGGCCCGGCCCCCGGGGGCCGGGGGGGCGAGCCCTTAGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGa                                                                                                                     <  2:5731305/150‑1 (MQ=255)
                                 ggCACCCGTGCCGGCGCCGCAGCGGGGGCGGCAATGGGCCCCGTGGTGCCGGCGACCGGCACCGCGGCGGGTGGCGCGGCCGGCGCCGCCCTGGGTGGCATGGGGGGCGAGTCGCTGGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGAc                                                                                                                    <  2:1201460/150‑1 (MQ=255)
                                    aCCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCACGCGATGGTTTGGTGACGAg                                                                                                                 >  2:5330416/1‑150 (MQ=255)
                                    aCCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGGTTGGTGACGAg                                                                                                                 >  2:4289401/1‑150 (MQ=255)
                                    aCCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGGTTGGGGACGAg                                                                                                                 >  1:2595964/1‑150 (MQ=255)
                                     cccTTGCCGGTGCTGCAGCGGGGGCGGCTAGTGGCGCCGTGGGGCCGGTGATCGGTCCCGCGGTGGGTGGCGCGGTCGGCGCTGCCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgc                                                                                                                <  2:3352825/150‑1 (MQ=255)
                                      ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgcc                                                                                                               >  1:5029880/1‑150 (MQ=255)
                                      ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgcc                                                                                                               >  1:3892103/1‑150 (MQ=255)
                                      ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGGTTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAgcc                                                                                                               >  1:4197358/1‑150 (MQ=255)
                                       cTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCACGCGAGGGTTTGGTGACGAgccg                                                                                                              >  2:3647219/1‑150 (MQ=255)
                                       cTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgccg                                                                                                              >  1:4642260/1‑150 (MQ=255)
                                           ccGGTGCTGCAGCGGGGGCGGCTATAGGCTCCGTGGTGCCGGTGATCGGAACCGCAGTGGGTGGCGCGGACGCCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGATTGGGGGGTGGCTGGGCAAGCGAAGGATTGGTGACGAgccggccg                                                                                                          >  2:6246449/1‑150 (MQ=255)
                                            cGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGGTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCcga                                                                                                         >  2:716768/1‑150 (MQ=255)
                                              gTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCcgacg                                                                                                       >  2:2691261/1‑150 (MQ=255)
                                              gTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCcgacg                                                                                                       >  2:6796557/1‑150 (MQ=255)
                                               ggctgcCGCGGGGGCGCCCCTGGCCCCCGGGGGGCCGGCGACCGCCACCGCGGGGGGGGCCGCGGTCCGCGCGGCCCTCGGTGGCCTGGGGGGCGAGGCGATTGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCcgacga                                                                                                      <  2:6042714/149‑1 (MQ=255)
                                                 ctgcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGCTGGTTGGGTGACGAGCCGGCCGACGATc                                                                                                    >  2:1083187/1‑150 (MQ=255)
                                                  tgcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCa                                                                                                   >  2:2961223/1‑150 (MQ=255)
                                                  tgcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCa                                                                                                   >  2:3770736/1‑150 (MQ=255)
                                                   gcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCaa                                                                                                  >  1:4604035/1‑150 (MQ=255)
                                                    cAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAag                                                                                                 >  1:6010836/1‑150 (MQ=255)
                                                         ggggCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAagccgag                                                                                            >  1:1561936/1‑150 (MQ=255)
                                                            gcggcTCTTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAagccgagccg                                                                                         >  2:687651/1‑150 (MQ=255)
                                                             cggcTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCCGCCGACGATCAAGCCGAGCCgg                                                                                        >  1:2416293/1‑150 (MQ=255)
                                                                cTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAAt                                                                                     >  1:566345/1‑150 (MQ=255)
                                                                 tATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATc                                                                                    >  2:484938/1‑150 (MQ=255)
                                                                 ccccggCCCCGGGGGGCCGGCGACCGGCCCCGCTGCGGGGGGCGCGGGCGGCGCGGCCCGGGGGGGCAGGGGGGGCGCGCCGCGCGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATc                                                                                    <  2:2563523/146‑1 (MQ=255)
                                                                  aTTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGGGTGGGGGGTGGCTGGGCCAGCGCTGGTTGGGGGTCGAGCCGGCCGACGATCCAGCCGAGCCGGGATcg                                                                                   >  1:5470487/1‑150 (MQ=255)
                                                                        tCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGccgccg                                                                             >  1:4803296/1‑150 (MQ=255)
                                                                                  cGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTa                                                                   >  2:1174381/1‑150 (MQ=255)
                                                                                  cGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTa                                                                   >  2:3295573/1‑150 (MQ=255)
                                                                                  cGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTa                                                                   >  2:5791984/1‑150 (MQ=255)
                                                                                   ggTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAAGCGCCGCCGGCACTGGGGaa                                                                  >  2:1534852/1‑150 (MQ=255)
                                                                                      gATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGc                                                               >  2:4201838/1‑150 (MQ=255)
                                                                                        tCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGt                                                             >  2:4497787/1‑150 (MQ=255)
                                                                                        tCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGGTTGGTGACGAGCCGGCCGACCATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGt                                                             >  2:6001928/1‑150 (MQ=255)
                                                                                         cGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGtt                                                            >  1:3410149/1‑150 (MQ=255)
                                                                                          ggTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTa                                                           >  1:3047617/1‑150 (MQ=255)
                                                                                           gTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATAGCCGCCGGCACAGGGTAAAGCGTTAc                                                          >  1:190908/1‑150 (MQ=255)
                                                                                            tACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAcg                                                         >  2:680067/1‑150 (MQ=255)
                                                                                              ccGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGAGCAAGCCGAGCCGGAATCGCCGCCGGCCCGGGGGAAAGCGTTAcgcg                                                       >  2:5421505/1‑150 (MQ=255)
                                                                                               cGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGt                                                      >  2:5220672/1‑150 (MQ=255)
                                                                                               cGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCCGCCGACGATCAAGCCGAGCCGGAATCGCCGCCCGCACGGGGGAAAGCGTTACGCGt                                                      >  1:1749966/1‑150 (MQ=255)
                                                                                                  tgtgGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAAc                                                   >  1:546039/1‑150 (MQ=255)
                                                                                                    gccgcgccgcgGACGGCGCCGTCCCCGGCGGCCCGGGGGCCGAGCCGTTTGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACgg                                                 <  2:6157765/143‑1 (MQ=255)
                                                                                                       gtggCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGa                                              >  1:5046612/1‑150 (MQ=255)
                                                                                                        tggCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGTTCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGAt                                             >  2:4604874/1‑150 (MQ=255)
                                                                                                            gcgGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACGGGGTAAAGCGTTACGCGTAACGGTGATgcca                                         >  1:5119807/1‑150 (MQ=255)
                                                                                                            gcgGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGGCGTTGGGGGGGTGGCGGGGCACGCGCGGGGTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGACTCGCCGCCGGCCCGGGGGAAAGCGGTACGCGTAACGGTGATgcca                                         >  1:958516/1‑150 (MQ=255)
                                                                                                                  ggCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGCAGGTTTGGTGACGAGCCGGCCGACGAACAAGCCCAGCCGGAATCGCCGCCCGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAg                                   >  2:3637651/1‑150 (MQ=255)
                                                                                                                      cTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAg                               >  1:3800461/1‑150 (MQ=255)
                                                                                                                       tGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCAGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAATCCGAGCCGGAATCGACGCCGGCACAGGCTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGaatt                              >  2:6553822/1‑148 (MQ=255)
                                                                                                                       tGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCCGCCCACGATCAAGCCGAGCCGGAATCGCCGCCGGCCCTGGGTAAAGCGTTACGCGTAACGGTGAAGCCAGCCAAGGAAgg                              >  2:1491426/1‑149 (MQ=255)
                                                                                                                         tCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATTGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCg                            >  2:5352029/1‑150 (MQ=255)
                                                                                                                          ccccggcggCACGGGGGGCGAGGCGTTGGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGa                           <  2:5641606/143‑1 (MQ=255)
                                                                                                                            ttGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGc                         >  1:6859864/1‑150 (MQ=255)
                                                                                                                            tgcgcggCCTGCGGGGCGAGGCGCTCGGGGGGCGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGc                         <  2:6235750/145‑1 (MQ=255)
                                                                                                                               gtggCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCCAGCGAGGGTGTGGTGACGAGCCGGCCGACGCGCAAGCCGAGCCGGAATCGCCGCCGGCACGGGGTCAAGCGGTACGCGGAACGGTGATGCCAGCCAAGGAAGTCGAGccgg                      >  2:907552/1‑148 (MQ=255)
                                                                                                                                 ggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGATCTGGCCGACGATCAAGCCCAGCCGGAATCCCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCCACGAAGTCGAGCCagag                    >  1:6537600/1‑150 (MQ=255)
                                                                                                                                    aTGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACCAGCCGGCCGACGATCAAGCCCAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCAGAGccc                 >  1:4767666/1‑150 (MQ=255)
                                                                                                                                    aTGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGGTTGGTGACGAGCCGGCCGAAGATCAAGCCGAGACCGAATAGCCGCCGGCACTTGGTAAAGCGTTACGCGTAACGGTGATGCCACCCAAGGAAGTCGATACTGAGccc                 >  2:1790144/1‑150 (MQ=255)
                                                                                                                                       gggggCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGGAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGcc              >  2:5938253/1‑150 (MQ=255)
                                                                                                                                           gCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACgg          >  1:3800899/1‑150 (MQ=255)
                                                                                                                                           gCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACgg          >  1:5417793/1‑150 (MQ=255)
                                                                                                                                           gCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCAGAGCCCGCCACgg          >  2:2424547/1‑150 (MQ=255)
                                                                                                                                                  gTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGt   >  1:2677596/1‑150 (MQ=255)
                                                                                                                                                   tttGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTc  >  2:3251833/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTC  >  NC_002947/3448373‑3448669

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