Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,683,764 T→A 5.6% 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,683,7640TA 5.6% 587.8 / 12.7 181L29* (TTG→TAG) ‡mltDmembrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base T (83/87);  new base A (5/6);  total (88/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

GACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAG  >  NC_002947/4683616‑4683912
                                                                                                                                                     |                                                                                                                                                    
gACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtg                                                                                                                                                     <  1:1008290/150‑1 (MQ=255)
 aCCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtgg                                                                                                                                                    >  2:3742634/1‑150 (MQ=255)
 aCCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtgg                                                                                                                                                    <  2:3570285/150‑1 (MQ=255)
  ccAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggt                                                                                                                                                   >  1:475918/1‑150 (MQ=255)
   cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                                  >  2:3467609/1‑150 (MQ=255)
   cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                                  <  2:2606264/150‑1 (MQ=255)
   cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                                  <  2:1134834/150‑1 (MQ=255)
    aGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtgg                                                                                                                                                 >  2:3969685/1‑150 (MQ=255)
       cccGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCt                                                                                                                                              >  1:2960669/1‑150 (MQ=255)
         cGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGc                                                                                                                                            <  2:4192658/150‑1 (MQ=255)
         cGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGc                                                                                                                                            <  2:904569/150‑1 (MQ=255)
          ggAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCAGTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGcc                                                                                                                                           <  2:4996993/150‑1 (MQ=255)
           gAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCa                                                                                                                                          <  1:5413063/150‑1 (MQ=255)
            aGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCAAGGTGGCTGCCagg                                                                                                                                        <  1:2314572/150‑2 (MQ=255)
             gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGTGGCTGCCagg                                                                                                                                        <  2:3633164/150‑2 (MQ=255)
             gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                        <  2:106447/150‑1 (MQ=255)
              cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAACTCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCagag                                                                                                                                       <  1:2514442/150‑1 (MQ=255)
              cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGACAGCCACCTTGGTGGGCTGCCagag                                                                                                                                       >  2:2044591/1‑150 (MQ=255)
               ggTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGGGc                                                                                                                                      <  1:3334029/150‑1 (MQ=255)
                     tAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                                >  2:2011680/1‑150 (MQ=255)
                       gAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCa                                                                                                                              <  2:481205/150‑1 (MQ=255)
                             ccGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                                        <  2:2272908/150‑1 (MQ=255)
                             ccGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                                        <  1:662634/150‑1 (MQ=255)
                              cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                       <  1:1685837/150‑1 (MQ=255)
                                  cTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                                   <  1:4821564/150‑1 (MQ=255)
                                   tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATc                                                                                                                  <  1:4984213/150‑1 (MQ=255)
                                   tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATc                                                                                                                  >  1:2467355/1‑150 (MQ=255)
                                   tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATc                                                                                                                  <  2:2745648/150‑1 (MQ=255)
                                         gCCCGGATACCTGTATCAGCCGATGTCTTCCCGAAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAg                                                                                                            >  2:4903843/1‑150 (MQ=255)
                                          cccGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGc                                                                                                           <  1:823398/150‑1 (MQ=255)
                                           ccGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGTGGCTGCCa                                                                                                                                          <  1:3243273/118‑1 (MQ=255)
                                           ccGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGTGGCTGCCa                                                                                                                                          >  2:3243273/1‑118 (MQ=255)
                                                  ccTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                   <  2:280971/150‑1 (MQ=255)
                                                   cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgcg                                                                                                  >  1:1852163/1‑150 (MQ=255)
                                                    tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCgcgcgc                                                                                                 >  2:1714121/1‑150 (MQ=255)
                                                    tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCgcgcgc                                                                                                 >  2:5157451/1‑150 (MQ=255)
                                                     gTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGca                                                                                                <  2:2627045/150‑1 (MQ=255)
                                                     gTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGca                                                                                                <  1:4430220/150‑1 (MQ=255)
                                                      tATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggt                                                                                                                                                   <  2:2065445/98‑1 (MQ=255)
                                                      tATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggt                                                                                                                                                   >  1:2065445/1‑98 (MQ=255)
                                                      tATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcac                                                                                               >  1:3321146/1‑150 (MQ=255)
                                                       aTCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCt                                                                                                                                              >  2:2976931/1‑102 (MQ=255)
                                                       aTCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCt                                                                                                                                              <  1:2976931/102‑1 (MQ=255)
                                                       aTCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcaca                                                                                              >  2:501866/1‑150 (MQ=255)
                                                          aGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACt                                                                                           >  2:478617/1‑150 (MQ=255)
                                                          aGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACt                                                                                           >  2:266873/1‑150 (MQ=255)
                                                            ccGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAc                                                                                         >  1:3870173/1‑150 (MQ=255)
                                                              gATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCa                                                                                       <  1:505905/150‑1 (MQ=255)
                                                                    ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACg                                                                                 >  2:4832475/1‑150 (MQ=255)
                                                                       ccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                                        <  1:3887077/108‑1 (MQ=255)
                                                                       ccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                                        >  2:3887077/1‑108 (MQ=255)
                                                                       ccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAt                                                                              <  1:4024898/150‑1 (MQ=255)
                                                                         gTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCa                                                                            <  2:246093/150‑1 (MQ=255)
                                                                         gTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCa                                                                            >  2:3594894/1‑150 (MQ=255)
                                                                          tAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCaa                                                                           <  2:313614/150‑1 (MQ=255)
                                                                          tAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCaa                                                                           <  1:3594891/150‑1 (MQ=255)
                                                                            gCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGc                                                                         <  2:3795111/150‑1 (MQ=255)
                                                                             ccGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGTGGCTGCCAGGGcca                                                                                                                                    <  1:4885753/90‑3 (MQ=21)
                                                                             ccGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGca                                                                        >  1:1305835/1‑150 (MQ=255)
                                                                                     ccTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCtt                                                                <  1:3467609/150‑1 (MQ=255)
                                                                                       tctcATTCCGTCGCCCTGCCGCGCCTGGCCCAAATCCGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATCGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                              <  1:414439/150‑1 (MQ=255)
                                                                                        ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCCCAAACCTTccc                                                             >  2:163616/1‑150 (MQ=255)
                                                                                        ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCATCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTccc                                                             >  1:2742687/1‑150 (MQ=255)
                                                                                         tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcccc                                                            >  1:1886188/1‑150 (MQ=255)
                                                                                           aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgc                                                          >  1:3161890/1‑150 (MQ=255)
                                                                                             tCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgctg                                                        >  2:5143174/1‑150 (MQ=255)
                                                                                                gTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGcca                                                                                                                                    >  1:5289664/1‑69 (MQ=255)
                                                                                                gTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGcca                                                                                                                                    <  2:5289664/71‑3 (MQ=255)
                                                                                                gTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGctgctg                                                     <  1:1999170/150‑1 (MQ=255)
                                                                                                gTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGctgctg                                                     <  1:918698/150‑1 (MQ=255)
                                                                                                     ccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACgg                                                                                <  1:128566/118‑1 (MQ=255)
                                                                                                     ccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACgg                                                                                >  2:128566/1‑118 (MQ=255)
                                                                                                     ccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCaa                                                <  2:4981805/150‑1 (MQ=255)
                                                                                                          cgcgCCTGGCCCAACTCACTGCGCTGGCCCTGGCAGCCACCATGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGAGAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccgg                                           <  2:901555/150‑1 (MQ=255)
                                                                                                           gcgcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAcc                                                                                        >  1:262575/1‑104 (MQ=255)
                                                                                                               cTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGaa                                      <  1:4538406/150‑1 (MQ=255)
                                                                                                                   cccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                                  <  1:5157254/150‑1 (MQ=255)
                                                                                                                   cccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                                  <  2:1110691/150‑1 (MQ=255)
                                                                                                                   cccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                                  >  1:2720814/1‑150 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                                     >  2:2541544/1‑82 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                                     <  1:2541544/82‑1 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgc                                 >  2:1353485/1‑150 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgc                                 >  1:5205746/1‑150 (MQ=255)
                                                                                                                     cAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgcg                                >  2:3408514/1‑150 (MQ=255)
                                                                                                                          cAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCa                           >  2:1027592/1‑150 (MQ=255)
                                                                                                                            gTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAgg                         <  1:3742634/150‑1 (MQ=255)
                                                                                                                             tGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGa                        <  2:2934377/150‑1 (MQ=255)
                                                                                                                                gcTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATCGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGt                     <  1:2752476/150‑1 (MQ=255)
                                                                                                                                  tGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCt                   <  2:1925370/150‑1 (MQ=255)
                                                                                                                                  tGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCt                   <  1:1202353/150‑1 (MQ=255)
                                                                                                                                      ccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa               <  2:575724/150‑1 (MQ=255)
                                                                                                                                        tGGCAGCCCACCAAGGTGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTg                  <  1:1068170/145‑1 (MQ=255)
                                                                                                                                        tGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                          >  2:2464094/1‑121 (MQ=255)
                                                                                                                                        tGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                          <  1:2464012/121‑1 (MQ=255)
                                                                                                                                        tGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgc             >  1:5002721/1‑150 (MQ=255)
                                                                                                                                        tGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgc             >  1:1357093/1‑150 (MQ=255)
                                                                                                                                          gCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgcgc           <  1:4331522/150‑1 (MQ=255)
                                                                                                                                               ccACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgc       <  2:475918/150‑1 (MQ=255)
                                                                                                                                                cACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcg      >  2:4429183/1‑150 (MQ=255)
                                                                                                                                                  ccTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCCCGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGcc    >  2:1033739/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAG  >  NC_002947/4683616‑4683912

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: