Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,624 G→A 15.2% I438I (ATC→ATT rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,6240GA15.2% 431.0 / 61.3 177I438I (ATC→ATTrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base G (74/76);  new base A (12/15);  total (86/91)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 186 total aligned reads displayed.

TCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAG  >  NC_002947/471475‑471773
                                                                                                                                                     |                                                                                                                                                     
tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGACTGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACgg                                                                                                                                                       >  2:4129258/1‑150 (MQ=255)
 cGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGTTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGa                                                                                                                                                      <  2:946108/150‑1 (MQ=255)
    ccAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATcc                                                                                                                                                   >  1:5049857/1‑150 (MQ=255)
    ccAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATcc                                                                                                                                                   >  1:1080631/1‑150 (MQ=255)
     cAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATcca                                                                                                                                                  <  2:4937954/150‑1 (MQ=255)
        ctctTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATccacca                                                                                                                                               <  1:3222032/150‑1 (MQ=255)
        ctctTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATccacca                                                                                                                                               <  1:272997/150‑1 (MQ=255)
           ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                                                                                                                            >  1:3786181/1‑150 (MQ=255)
            tCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTg                                                                                                                                           <  1:4662507/150‑1 (MQ=255)
               ggAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCa                                                                                                                                        <  2:2726626/150‑1 (MQ=255)
                 aGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATa                                                                                                                                      <  1:1426905/150‑1 (MQ=255)
                  gTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAg                                                                                                                                     >  1:4599518/1‑150 (MQ=255)
                     ggTTCGCGACCCATTTCCTGCAGCATCTGGCGTGAAATACGGTTGATCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCCTGCCTTGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTc                                                                                                                                  <  2:4009411/150‑1 (MQ=255)
                      gTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCg                                                                                                                                 >  2:2325116/1‑150 (MQ=255)
                      gTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCg                                                                                                                                 >  1:3468932/1‑150 (MQ=255)
                        tCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCgag                                                                                                                               <  2:2930646/150‑1 (MQ=255)
                          gcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGaa                                                                                                                             <  2:2180698/150‑1 (MQ=255)
                             aCCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACtg                                                                                                                          >  1:352808/1‑149 (MQ=255)
                             aCCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACtt                                                                                                                          <  2:3805428/150‑1 (MQ=255)
                             aCCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACtt                                                                                                                          <  1:3252758/150‑1 (MQ=255)
                                gaTTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTa                                                                                                                       <  1:815507/149‑1 (MQ=255)
                                 aTTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAg                                                                                                                      >  2:3598059/1‑150 (MQ=255)
                                     ccTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCcgc                                                                                                                  <  1:2183346/150‑1 (MQ=255)
                                       tGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGa                                                                                                                >  2:5200110/1‑150 (MQ=255)
                                         cagcaTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACg                                                                                                              <  2:767274/150‑1 (MQ=255)
                                         cagcaTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACg                                                                                                              >  1:324741/1‑150 (MQ=255)
                                         cagcaTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACg                                                                                                              <  2:214171/150‑1 (MQ=255)
                                             aTCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGatat                                                                                                          >  2:1220327/1‑150 (MQ=255)
                                              tCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAtt                                                                                                         >  2:1998446/1‑150 (MQ=255)
                                                 ggCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGCGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGa                                                                                                      <  1:768235/150‑1 (MQ=255)
                                                   cGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc                                                                                                    <  1:602033/150‑1 (MQ=255)
                                                   cGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc                                                                                                    <  2:2718994/150‑1 (MQ=255)
                                                     ggAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACtt                                                                                                  <  1:4611334/150‑1 (MQ=255)
                                                          tACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCa                                                                                             <  2:4761640/150‑1 (MQ=255)
                                                           aCGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCAc                                                                                            >  2:5279963/1‑150 (MQ=255)
                                                              gTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGc                                                                                         >  1:2620376/1‑150 (MQ=255)
                                                                 gatcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCctt                                                                                      >  2:2697441/4‑150 (MQ=255)
                                                                  atcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGt                                                                                                                        >  1:715739/3‑115 (MQ=255)
                                                                  atcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGt                                                                                                                        <  2:715739/113‑1 (MQ=255)
                                                                    cTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcat                                                                                   <  1:4206843/150‑1 (MQ=255)
                                                                     ttgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatc                                                                                  >  1:4324989/1‑150 (MQ=255)
                                                                      tgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatca                                                                                 <  1:4415542/150‑1 (MQ=255)
                                                                       gttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAg                                                                                >  2:1742702/1‑150 (MQ=255)
                                                                          gATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAcc                                                                             >  1:555078/1‑150 (MQ=255)
                                                                              gTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATg                                                                         >  1:1474087/1‑150 (MQ=255)
                                                                              gTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATg                                                                         >  1:1546532/1‑150 (MQ=255)
                                                                              gTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATg                                                                         >  1:2944854/1‑150 (MQ=255)
                                                                               tctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTTTTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                                        <  2:4604970/150‑1 (MQ=255)
                                                                                ttcGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATAtt                                                                       <  2:4277742/149‑1 (MQ=255)
                                                                                          tgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                                             <  2:626336/150‑1 (MQ=255)
                                                                                          tgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                                             >  1:389662/1‑150 (MQ=255)
                                                                                            cACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGAt                                                           >  2:4430935/1‑150 (MQ=255)
                                                                                             aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATc                                                          <  1:3642127/150‑1 (MQ=255)
                                                                                             aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATc                                                          >  1:2378438/1‑150 (MQ=255)
                                                                                              ccGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCa                                                         <  1:3856609/150‑1 (MQ=255)
                                                                                              ccGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCa                                                         <  2:1580729/150‑1 (MQ=255)
                                                                                              ccGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCa                                                         >  2:2519252/1‑150 (MQ=255)
                                                                                               cGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAg                                                        >  1:3184156/1‑150 (MQ=255)
                                                                                                  aaTGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTc                                                     <  2:4324989/150‑1 (MQ=255)
                                                                                                    tacGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGa                                                   <  2:2120921/148‑1 (MQ=255)
                                                                                                    tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                   <  2:1080631/150‑1 (MQ=255)
                                                                                                    tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                   >  1:4551397/1‑150 (MQ=255)
                                                                                                    tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                   <  1:4029762/150‑1 (MQ=255)
                                                                                                    tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                   <  1:5279963/150‑1 (MQ=255)
                                                                                                    tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                   <  2:2475398/150‑1 (MQ=255)
                                                                                                     gCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAg                                                  <  2:1306317/150‑1 (MQ=255)
                                                                                                     gCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAg                                                  <  2:2635398/150‑1 (MQ=255)
                                                                                                     gCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAg                                                  <  2:921313/150‑1 (MQ=255)
                                                                                                     gCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAg                                                  <  2:3105392/150‑1 (MQ=255)
                                                                                                       ggATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGa                                                >  2:1775890/1‑150 (MQ=255)
                                                                                                        gATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGaa                                               <  1:731387/150‑1 (MQ=255)
                                                                                                         aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                              <  1:4400983/150‑1 (MQ=255)
                                                                                                         aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                              <  2:3507806/150‑1 (MQ=255)
                                                                                                         aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                              >  1:4708495/1‑150 (MQ=255)
                                                                                                              gcgtgcCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAg                                         <  1:2157765/150‑1 (MQ=255)
                                                                                                                       gTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtc                                                                                     <  2:1239666/97‑1 (MQ=255)
                                                                                                                       gTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGtt                                >  1:1098123/1‑150 (MQ=255)
                                                                                                                       atCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGtt                                <  1:3353532/149‑1 (MQ=255)
                                                                                                                          ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTTGCCGCGTCGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGAGGTTGCCTTCCTGGATCAGGTCGTGGACCTGCCGGCCTCGGTTGgt                             >  2:1414089/1‑150 (MQ=255)
                                                                                                                          ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgt                             >  1:4448453/1‑150 (MQ=255)
                                                                                                                           gCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAGTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgtg                            >  2:3497642/1‑150 (MQ=255)
                                                                                                                            cGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgtgt                           >  2:883071/1‑150 (MQ=255)
                                                                                                                             gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                          <  1:3467782/150‑1 (MQ=255)
                                                                                                                              aTCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTAc                         >  1:5137688/1‑150 (MQ=255)
                                                                                                                                  aaCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtct                     >  2:2439938/1‑150 (MQ=255)
                                                                                                                                    cGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTg                   <  1:2537292/150‑1 (MQ=255)
                                                                                                                                      ggTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGc                 >  2:1926148/1‑150 (MQ=255)
                                                                                                                                          aTCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATg             <  2:2335006/150‑1 (MQ=255)
                                                                                                                                              ccTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGaga         <  1:2226926/150‑3 (MQ=255)
                                                                                                                                              ccTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGaga         <  1:1433117/150‑3 (MQ=255)
                                                                                                                                                 gACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa      >  1:2820951/1‑150 (MQ=255)
                                                                                                                                                 gACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa      >  1:4739836/1‑150 (MQ=255)
                                                                                                                                                 gACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAtaa      >  1:289984/1‑148 (MQ=255)
                                                                                                                                                 gACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAtaa      <  1:2661586/150‑3 (MQ=255)
                                                                                                                                                   cGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCAcc    <  2:4297909/150‑1 (MQ=255)
                                                                                                                                                    ggATCCACCTGGTTGAATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa   <  1:1884146/150‑1 (MQ=255)
                                                                                                                                                    ggATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa   <  1:407821/150‑1 (MQ=255)
                                                                                                                                                    ggATCCACCAGGTGGCATAGGTCGAGAACTTGGAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa   <  2:3393797/150‑1 (MQ=255)
                                                                                                                                                     gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg  >  2:529977/1‑150 (MQ=255)
                                                                                                                                                     gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg  >  1:4513983/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAG  >  NC_002947/471475‑471773

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: