Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,686,791 T→A 5.3% L29* (TTG→TAG) ‡ mltD → membrane‑bound lytic murein transglycosylase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,686,7910TA 5.3% 1072.3 / 34.0 322L29* (TTG→TAG) ‡mltDmembrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base T (157/147);  new base A (10/7);  total (168/154)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 350 total aligned reads displayed.

TGACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAG  >  NC_002947/4686642‑4686939
                                                                                                                                                      |                                                                                                                                                    
  aCCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtgg                                                                                                                                                    <  1:51902/150‑1 (MQ=255)
    cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                                  <  1:1658750/150‑1 (MQ=255)
    cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                                  <  1:3814748/150‑1 (MQ=255)
    cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                                  <  2:6520255/150‑1 (MQ=255)
    cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                                  <  1:3479778/150‑1 (MQ=255)
     aGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtgg                                                                                                                                                 <  1:1005411/150‑1 (MQ=255)
     aGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtgg                                                                                                                                                 <  1:811266/150‑1 (MQ=255)
      ggCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTggg                                                                                                                                                <  2:3734237/150‑1 (MQ=255)
        cccGGAGCGGTGTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCt                                                                                                                                              <  2:4853589/150‑1 (MQ=255)
          cGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGc                                                                                                                                            <  1:4029751/150‑1 (MQ=255)
             aGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCag                                                                                                                                         <  2:1601999/150‑1 (MQ=255)
              gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGTCACCTTGGTGGGCTGCCaga                                                                                                                                        <  2:1940619/150‑1 (MQ=255)
              gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                        <  1:6191629/150‑1 (MQ=255)
              gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                        <  2:5088480/150‑1 (MQ=255)
                  tttCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAc                                                                                                                                    <  2:5496018/150‑1 (MQ=255)
                   ttCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAcc                                                                                                                                   <  1:841545/150‑1 (MQ=255)
                   ttCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAcc                                                                                                                                   >  1:13519/1‑150 (MQ=255)
                    tCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                                  >  2:5065683/1‑150 (MQ=255)
                     cTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCg                                                                                                                                 >  2:1694326/1‑150 (MQ=255)
                      tAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                                >  2:61965/1‑150 (MQ=255)
                      tAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCCGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                                <  2:2993959/150‑1 (MQ=255)
                        gAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAACCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                                <  1:294298/148‑1 (MQ=255)
                         aaTCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAg                                                                                                                             >  1:4548822/1‑150 (MQ=255)
                         aaTCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAg                                                                                                                             >  2:2701583/1‑150 (MQ=255)
                         aaTCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAg                                                                                                                             <  1:1768643/150‑1 (MQ=255)
                         aaTCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAg                                                                                                                             >  1:5867713/1‑150 (MQ=255)
                            cGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTc                                                                                                                          <  2:5654646/150‑1 (MQ=255)
                            cGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTc                                                                                                                          <  2:1163738/150‑1 (MQ=255)
                              ccGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                                        >  2:7084711/1‑150 (MQ=255)
                              ccGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                                        >  1:7172004/1‑150 (MQ=255)
                               cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                       <  2:2853075/150‑1 (MQ=255)
                               cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCAGTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                       <  2:6474548/150‑1 (MQ=255)
                                gAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgacg                                                                                                                      >  2:4231616/1‑150 (MQ=255)
                                 aaCTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgacga                                                                                                                     <  2:1827780/150‑1 (MQ=255)
                                    tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATc                                                                                                                  >  1:260792/1‑150 (MQ=255)
                                      tttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCg                                                                                                                >  2:3075579/1‑150 (MQ=255)
                                      tttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCg                                                                                                                >  2:3163750/1‑150 (MQ=255)
                                         cGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATa                                                                                                             <  1:2506101/150‑1 (MQ=255)
                                         cGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATa                                                                                                             <  2:3965421/150‑1 (MQ=255)
                                          gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAg                                                                                                            >  1:5266851/1‑150 (MQ=255)
                                                      gTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCgcgcgcc                                                                                                <  1:61363/150‑2 (MQ=255)
                                                      gTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGca                                                                                                <  2:5867713/150‑1 (MQ=255)
                                                           aGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACt                                                                                           >  1:4218905/1‑150 (MQ=255)
                                                              cGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAcc                                                                                        <  1:3199962/150‑1 (MQ=255)
                                                               gATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCa                                                                                       <  1:2464534/150‑1 (MQ=255)
                                                                     ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACg                                                                                 >  1:5571046/1‑150 (MQ=255)
                                                                     ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACg                                                                                 >  1:956038/1‑150 (MQ=255)
                                                                     ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACg                                                                                 >  1:6960377/1‑150 (MQ=255)
                                                                     ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACg                                                                                 >  1:634608/1‑150 (MQ=255)
                                                                     ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACg                                                                                 >  2:5456412/1‑150 (MQ=255)
                                                                      tCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                       <  2:2621694/111‑1 (MQ=255)
                                                                      tCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACgg                                                                                >  2:3703831/1‑150 (MQ=255)
                                                                        ccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAt                                                                              >  2:5619433/1‑150 (MQ=255)
                                                                        ccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAt                                                                              <  2:628447/150‑1 (MQ=255)
                                                                           tAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCaa                                                                           <  2:4124203/150‑1 (MQ=255)
                                                                            aGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAg                                                                          <  2:825244/150‑1 (MQ=255)
                                                                            aGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAg                                                                          <  2:938070/150‑1 (MQ=255)
                                                                              ccGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGca                                                                        <  2:2148487/150‑1 (MQ=255)
                                                                              ccGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGca                                                                        <  1:5038520/150‑1 (MQ=255)
                                                                                gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGcaca                                                                      >  2:5372039/1‑150 (MQ=255)
                                                                                  aGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACaaa                                                                    <  1:728586/150‑1 (MQ=255)
                                                                                   gAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAc                                                                   <  1:4763883/150‑1 (MQ=255)
                                                                                      ccTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCtt                                                                >  2:5177058/1‑150 (MQ=255)
                                                                                       ctctCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                               <  2:2048997/150‑1 (MQ=255)
                                                                                        tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                              <  1:3965681/150‑1 (MQ=255)
                                                                                        tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                              <  2:334376/150‑1 (MQ=255)
                                                                                        tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCAACTTGGTGGGCTGCCAGAGAACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAATCCTTcc                                                              >  1:3303285/1‑150 (MQ=255)
                                                                                          tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcccc                                                            >  2:5560668/1‑150 (MQ=255)
                                                                                           cATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCg                                                           <  1:6566257/150‑1 (MQ=255)
                                                                                              tCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgctg                                                        <  1:2981752/150‑1 (MQ=255)
                                                                                                  tCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTgg                                                    >  1:4889870/1‑150 (MQ=255)
                                                                                                       cTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCGGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAg                                               <  2:7396566/150‑1 (MQ=255)
                                                                                                         gACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTTGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGTTAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgcc                                             <  2:4935110/150‑1 (MQ=255)
                                                                                                            gcgcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                          >  2:4573946/1‑150 (MQ=255)
                                                                                                              gcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggccg                                        <  2:6578581/150‑1 (MQ=255)
                                                                                                                  ggCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACa                                    <  1:2015953/150‑1 (MQ=255)
                                                                                                                      cAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgcg                                >  2:2333068/1‑150 (MQ=255)
                                                                                                                        aaTCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgcc                              >  2:3271062/1‑150 (MQ=255)
                                                                                                                         aTCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgccg                             >  1:4629522/1‑150 (MQ=255)
                                                                                                                               gcgcTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGAt                       <  2:677850/150‑1 (MQ=255)
                                                                                                                                 gcTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCGTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGt                     >  2:3825192/1‑150 (MQ=255)
                                                                                                                                   tGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCt                   <  2:5276496/150‑1 (MQ=255)
                                                                                                                                   tGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCt                   <  1:4435215/150‑1 (MQ=255)
                                                                                                                                    ggCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTg                  <  1:1081163/150‑1 (MQ=255)
                                                                                                                                     gCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg                 <  2:2907798/150‑1 (MQ=255)
                                                                                                                                       tcTGGCAGCCCACCAAGGTGGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               >  1:2131620/2‑150 (MQ=255)
                                                                                                                                       tcTGGCAGCCCACCAAGGTGGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               >  1:7316677/2‑150 (MQ=255)
                                                                                                                                       tcTGGCAGCCCACCAAGGTGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               >  2:3447640/2‑150 (MQ=255)
                                                                                                                                       ccTGGCAGCCACCTTGGTGGGCTGCCGGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               >  1:2730985/1‑150 (MQ=255)
                                                                                                                                       ccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACCGGCGCCGCAGGATGTCTGGGa               <  2:2219085/150‑1 (MQ=255)
                                                                                                                                       ccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               >  1:2914137/1‑150 (MQ=255)
                                                                                                                                       ccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               >  1:6025885/1‑150 (MQ=255)
                                                                                                                                       ccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               <  2:7357167/150‑1 (MQ=255)
                                                                                                                                         tGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgc             <  1:2826139/150‑1 (MQ=255)
                                                                                                                                           gCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgcgc           >  2:7301139/1‑150 (MQ=255)
                                                                                                                                                 cACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcg      >  2:6005225/1‑150 (MQ=255)
                                                                                                                                                 cACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcg      >  1:1926758/1‑150 (MQ=255)
                                                                                                                                                  aCCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcgc     >  2:4398088/1‑150 (MQ=255)
                                                                                                                                                    cTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCa   >  1:3600689/1‑150 (MQ=255)
                                                                                                                                                     ttGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAg  >  1:5365638/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAG  >  NC_002947/4686642‑4686939

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: