Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,716 G→A 5.7% D241D (GAC→GAT PP_3503 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,7160GA 5.7% 935.9 / 29.4 296D241D (GAC→GATPP_3503sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (146/133);  new base A (10/7);  total (156/140)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 301 total aligned reads displayed.

AACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGC  >  NC_002947/3975568‑3975865
                                                                                                                                                    |                                                                                                                                                     
    cccTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGcc                                                                                                                                                  <  1:8169414/150‑1 (MQ=255)
     ccTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCg                                                                                                                                                 >  1:3394926/1‑150 (MQ=255)
           tCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtca                                                                                                                                           <  1:6824917/150‑1 (MQ=255)
           tCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtca                                                                                                                                           >  1:1936789/1‑150 (MQ=255)
           tCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtca                                                                                                                                           >  2:5264911/1‑150 (MQ=255)
           tCGAAGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGAAAACGTAAATCCACCGGTATGGGCAGGCTGCCACCCACCCGCTCGATGTGCTTATACTGCAAGAAACGCAGCAAGTTGGACTGCAATTCAAACGGCAGGTAGCCGATTtca                                                                                                                                           >  2:7381522/1‑150 (MQ=255)
            cGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcat                                                                                                                                          <  2:5632323/150‑1 (MQ=255)
             gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatc                                                                                                                                         >  2:8912848/1‑150 (MQ=255)
             gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatc                                                                                                                                         <  1:8164247/150‑1 (MQ=255)
             gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatc                                                                                                                                         <  2:2113349/150‑1 (MQ=255)
             gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatc                                                                                                                                         <  2:6437138/150‑1 (MQ=255)
             gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcat                                                                                                                                          >  1:3231457/1‑149 (MQ=255)
              aTGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatca                                                                                                                                        <  1:3978927/150‑1 (MQ=255)
               tGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCaa                                                                                                                                       <  1:1598890/150‑1 (MQ=255)
                 gCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGa                                                                                                                                     <  2:711767/150‑1 (MQ=255)
                      cTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACaa                                                                                                                                >  1:8564761/1‑150 (MQ=255)
                       tCAAGGTCCTCATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAAt                                                                                                                               <  1:7401313/150‑1 (MQ=255)
                        cAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAAtg                                                                                                                              >  2:8189113/1‑150 (MQ=255)
                         aaGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAAtgt                                                                                                                             >  1:8804331/1‑150 (MQ=255)
                            gTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTgcc                                                                                                                          <  2:3296266/150‑1 (MQ=255)
                              ccACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTgccgc                                                                                                                        <  2:4358437/150‑1 (MQ=255)
                               cacaTGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTgccgcc                                                                                                                       >  2:1069282/1‑150 (MQ=255)
                                acaTGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCa                                                                                                                      >  2:6723429/1‑150 (MQ=255)
                                acaTGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCa                                                                                                                      >  1:3797931/1‑150 (MQ=255)
                                  aTGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATc                                                                                                                    <  2:1125899/150‑1 (MQ=255)
                                   tGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCg                                                                                                                   <  2:3054715/150‑1 (MQ=255)
                                    gggTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCgg                                                                                                                  <  1:5118380/150‑1 (MQ=255)
                                     ggTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGc                                                                                                                 <  1:8971521/150‑1 (MQ=255)
                                         gccgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCt                                                                                                             >  1:9015509/1‑150 (MQ=255)
                                          ccgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCtt                                                                                                            <  1:7139278/150‑1 (MQ=255)
                                           cgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTc                                                                                                           >  2:6849247/1‑150 (MQ=255)
                                             ccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGa                                                                                                         <  2:269730/150‑1 (MQ=255)
                                                  aCCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGc                                                                                                    <  2:4507831/150‑1 (MQ=255)
                                                    ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                                  >  2:4670608/1‑150 (MQ=255)
                                                    ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                                  >  2:1245246/1‑150 (MQ=255)
                                                    ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                                  >  2:3004359/1‑150 (MQ=255)
                                                    ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                                  <  1:6282618/150‑1 (MQ=255)
                                                       tACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGt                                                                                               <  2:3814365/150‑1 (MQ=255)
                                                       tACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGt                                                                                               <  1:4908920/150‑1 (MQ=255)
                                                             cACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTgcg                                                                                         >  1:4370520/1‑150 (MQ=255)
                                                                 ggTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTgg                                                                                     <  1:6658871/150‑1 (MQ=255)
                                                                  gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                                    <  1:5257108/150‑1 (MQ=255)
                                                                    aTGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGa                                                                                  <  1:7354805/150‑1 (MQ=255)
                                                                       ggctggctGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCa                                                                               <  1:1079321/150‑1 (MQ=255)
                                                                          tggctgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCg                                                                            >  2:771854/1‑150 (MQ=255)
                                                                              tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaa                                                                        >  2:5379586/1‑150 (MQ=255)
                                                                              tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaa                                                                        >  2:1448328/1‑150 (MQ=255)
                                                                               gCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaaa                                                                       <  1:6426503/150‑1 (MQ=255)
                                                                                ccacccccCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCACTTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCGCCGGTAAAg                                                                      >  2:5063613/1‑150 (MQ=255)
                                                                                  acccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgc                                                                    <  1:902142/150‑1 (MQ=255)
                                                                                    ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                  >  1:4213606/1‑150 (MQ=255)
                                                                                    ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                  >  2:7467150/1‑150 (MQ=255)
                                                                                    ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                  >  2:2323094/1‑150 (MQ=255)
                                                                                      acccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGccc                                                                >  2:2517335/1‑150 (MQ=255)
                                                                                         cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                             <  1:2632476/150‑1 (MQ=255)
                                                                                         cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                             <  2:3574630/150‑1 (MQ=255)
                                                                                         cGCTCGAGGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                             <  2:7309969/150‑1 (MQ=255)
                                                                                           cTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtctc                                                           <  1:3098447/150‑1 (MQ=255)
                                                                                           cTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtctc                                                           <  1:5397973/150‑1 (MQ=255)
                                                                                           cTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtctc                                                           <  1:3624446/150‑1 (MQ=255)
                                                                                            tCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCa                                                          <  2:8356026/150‑1 (MQ=255)
                                                                                             cGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCAt                                                         <  1:6555247/150‑1 (MQ=255)
                                                                                              cctgtgCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCACTTCTCATg                                                        >  2:3643012/3‑150 (MQ=255)
                                                                                                    cTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGaa                                                  <  2:4608688/150‑1 (MQ=255)
                                                                                                    cTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGaa                                                  >  2:8427292/1‑150 (MQ=255)
                                                                                                      tctcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACa                                                <  2:3715167/150‑1 (MQ=255)
                                                                                                        tcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACatt                                              <  1:7686272/150‑3 (MQ=255)
                                                                                                          cTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTc                                            <  2:8000028/150‑1 (MQ=255)
                                                                                                          cTGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTc                                            <  2:4217277/150‑1 (MQ=255)
                                                                                                           tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCACGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGACCAGGTCg                                           >  2:6984955/1‑150 (MQ=255)
                                                                                                           tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGc                                                                                                    <  1:8225438/93‑1 (MQ=255)
                                                                                                           tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                           >  1:8104982/1‑150 (MQ=255)
                                                                                                             cAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGa                                         <  1:7213503/150‑1 (MQ=255)
                                                                                                             cAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGa                                         <  1:2541753/150‑1 (MQ=255)
                                                                                                             cAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGa                                         >  1:8720388/1‑150 (MQ=255)
                                                                                                              aaGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGAAAAGGCGCCCTTCTCATGGCCGAACAGATCGGAc                                        >  1:5852322/1‑150 (MQ=255)
                                                                                                                 aaaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTgg                                     >  2:606387/1‑150 (MQ=255)
                                                                                                                  aaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGa                                    >  1:578901/1‑150 (MQ=255)
                                                                                                                   aCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGAt                                   >  2:1320549/1‑150 (MQ=255)
                                                                                                                      cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                <  1:145573/150‑1 (MQ=255)
                                                                                                                       agcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAg                               >  2:8561027/1‑150 (MQ=255)
                                                                                                                       agcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAg                               <  1:97502/150‑1 (MQ=255)
                                                                                                                        gcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGt                              >  2:5699667/1‑150 (MQ=255)
                                                                                                                            gTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTc                          <  2:6651561/150‑1 (MQ=255)
                                                                                                                            gTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCGTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTc                          <  1:1342319/150‑1 (MQ=255)
                                                                                                                              tGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCgg                        <  1:2767304/150‑1 (MQ=255)
                                                                                                                               ggCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCggg                       >  1:1220054/1‑150 (MQ=255)
                                                                                                                                 ccTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGa                     >  2:7603920/1‑150 (MQ=255)
                                                                                                                                  cTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCgggg                      >  1:4187767/1‑148 (MQ=255)
                                                                                                                                  cTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGAt                    >  1:5050541/1‑150 (MQ=255)
                                                                                                                                   tGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATc                   >  1:2365648/1‑150 (MQ=255)
                                                                                                                                    gCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCg                  >  2:8274883/1‑150 (MQ=255)
                                                                                                                                        ttCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcccc              <  2:2429185/150‑1 (MQ=255)
                                                                                                                                         tCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCg             >  2:1873221/1‑150 (MQ=255)
                                                                                                                                           cAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCa           <  2:1723083/150‑1 (MQ=255)
                                                                                                                                           cAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCa           <  2:3394926/150‑1 (MQ=255)
                                                                                                                                            aaCGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAg          >  2:8710974/1‑150 (MQ=255)
                                                                                                                                               ggCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTg       >  1:8338064/1‑150 (MQ=255)
                                                                                                                                                    gTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATcgc  <  1:2111153/150‑1 (MQ=255)
                                                                                                                                                    gTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATcgc  >  2:2032365/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGC  >  NC_002947/3975568‑3975865

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: