Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,448,521 T→G 13.1% 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,5210TG13.1% 349.3 / 6.8 174F763C (TTT→TGT) PP_3062tail length determination protein
Reads supporting (aligned to +/- strand):  ref base T (86/55);  new base G (15/12);  total (101/73)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

CGAGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTC  >  NC_002947/3448372‑3448669
                                                                                                                                                     |                                                                                                                                                    
cGAGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGtt                                                                                                                                                      >  2:3277358/1‑150 (MQ=255)
cGAGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGtt                                                                                                                                                      <  1:5108533/150‑1 (MQ=255)
cGAGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCGATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGtt                                                                                                                                                      >  2:2829646/1‑150 (MQ=255)
  aGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTg                                                                                                                                                    <  1:72967/150‑1 (MQ=255)
   gggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgg                                                                                                                                                   <  1:4802729/150‑1 (MQ=255)
   gggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgg                                                                                                                                                   <  2:4048355/150‑1 (MQ=255)
    ggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggg                                                                                                                                                  <  2:288214/150‑1 (MQ=255)
    ggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggg                                                                                                                                                  <  1:641127/150‑1 (MQ=255)
    ggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggg                                                                                                                                                  <  1:171620/150‑1 (MQ=255)
    ggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggg                                                                                                                                                  <  1:1974203/150‑1 (MQ=255)
     gCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgggg                                                                                                                                                 <  1:1084371/150‑1 (MQ=255)
     gCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgggg                                                                                                                                                 <  1:5387219/150‑1 (MQ=255)
       tACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgggggg                                                                                                                                               >  2:334501/1‑150 (MQ=255)
       tACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgggggg                                                                                                                                               >  1:4329850/1‑150 (MQ=255)
        aCGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCTTGTGATGATATTTGCTCAGTGGTGCAGGTGATCGGTACAGCTTTGGGTGGCGCGGTCTGCGCTGTCCTTGGTGGCATGGGGGGCGATTCGTTTgttgtgt                                                                                                                                              <  2:177116/150‑7 (MQ=255)
              gcgcTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCTTTTGGGGGGTGGCTgg                                                                                                                                        <  2:3188137/150‑1 (MQ=255)
                gcTGCCGGCAGCCGGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGAGCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGc                                                                                                                                      <  1:2091966/150‑1 (MQ=255)
                  tGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCaa                                                                                                                                    >  2:5630411/1‑150 (MQ=255)
                  tGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCaa                                                                                                                                    >  1:3892427/1‑150 (MQ=255)
                      ggCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGa                                                                                                                                <  2:1865484/150‑1 (MQ=255)
                      ggCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGa                                                                                                                                >  1:3736393/1‑150 (MQ=255)
                       gCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGAt                                                                                                                               >  1:3223612/1‑150 (MQ=255)
                        cAGCCGGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGCTACCGCGGTGGGTGCCGCGGTCGGCGCTGCCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATg                                                                                                                              <  2:2388189/150‑1 (MQ=255)
                         aGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGCTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGAGGCTGGGCAAGCGATgg                                                                                                                             <  2:3707674/150‑1 (MQ=255)
                          gCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGCTCGGTACCGCTGTGGGTGGCGCGGCCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGGGGCTGGGCAAGCGATGGt                                                                                                                            <  2:612740/150‑1 (MQ=255)
                             tGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTGGGGGGGTGGCTGGGCAAGCGATGGTTTg                                                                                                                         >  2:1135750/1‑150 (MQ=255)
                             gggCGGGCCCCCGGGCCGGGGCGGCTGCGGGGGCGGCGCCGGGCGCCGTGGTGCCGGGGCCCGGGCCCGCGGTGGGGGGCGCGGTCGGCGCGGGCCGGGGGGGCCGGGGGGGCGCGTCGTGGGGGGGGGGGCTGGGCAAGCGATGGTTTg                                                                                                                         <  2:5173708/149‑1 (MQ=255)
                                cGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTg                                                                                                                      >  1:338715/1‑150 (MQ=255)
                                 gggCACCCTGGCCGGTGCTGCAGCGGGGGCGGCTATCGGCTCCGTGGGGCCGGTGATCGGTCCCGCTGTGGGTGGCGCGGTCGGCGCGGTCCTTGGTGGCATGGGGGGCGAGTCGTTGGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGa                                                                                                                     <  2:1387273/150‑1 (MQ=255)
                                    cACCCTTCCCGGTGCTGCAGCGGGGGCGGCTATGGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGCCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGa                                                                                                                  <  2:4939656/150‑1 (MQ=255)
                                    cACCCGTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGa                                                                                                                  <  2:699758/150‑1 (MQ=255)
                                     ccccgggCCGGTGCGGCCGCGGGGGCGGCGCGGGGCTCCGTGGTGCCGGTGCGCGGTCCCGCGGTGGGGGGCGCGGTCGGCGCGGTCCTTGGGGGCATGGGGGGCGCGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAg                                                                                                                 <  2:2786385/144‑1 (MQ=255)
                                       ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgcc                                                                                                               >  1:4791556/1‑150 (MQ=255)
                                       ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCCAGCGATGGGTTGGTGACGAgcc                                                                                                               >  2:2474728/1‑150 (MQ=255)
                                          tGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGATTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgccggc                                                                                                            <  2:4189059/150‑1 (MQ=255)
                                           cccGGCGCCGCAGCGCGGGCCCCTACTCGCCCCGTGGTGCCGGCGCCCGGGACCGCCCGGGGGGGCGCGGGCGGCGCGGCCCCGGCTGCCACGGGGGGCGAGACGATTGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAgccggcc                                                                                                           <  2:1889929/149‑1 (MQ=255)
                                             cGGTGCGGCAGCGGGGGCGGCTATCGCCCCCGTGGTGCCGGTGACCGGTACCGCTGTGGGTGGCGCGGCCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCcga                                                                                                         <  2:2008074/150‑1 (MQ=255)
                                               gTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCcgacg                                                                                                       >  1:5054767/1‑150 (MQ=255)
                                                   tgcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTTGTGACGAGCCGGCCGACGATCa                                                                                                   >  1:1676026/1‑150 (MQ=255)
                                                   tgcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCa                                                                                                   >  1:3018949/1‑150 (MQ=255)
                                                   tgcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCa                                                                                                   >  1:265989/1‑150 (MQ=255)
                                                    gccccggcggCGCCTATGCGCACCCTGGCCCCGGTGAGCGGTCCCGCCGGGGGGGGCGCGGGCGCCCCCGTCCTGGGTCCCACGGGGGGCGCGTCGCCGGGGCGCGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCaa                                                                                                  <  2:5410115/146‑1 (MQ=255)
                                                    gcAGCGGGGGCGGCTCTTGGCTCCGTGGTGCCGGTGCTCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGACCTTGGTGGCATGGGGGGCGAGTCGATTGTGGGGTGGCTGGGCAAGGGATGGGGTGGTGACGAGCCGGCCCACGCTCaa                                                                                                  >  2:4854364/1‑150 (MQ=255)
                                                    gcAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGGTGGTTTGGTGACGAGCCGGCCGACGATCaa                                                                                                  >  2:802826/1‑150 (MQ=255)
                                                     cAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGCTCGGTCCCGCTGGGGGTGGCGCGGTCGGCGCTGACCTGGGTGGCAGGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAag                                                                                                 <  2:1297645/150‑1 (MQ=255)
                                                      aGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAagc                                                                                                >  1:3632691/1‑150 (MQ=255)
                                                      aGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAagc                                                                                                >  2:1564389/1‑150 (MQ=255)
                                                      aGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGGGGGGGGTGGCGGGGCAAGCGCGGGGTTGGTGACGAGCCGGCCGACGATCAagc                                                                                                >  2:1785728/1‑150 (MQ=255)
                                                       gCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCGGGGCAAGCGCTGGGGTGGTGACGCGCCGGCCGACGATCAagcc                                                                                               >  2:2957561/1‑150 (MQ=255)
                                                                 cgcgcggCGCCGGGGGGCCGGTGAGCGGTCCCGCGGGGGGGGGCGCGGGCGGCGCGGGCCGGGGGGGCCGGGGGGGCGCGGCGCGTGGGGGGGGGCGGGGCAAGCGAGGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAAt                                                                                     <  2:5277392/145‑1 (MQ=255)
                                                                  tagtGGCCCCGCGGTCCCGGTGCCCGGTACCGCGGTGGGGGGCGCGGTCGGCGCGGTCCTTGGTGGCATGGGGGGCGAGCCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATc                                                                                    <  2:407861/147‑1 (MQ=255)
                                                                  gaggggCGCCGCGGCGCCGGTGATCGGGACCGCGGTGGGGGGCGCGGCCGGCGCGGCCCCGGGGGGCAGGGGGGGCGAGTCGTCGGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATc                                                                                    <  2:3934848/146‑1 (MQ=255)
                                                                    ttGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATcgc                                                                                  >  1:5499540/1‑150 (MQ=255)
                                                                    ttGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAAAcgc                                                                                  >  2:261007/1‑150 (MQ=255)
                                                                       gCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATcgccgc                                                                               >  1:817114/1‑150 (MQ=255)
                                                                         tCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAATCCGAGCCGGAATCGccgccg                                                                             >  1:2326923/1‑150 (MQ=255)
                                                                         tCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGccgccg                                                                             >  1:3190772/1‑150 (MQ=255)
                                                                               gtgCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATcgccgc                                                                               <  2:4150732/142‑1 (MQ=255)
                                                                               gtgCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATcgccgc                                                                               >  1:4150732/1‑142 (MQ=255)
                                                                               gtgCCGGTGAGCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGGGTGGGGGGTGGCTGGGCAAGCGAGGGTTTGGTGACGAGCCGGCCGCCGATCAAGCCGAGCCGGAATCGCCGCCGGCACTg                                                                       >  2:2829849/1‑150 (MQ=255)
                                                                                ggCCGGCGAGCGGCCCCGCCGGGGGTGGCGCGGCCGGCGCGGTCCCGGGGGGCCGGGGGGGCGAGCCGATCGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTgg                                                                      <  2:965148/149‑1 (MQ=255)
                                                                                   cGGTGAGCGGTACCGCGGTGGGTGGCGCGGTCGGCGCTGTCCGTGGTGGCCTGGGGGGCGAGGCGTGTGGGGGGTGGCTGGGCACGCGATGGTTTGGTGACGAGCCGGCCGCCGATCACGCCGAGCCGGCATCGCCGCCGGCACggggga                                                                   >  2:4978920/1‑148 (MQ=255)
                                                                                    ggTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTaa                                                                  >  1:2856922/1‑150 (MQ=255)
                                                                                     gTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGAGGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTcaa                                                                 >  2:3386936/1‑150 (MQ=255)
                                                                                      ggATCGCCACCGCTCCGGGTGGCGCGGTCGGCGCCGACCCGGGTGGCCTGGCGGGCGAGCCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAg                                                                <  2:4206747/149‑1 (MQ=255)
                                                                                        aTCGGTCCCGCCGTGGGCGGCGCGGCCGGCGCGGCCCCTGGGGGCATAGGGGGCGAGACGACGGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCg                                                              <  1:1451782/150‑1 (MQ=255)
                                                                                           ggTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCACGCCAGGGTTTGGTGACGAGCCGGCCCACGATCAATCCGATCCGGAATCGCCCCCGCCACGGGGGAAAGCCTTa                                                           >  1:487794/1‑150 (MQ=255)
                                                                                            gccccGCTGGGGGGGGCGCGGGCGGCGCCGCCCGTGGTGGCAGGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAc                                                          <  2:4689994/147‑1 (MQ=255)
                                                                                            gTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAc                                                          >  1:787794/1‑150 (MQ=255)
                                                                                            gTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGGTTGGGGGGTGGCTGGGCAACCGATGGTTTGGTGTCGAGCCGGCCGACGATCAACCCGAGCCGGAAACGCCGCCGGCACTGGGTAAAGCGTTAc                                                          >  1:1051736/1‑150 (MQ=255)
                                                                                             tACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGGAAAGCGTTAcg                                                         >  2:5387557/1‑150 (MQ=255)
                                                                                              cccGCTGCGGGTGGCGCGGCCGCCGCTCACCCCGGCCGCCCGCGGGGCGAGTCGAAGGGGGGGCGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAcgc                                                        <  1:921034/149‑1 (MQ=255)
                                                                                              aCCGCTGTGGGTGGCGCGGTCGGCGCTGGCCTTGGTGGCATGGGGGGCGAGCCGCTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAcgc                                                        <  2:449403/150‑1 (MQ=255)
                                                                                                 gCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTa                                                     >  1:4616097/1‑150 (MQ=255)
                                                                                                 gCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTa                                                     >  1:4616130/1‑150 (MQ=255)
                                                                                                  cTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTaa                                                    >  2:2902510/1‑150 (MQ=255)
                                                                                                  cTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTaa                                                    <  2:4485367/150‑1 (MQ=255)
                                                                                                     tgggtggCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGGAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACgg                                                 >  2:1211466/1‑150 (MQ=255)
                                                                                                      gggcggCGCGGTCGGCGCCGCCCTCGGTGGCATGGGGGCCGAGTCGGTTGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGt                                                <  2:1334032/150‑1 (MQ=255)
                                                                                                           gcgcGGCCGGCGCTCCCCTCGCCGGCCGGGGGCGCGACGCGCGTGGGGGGCGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATgc                                           <  2:559489/150‑1 (MQ=255)
                                                                                                             gcgGCCGCCGCGGTCCGTGGGGGCAGGGGGGGCGAGCCGGTTGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATgcca                                         <  2:4558734/150‑1 (MQ=255)
                                                                                                                     cgcTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGa                                 >  2:5110211/1‑150 (MQ=255)
                                                                                                                         gtccggggtggAATGGGGGGCGTCTCGTTGGGGGGGTTGCTGGGCAAGCGATGGTTTTGTGACTAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCATCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTc                             <  2:284381/144‑1 (MQ=255)
                                                                                                                          tCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGGTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCg                            >  2:4307221/1‑150 (MQ=255)
                                                                                                                           ccTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGa                           <  2:3537784/150‑1 (MQ=255)
                                                                                                                            cTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCCAGGAAGTCGAg                          >  1:254800/1‑150 (MQ=255)
                                                                                                                            cTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCCGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAg                          >  1:3893724/1‑150 (MQ=255)
                                                                                                                            cTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCGTGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCCAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAg                          >  2:1890843/1‑150 (MQ=255)
                                                                                                                             ttGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCCGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGGTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGc                         >  2:2512066/1‑150 (MQ=255)
                                                                                                                             ttGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCCACCCACGATAAATCCGAGCCCGAATAGCCGCCGGCACTGGGTAAAGCGTTACGCGTAAACGTGATGCCAGCCAAGGAAGTCGAGc                         >  1:1964061/1‑150 (MQ=255)
                                                                                                                                  ggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAg                    >  2:1327541/1‑150 (MQ=255)
                                                                                                                                   gCATGGGGGGCGAGTCGTTTTGGGGGTGGCTGGGCAAGCGATTGTTTGGTGACGAGCCCGCCGACGATAAATCCGAGCCGGAATCGCCGCCGGCACTGGGTAAATCGGTACGCGTAAAGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGc                   >  2:2344323/1‑150 (MQ=255)
                                                                                                                                   gCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACCAGCCCGCCGACGAACAACCCCAGCCCGAATCGCCGCCCGCACAGGGTAAAGCGTTACGCGTAACCGTTATTCCAGCCAATGAACTCGATCCTGAGc                   >  1:5227288/1‑150 (MQ=255)
                                                                                                                                      tGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAAAAGCCGCCCGCACTGGGTAAAGCGTTACGCGTAACCGTGATGCCACCCAAGGAAGTCGAGCCTGAGcccc                >  1:287889/1‑149 (MQ=255)
                                                                                                                                         ggggCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCa             >  2:2526389/1‑150 (MQ=255)
                                                                                                                                              gAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCt        >  2:2725718/1‑150 (MQ=255)
                                                                                                                                              gAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCt        >  2:4147682/1‑150 (MQ=255)
                                                                                                                                                gTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTcc      >  2:358718/1‑150 (MQ=255)
                                                                                                                                                    tttGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTCAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCCCGGCTCCGGTc  >  2:4289746/1‑150 (MQ=255)
                                                                                                                                                    tttGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTc  >  2:3159036/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTC  >  NC_002947/3448372‑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: