Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,438,169 T→A 7.2% W47R (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,1690TA 7.2% 786.4 / 35.1 236W47R (TGG→AGG) gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base T (104/115);  new base A (7/10);  total (111/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.77e-01

Only 100 of 238 total aligned reads displayed.

AAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATT  >  NC_002947/3438020‑3438316
                                                                                                                                                      |                                                                                                                                                   
 aGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTg                                                                                                                                                    <  2:3198073/150‑1 (MQ=255)
  ggtggtGACATGAGCCTGTTACCGCACCACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACCTTAGCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                                   <  2:3054350/150‑1 (MQ=255)
  ggtggtGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                                   <  1:987029/150‑1 (MQ=255)
  ggtggtGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                                   <  2:2319250/150‑1 (MQ=255)
  ggtggtGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                                   <  1:4489536/150‑1 (MQ=255)
    tggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGc                                                                                                                                                 >  2:4410080/1‑150 (MQ=255)
     ggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCa                                                                                                                                                >  2:3810868/1‑150 (MQ=255)
      gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCAt                                                                                                                                               <  1:2514476/150‑1 (MQ=255)
      gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCAt                                                                                                                                               <  1:674497/150‑1 (MQ=255)
      gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCAt                                                                                                                                               <  1:2105277/150‑1 (MQ=255)
        gACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                             >  2:3516083/1‑150 (MQ=255)
          cATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGaa                                                                                                                                           >  2:4117998/1‑150 (MQ=255)
              aGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTc                                                                                                                                       <  2:1564981/150‑1 (MQ=38)
              aGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTc                                                                                                                                       <  1:2155677/150‑1 (MQ=255)
                 cTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAg                                                                                                                                    <  2:2601899/150‑1 (MQ=255)
                  tGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAgg                                                                                                                                   >  1:3699340/1‑150 (MQ=255)
                   ggtACCGCCCACCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTGGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAggg                                                                                                                                  <  2:1341002/148‑1 (MQ=255)
                    ttACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGc                                                                                                                                 >  2:3935394/1‑150 (MQ=255)
                       ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                              >  2:3231558/1‑150 (MQ=255)
                         gCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGa                                                                                                                            <  1:3486319/150‑1 (MQ=255)
                          cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                             >  2:145381/1‑132 (MQ=255)
                          cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                             <  1:145381/132‑1 (MQ=255)
                           acaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGAGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAgg                                                                                                                          >  2:2747885/1‑150 (MQ=255)
                            caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                                         >  2:1466573/1‑150 (MQ=255)
                            caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                                         <  1:1179400/150‑1 (MQ=255)
                             aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcg                                                                                                                        >  1:1623248/1‑150 (MQ=255)
                                  cACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTa                                                                                                                   >  1:3367130/1‑150 (MQ=255)
                                   aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTac                                                                                                                  <  1:4325725/150‑1 (MQ=255)
                                   aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTac                                                                                                                  <  2:623090/150‑1 (MQ=255)
                                   aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTACCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGTTGGGAGGCGGCTTac                                                                                                                  <  1:1784931/150‑1 (MQ=255)
                                    cccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                                 <  2:2862674/150‑1 (MQ=255)
                                              gcgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGa                                                                                                                            <  1:3735913/129‑1 (MQ=255)
                                               cgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAg                                                                                                      <  2:1770892/150‑1 (MQ=255)
                                                 cgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAgcg                                                                                                    <  1:3183001/150‑1 (MQ=255)
                                                 cgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAgcg                                                                                                    <  1:2600661/150‑1 (MQ=255)
                                                     cTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgccgc                                                                                                <  1:3454295/150‑1 (MQ=255)
                                                     cTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGACGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGGCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgccgc                                                                                                >  2:67464/1‑150 (MQ=255)
                                                      tGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcgcg                                                                                               >  2:4179295/1‑150 (MQ=255)
                                                        gAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcgcgcg                                                                                             >  1:4346204/1‑150 (MQ=255)
                                                             ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGCCTC‑‑‑‑CCAGCCCTCGACCTTCCTGCCCTGGCTGGCCTgg                                                                                                                                                   <  2:1025136/85‑1 (MQ=16)
                                                             ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGAt                                                                                        >  1:565321/1‑150 (MQ=255)
                                                                 gccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                                    <  2:1632228/150‑1 (MQ=255)
                                                                 gccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                                    >  2:1560793/1‑150 (MQ=255)
                                                                   cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                                  <  1:1825977/150‑1 (MQ=255)
                                                                   cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                                  >  1:360235/1‑150 (MQ=255)
                                                                    gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                                 <  1:2014517/150‑1 (MQ=255)
                                                                        ttGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCa                                                                             >  1:150134/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:2005110/1‑150 (MQ=255)
                                                                                gACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                     <  2:1289537/150‑1 (MQ=255)
                                                                                  cccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                                   >  1:4489420/1‑150 (MQ=255)
                                                                                   ccGGACATTATCCGGGGCGTGGCCGACTGGTCACGCTGCCCGCCGAACTACCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                                  >  1:337596/1‑150 (MQ=255)
                                                                                   ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                                  <  2:1586509/150‑1 (MQ=255)
                                                                                     ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                                >  2:1335980/1‑150 (MQ=255)
                                                                                     ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                                <  1:1691542/150‑1 (MQ=255)
                                                                                     ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                                <  1:720678/150‑1 (MQ=255)
                                                                                          ttATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCa                                                           >  1:104854/1‑150 (MQ=255)
                                                                                            aTCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAg                                                         <  1:2797073/150‑1 (MQ=255)
                                                                                            aTCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAg                                                         <  1:2717240/150‑1 (MQ=255)
                                                                                             tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                                        <  1:4264922/150‑1 (MQ=255)
                                                                                             tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                                        <  2:3873337/150‑1 (MQ=255)
                                                                                             tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                                        >  2:2441112/1‑150 (MQ=255)
                                                                                              ccGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCCTTCCGGTTCACAAGACCAAggg                                                       >  1:4170052/1‑150 (MQ=255)
                                                                                                     gTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTc                                                <  2:1498262/150‑1 (MQ=255)
                                                                                                         ccGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                            <  1:2180839/150‑1 (MQ=255)
                                                                                                          cGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcgc                                           <  1:3261438/150‑1 (MQ=255)
                                                                                                             cTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGAt                                        <  1:100802/150‑1 (MQ=255)
                                                                                                                 tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcgg                                    >  2:2135066/1‑150 (MQ=255)
                                                                                                                  cACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggc                                   >  2:1846693/1‑150 (MQ=255)
                                                                                                                    cGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                 <  2:3636830/150‑1 (MQ=255)
                                                                                                                    cGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                 <  1:3810868/150‑1 (MQ=255)
                                                                                                                     gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCCTTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCggg                                >  2:2908303/1‑150 (MQ=255)
                                                                                                                      cTGCCCGCCGAACGTCCTGCCCTGGCTGGCCTGGGCAAGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGt                               <  2:3380588/150‑1 (MQ=255)
                                                                                                                       cagccctcgACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCATGGTCGGCGCGATTCGGCGGGTg                              >  2:862106/8‑150 (MQ=25)
                                                                                                                        gCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGc                             <  2:521841/150‑1 (MQ=255)
                                                                                                                          ccgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                           <  1:2878846/150‑1 (MQ=255)
                                                                                                                          ccgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGACATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                           <  2:3328129/150‑1 (MQ=255)
                                                                                                                           cgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAGTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCa                          >  2:4086678/1‑150 (MQ=255)
                                                                                                                            gccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCaa                         <  1:2094131/150‑1 (MQ=255)
                                                                                                                            gccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCaa                         <  1:2223064/150‑1 (MQ=255)
                                                                                                                            ctcgACCTTCCTGCCCTGGCTGGCCTGGGCCTGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCCTTCCGGTTCACAAGACCAAGGGCCCGGTCGGCGCGAGTCGGCGGGTGCTCaa                         >  2:1692558/3‑150 (MQ=255)
                                                                                                                                 aCTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCgg                    <  1:2572002/150‑1 (MQ=255)
                                                                                                                                  cgtCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGt                   <  1:4074483/148‑1 (MQ=255)
                                                                                                                                  cTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGt                   <  1:956770/150‑1 (MQ=255)
                                                                                                                                   ttccatGCCCAGGCCAGCCAGGGCAGGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGt                   >  1:3071074/6‑150 (MQ=255)
                                                                                                                                   ttCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGtt                  >  1:1321095/1‑150 (MQ=255)
                                                                                                                                    tCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTc                 <  2:3699340/150‑1 (MQ=255)
                                                                                                                                     ccTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg                >  2:2827565/1‑150 (MQ=255)
                                                                                                                                      cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg               <  2:2203282/150‑1 (MQ=255)
                                                                                                                                      cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg               <  2:2374302/150‑1 (MQ=255)
                                                                                                                                      cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg               <  1:1692965/150‑1 (MQ=255)
                                                                                                                                       atGCCCAGGCCAGCCAGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg               >  2:797518/2‑150 (MQ=255)
                                                                                                                                        tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCggg              >  1:1828703/1‑150 (MQ=255)
                                                                                                                                         gCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGt             >  1:4219642/1‑150 (MQ=255)
                                                                                                                                         gCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGt             >  2:2525625/1‑150 (MQ=255)
                                                                                                                                          cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTc            >  1:240291/1‑150 (MQ=255)
                                                                                                                                                  ggcCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGa    >  2:3456776/1‑150 (MQ=255)
                                                                                                                                                    ccTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGAtt  >  1:4016116/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATT  >  NC_002947/3438020‑3438316

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: