Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,438,160 T→A 7.2% W44R (TGG→AGG)  gpI → tail formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,438,1600TA 7.2% 754.0 / 35.3 235W44R (TGG→AGG) gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base T (103/114);  new base A (7/10);  total (110/125)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 249 total aligned reads displayed.

CTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTC  >  NC_002947/3438011‑3438306
                                                                                                                                                      |                                                                                                                                                  
 tGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCctg                                                                                                                                                   >  1:1998156/1‑150 (MQ=255)
   aCGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGATATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCctggc                                                                                                                                                 >  1:2203155/1‑150 (MQ=255)
        gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGcc                                                                                                                                            <  1:2661696/150‑1 (MQ=255)
         aaGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCt                                                                                                                                           >  2:3657245/1‑150 (MQ=255)
          aGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTg                                                                                                                                          <  2:3768745/150‑1 (MQ=255)
           ggtggtGACATGAGCCTGTTACCGCACCACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACCTTAGCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                         <  2:3054350/150‑1 (MQ=255)
           ggtggtGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                         <  1:4489536/150‑1 (MQ=255)
             tggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGc                                                                                                                                       >  2:4410080/1‑150 (MQ=255)
               gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCAt                                                                                                                                     <  1:2514476/150‑1 (MQ=255)
               gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCAt                                                                                                                                     <  1:2105269/150‑1 (MQ=255)
                  aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCAGGCCAGCCAGGGCaggaa                                                                                                                                  <  2:4497642/150‑5 (MQ=255)
                   cATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGaa                                                                                                                                 >  2:4117998/1‑150 (MQ=255)
                   cATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGaa                                                                                                                                 <  1:863872/150‑1 (MQ=255)
                      gAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGt                                                                                                                              >  2:3158669/1‑150 (MQ=255)
                       aGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTc                                                                                                                             <  2:1565606/150‑1 (MQ=38)
                         ccTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGa                                                                                                                           <  1:998337/150‑1 (MQ=255)
                          cTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCTCGACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCggcg                                                                                                                         <  2:2375588/150‑4 (MQ=25)
                           tGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAgg                                                                                                                         >  1:3699340/1‑150 (MQ=255)
                             ttACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGc                                                                                                                       >  2:3935394/1‑150 (MQ=255)
                                ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                    >  2:3231558/1‑150 (MQ=255)
                                ccGCACAACGCCACCCAGGTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCGGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                    <  1:2029448/150‑1 (MQ=255)
                                  gCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGa                                                                                                                  <  1:3486319/150‑1 (MQ=255)
                                   cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                   >  2:145381/1‑132 (MQ=255)
                                     caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                               >  2:1466587/1‑150 (MQ=255)
                                     caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                               <  1:1179400/150‑1 (MQ=255)
                                      aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcg                                                                                                              >  1:1623259/1‑150 (MQ=255)
                                        cGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                            <  2:4004738/150‑1 (MQ=255)
                                        cGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                            <  1:3875499/150‑1 (MQ=255)
                                           cACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTa                                                                                                         >  1:661/1‑150 (MQ=255)
                                            aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTac                                                                                                        >  2:2769962/1‑150 (MQ=255)
                                            aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTac                                                                                                        <  1:4325725/150‑1 (MQ=255)
                                            aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTac                                                                                                        <  2:623090/150‑1 (MQ=255)
                                                       gcgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGa                                                                                                                  <  1:3735913/129‑1 (MQ=255)
                                                        cgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAg                                                                                            <  2:1770538/150‑1 (MQ=255)
                                                          cgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAgcg                                                                                          <  1:2600661/150‑1 (MQ=255)
                                                            cccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgcc                                                                                        <  2:2600577/150‑1 (MQ=255)
                                                              cTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGACGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGGCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgccgc                                                                                      >  2:67464/1‑150 (MQ=255)
                                                                 gAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcgcgcg                                                                                   >  1:4346204/1‑150 (MQ=255)
                                                                  cggcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                                  <  1:312888/149‑1 (MQ=255)
                                                                    gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTTGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTg                                                                                >  2:4191746/1‑150 (MQ=255)
                                                                    gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTg                                                                                >  1:1700309/1‑150 (MQ=255)
                                                                      ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGCCTC‑‑‑‑CCAGCCCTCGACCTTCCTGCCCTGGCTGGCCTgg                                                                                                                                         <  2:1025136/85‑1 (MQ=16)
                                                                      ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGAt                                                                              >  1:565692/1‑150 (MQ=255)
                                                                       gccgccGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATc                                                                             >  2:1603176/1‑150 (MQ=255)
                                                                          gccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                          >  2:1560795/1‑150 (MQ=255)
                                                                            cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                        >  1:360235/1‑150 (MQ=255)
                                                                                 ttGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCa                                                                   >  1:150134/1‑150 (MQ=255)
                                                                                 ttGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCa                                                                   >  2:3462710/1‑150 (MQ=255)
                                                                                     ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTcc                                                               >  1:132929/1‑150 (MQ=255)
                                                                                      cTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                              >  1:1902274/1‑150 (MQ=255)
                                                                                      cTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                              >  2:4253627/1‑150 (MQ=255)
                                                                                        tGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGt                                                            >  2:2005316/1‑150 (MQ=255)
                                                                                        tGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGt                                                            >  2:2005093/1‑150 (MQ=255)
                                                                                         gACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                           <  2:1289537/150‑1 (MQ=255)
                                                                                          aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGc                                                                                                                                       <  1:2867755/73‑1 (MQ=255)
                                                                                          aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGc                                                                                                                                       >  2:2867755/1‑73 (MQ=255)
                                                                                           cccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                         >  1:4489420/1‑150 (MQ=255)
                                                                                           cccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                         <  1:1469900/150‑1 (MQ=255)
                                                                                            ccGGACATTATCCGGGGCGTGGCCGACTGGTCACGCTGCCCGCCGAACTACCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                        >  1:337596/1‑150 (MQ=255)
                                                                                            ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGACGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                        >  2:2725671/1‑150 (MQ=255)
                                                                                            ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                        >  1:2263532/1‑150 (MQ=255)
                                                                                              ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                      <  1:720678/150‑1 (MQ=255)
                                                                                              ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                      <  1:2548326/150‑1 (MQ=255)
                                                                                              ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                      >  2:1335980/1‑150 (MQ=255)
                                                                                              ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                      <  1:1691542/150‑1 (MQ=255)
                                                                                                  aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAcc                                                  >  1:2970170/1‑150 (MQ=255)
                                                                                                   ttATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCa                                                 >  1:104854/1‑150 (MQ=255)
                                                                                                      tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                              <  1:4264834/150‑1 (MQ=255)
                                                                                                       ccGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCCTTCCGGTTCACAAGACCAAggg                                             >  1:4170052/1‑150 (MQ=255)
                                                                                                       ccGGGGCGCGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAggg                                             <  2:1397355/150‑1 (MQ=255)
                                                                                                              gTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTc                                      <  2:1498262/150‑1 (MQ=255)
                                                                                                                  ccGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                  <  1:2180839/150‑1 (MQ=255)
                                                                                                                  ccGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                  <  1:334478/150‑1 (MQ=255)
                                                                                                                   cGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcgc                                 <  1:1793971/150‑1 (MQ=255)
                                                                                                                      cTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGAt                              <  1:100802/150‑1 (MQ=255)
                                                                                                                       tCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGAtt                             <  1:4268230/150‑1 (MQ=255)
                                                                                                                          tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGTGGGCGGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcgg                          >  2:2323458/1‑150 (MQ=255)
                                                                                                                          tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcgg                          >  1:2211731/1‑150 (MQ=255)
                                                                                                                            aCGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcg                        >  1:3243326/1‑150 (MQ=255)
                                                                                                                             cGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                       <  1:3810868/150‑1 (MQ=255)
                                                                                                                               cTGCCCGCCGAACGTCCTGCCCTGGCTGGCCTGGGCAAGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGt                     <  2:3380588/150‑1 (MQ=255)
                                                                                                                                cagccctcgACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCATGGTCGGCGCGATTCGGCGGGTg                    >  2:862106/8‑150 (MQ=25)
                                                                                                                                 gCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGc                   >  1:4119418/1‑150 (MQ=255)
                                                                                                                                   ccgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGACATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                 <  2:3328129/150‑1 (MQ=255)
                                                                                                                                    cgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAGTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCa                >  2:4086678/1‑150 (MQ=255)
                                                                                                                                           cTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGt         <  1:956770/150‑1 (MQ=255)
                                                                                                                                            ttCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGtt        >  1:1321095/1‑150 (MQ=255)
                                                                                                                                             tccatGCCCAGGCCAGCCAGGGCAGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTc       >  2:1808037/5‑150 (MQ=255)
                                                                                                                                             tCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTc       <  2:2522551/150‑1 (MQ=255)
                                                                                                                                              ccTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg      >  2:2827576/1‑150 (MQ=255)
                                                                                                                                               cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     <  2:3842313/150‑1 (MQ=255)
                                                                                                                                               cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     <  1:2133860/150‑1 (MQ=255)
                                                                                                                                               cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     <  2:2203350/150‑1 (MQ=255)
                                                                                                                                               cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     <  2:2203155/150‑1 (MQ=255)
                                                                                                                                               cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     <  2:3565771/150‑1 (MQ=255)
                                                                                                                                                atGCCCAGGCCAGCCAGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     >  2:797518/2‑150 (MQ=255)
                                                                                                                                                 tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCggg    <  1:2747885/150‑1 (MQ=255)
                                                                                                                                                 tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCggg    >  2:1651575/1‑150 (MQ=255)
                                                                                                                                                 tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCggg    >  1:1828703/1‑150 (MQ=255)
                                                                                                                                                 tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCggg    <  1:3426796/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTC  >  NC_002947/3438011‑3438306

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