Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,681,795 T→A 5.6% L29M (TTG→ATG) ‡ 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,681,7950TA 5.6% 220.2 / 2.4 71L29M (TTG→ATG) ‡mltDmembrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

GACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681648‑4681941
                                                                                                                                                     |                                                                                                                                                  
gACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCtt                                                                                                                                                   >  8:279734/1‑149 (MQ=255)
 aCCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtg                                                                                                                                                  >  8:312063/1‑149 (MQ=255)
 aCCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtg                                                                                                                                                  <  4:241346/149‑1 (MQ=255)
   cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACGTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggt                                                                                                                                                >  6:206783/1‑149 (MQ=255)
        ccGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCt                                                                                                                                           <  6:250100/149‑1 (MQ=255)
        ccGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCt                                                                                                                                           >  4:260378/1‑149 (MQ=255)
          ggAGAGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGc                                                                                                                                         <  4:35133/149‑1 (MQ=255)
           gAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGcc                                                                                                                                        <  7:116823/149‑1 (MQ=255)
             gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGTGGCTGCCag                                                                                                                                      <  1:239591/149‑1 (MQ=255)
              cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCGGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                     <  6:34384/149‑1 (MQ=255)
              cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                     <  7:279734/149‑1 (MQ=255)
              cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                     <  7:18131/149‑1 (MQ=255)
                    cTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAcgc                                                                                                                               >  5:324006/1‑147 (MQ=255)
                    cTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                               >  1:113888/1‑149 (MQ=255)
                        aaTCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCa                                                                                                                           <  5:206783/149‑1 (MQ=255)
                            gCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTc                                                                                                                       >  8:281825/1‑149 (MQ=255)
                                aaCTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    >  4:166257/1‑148 (MQ=255)
                                  cTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGaa                                                                                                                 <  2:87616/149‑1 (MQ=255)
                                        cGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGa                                                                                                            <  7:312063/148‑1 (MQ=255)
                                         gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATa                                                                                                          <  4:302454/149‑1 (MQ=255)
                                                   cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGc                                                                                                                                         <  3:32602/108‑1 (MQ=255)
                                                   cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGc                                                                                                                                         >  4:32602/1‑108 (MQ=255)
                                                   cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCAGTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                >  3:126700/1‑149 (MQ=255)
                                                        tCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcac                                                                                            >  3:38169/1‑148 (MQ=255)
                                                        tCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcac                                                                                            >  8:283919/1‑148 (MQ=255)
                                                               aTGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCa                                                                                    <  6:324006/149‑1 (MQ=255)
                                                               aTGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCa                                                                                    <  2:42015/149‑1 (MQ=255)
                                                                  tCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTACCCGAATCCGATAGCGTTCGCGCGCACCACTACCAgg                                                                                  >  2:360576/1‑148 (MQ=255)
                                                                    ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                               >  3:131075/1‑149 (MQ=255)
                                                                    ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                               >  3:35281/1‑149 (MQ=255)
                                                                      cccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                               <  2:225468/147‑1 (MQ=255)
                                                                      cccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACg                                                                              <  7:283919/148‑1 (MQ=255)
                                                                            gCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAg                                                                       <  7:281825/149‑1 (MQ=255)
                                                                             ccGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGc                                                                      <  2:113888/149‑1 (MQ=255)
                                                                               gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCAGGGCAGCCACCTTGGTGGGCAGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGca                                                                     >  1:105408/1‑149 (MQ=255)
                                                                                     ccTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAcc                                                               >  8:98697/1‑148 (MQ=255)
                                                                                       tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                            >  6:97886/1‑149 (MQ=255)
                                                                                        ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                           >  1:316280/1‑149 (MQ=255)
                                                                                         tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGTGGCTGCCagg                                                                                                                                     <  6:9900/74‑2 (MQ=25)
                                                                                         tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGTGGCTGCCagg                                                                                                                                     >  5:9900/1‑73 (MQ=25)
                                                                                         tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTccc                                                          >  1:342340/1‑149 (MQ=255)
                                                                                         tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                            <  3:166257/147‑1 (MQ=255)
                                                                                           aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                     >  5:322796/1‑72 (MQ=255)
                                                                                           aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCaga                                                                                                                                     <  6:322796/72‑1 (MQ=255)
                                                                                           aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTccc                                                          <  5:97886/147‑1 (MQ=255)
                                                                                           aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                           <  1:219924/146‑1 (MQ=255)
                                                                                            ttCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgc                                                       >  3:321784/1‑149 (MQ=255)
                                                                                             tCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgct                                                      >  8:313580/1‑149 (MQ=255)
                                                                                                gTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgctgct                                                   >  3:60299/1‑149 (MQ=255)
                                                                                                   gCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTgg                                                 <  8:260085/148‑1 (MQ=255)
                                                                                                    cccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTc                                                                                                    <  4:337049/96‑1 (MQ=255)
                                                                                                    cccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTc                                                                                                    >  3:337049/1‑96 (MQ=255)
                                                                                                     ccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTc                                               >  6:341459/1‑148 (MQ=255)
                                                                                                       tGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGTCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAg                                            >  8:211510/1‑149 (MQ=255)
                                                                                                            cgcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                       <  4:131075/149‑1 (MQ=255)
                                                                                                            cgcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                       <  4:321784/149‑1 (MQ=255)
                                                                                                            cgcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCAAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                       <  7:312531/149‑1 (MQ=255)
                                                                                                                   cccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAg                                <  4:35281/149‑1 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCATGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACa                                 >  2:315866/1‑147 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCATGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACa                                 <  1:315866/147‑1 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCACGCACAAACCTTCCCCGCTGCTGGTCCAGCCGGCCGAACAgg                               >  8:54019/1‑149 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                               >  1:59074/1‑149 (MQ=255)
                                                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAACTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                               >  1:102248/1‑149 (MQ=255)
                                                                                                                      aaaTCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgcg                             >  1:200167/1‑149 (MQ=255)
                                                                                                                            gTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAg                       <  4:126700/149‑1 (MQ=255)
                                                                                                                               cgcTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGAt                    <  3:260378/149‑1 (MQ=255)
                                                                                                                                     ctcTGGCAGCCCACCAAGGTGGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              >  6:325276/3‑149 (MQ=255)
                                                                                                                                     cccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              <  3:48831/149‑1 (MQ=255)
                                                                                                                                       cTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa            <  2:105408/149‑1 (MQ=255)
                                                                                                                                                 cACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgc    >  4:210722/1‑149 (MQ=255)
                                                                                                                                                  aCCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcg   >  5:47984/1‑149 (MQ=255)
                                                                                                                                                   ccTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcgc  >  7:254233/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681648‑4681941

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: