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

ACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681649‑4681941
                                                                                                                                                   |                                                                                                                                                  
aCCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCCCCCGGATACCTGTATCAGCCGATCTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCtt                                                                                                                                                   <  5:329565/148‑1 (MQ=255)
 ccAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtgg                                                                                                                                                 >  2:375828/1‑149 (MQ=255)
 ccAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtg                                                                                                                                                  <  3:313050/148‑1 (MQ=255)
  cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggt                                                                                                                                                <  8:143007/149‑1 (MQ=255)
  cAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggt                                                                                                                                                >  3:105193/1‑149 (MQ=255)
   aGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtggtg                                                                                                                                               <  4:91917/149‑1 (MQ=255)
     gCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCATGGTggg                                                                                                                                             >  5:55855/1‑149 (MQ=255)
         ggAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGc                                                                                                                                         <  4:361930/149‑1 (MQ=255)
          gAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGcc                                                                                                                                        <  6:55855/149‑1 (MQ=255)
           aGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCa                                                                                                                                       <  5:169607/149‑1 (MQ=255)
            gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGTGGCTGCCag                                                                                                                                      <  8:254880/149‑1 (MQ=255)
             cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCa                                                                                                                                       >  2:208119/1‑147 (MQ=255)
              ggTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCagag                                                                                                                                    >  7:250795/1‑149 (MQ=255)
               gTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGTGGCTGCCAGGGc                                                                                                                                   <  7:137649/149‑1 (MQ=255)
               gTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtg                                                                                                                                                  <  8:284254/134‑1 (MQ=255)
               gTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTtg                                                                                                                                                  >  7:284254/1‑134 (MQ=255)
                      gAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                               <  6:171202/146‑1 (MQ=255)
                      gAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                               <  1:408857/146‑1 (MQ=255)
                         tCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGc                                                                                                                         >  5:126524/1‑149 (MQ=255)
                             cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                                     >  4:233114/1‑149 (MQ=255)
                              gAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    >  5:425818/1‑149 (MQ=255)
                              gAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    >  1:335782/1‑149 (MQ=255)
                                  tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                                >  6:87759/1‑149 (MQ=255)
                                     ttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATcc                                                                                                              <  1:375828/148‑1 (MQ=255)
                                       cGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGa                                                                                                            >  3:396642/1‑148 (MQ=255)
                                        gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGATCCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGAt                                                                                                           <  1:187746/148‑1 (MQ=255)
                                        gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTGGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATa                                                                                                          <  3:209425/149‑1 (MQ=255)
                                             gATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGTGGCTGCCAGGGcca                                                                                                                                 <  7:140358/121‑3 (MQ=37)
                                             gATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGTGGCTGCCAGGGcca                                                                                                                                 >  8:140358/1‑119 (MQ=37)
                                             gATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGtt                                                                                                     >  5:424163/1‑149 (MQ=255)
                                                  cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGCCGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                >  3:85336/1‑149 (MQ=255)
                                                  cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                >  1:313090/1‑149 (MQ=255)
                                                  cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                <  3:374731/149‑1 (MQ=255)
                                                   tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgcg                                                                                               >  7:270276/1‑149 (MQ=255)
                                                       tCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcac                                                                                            <  2:335782/148‑1 (MQ=255)
                                                         aGCCGATGTCTTCCTGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACaa                                                                                          <  4:78578/148‑1 (MQ=255)
                                                              aTGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCa                                                                                    <  8:270276/149‑1 (MQ=255)
                                                                gTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGtt                                                                                                     <  2:153408/130‑1 (MQ=255)
                                                                gTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGtt                                                                                                     >  1:153408/1‑130 (MQ=255)
                                                                  cTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCa                                                                                >  3:79152/1‑149 (MQ=255)
                                                                   ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGCGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                               >  7:189005/1‑149 (MQ=255)
                                                                   ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                               >  3:106682/1‑149 (MQ=255)
                                                                           gCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCaa                                                                        <  3:414894/148‑1 (MQ=255)
                                                                           gCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAg                                                                       >  5:285474/1‑149 (MQ=255)
                                                                             cGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGca                                                                     <  2:313090/149‑1 (MQ=255)
                                                                              gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGca                                                                     >  3:398622/1‑148 (MQ=255)
                                                                               cAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGcac                                                                    >  4:189963/1‑148 (MQ=255)
                                                                                     ctctCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCtt                                                             <  6:425818/149‑1 (MQ=255)
                                                                                      tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCCTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                            >  3:339714/1‑149 (MQ=255)
                                                                                        tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                            <  6:126524/147‑1 (MQ=255)
                                                                                           ttCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgc                                                       >  7:190230/1‑149 (MQ=255)
                                                                                                   cccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGTACAAACCTTCCCCGCTGCTGGTc                                               <  1:208119/149‑1 (MQ=255)
                                                                                                        aCGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGGTGACAGAGCACCCGTCAGCTCGACGAATCCGATTGCGGTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgcc                                          <  3:112169/149‑1 (MQ=255)
                                                                                                                 gCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGaa                                   <  6:424163/147‑1 (MQ=255)
                                                                                                                 gCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACa                                 <  4:79152/149‑1 (MQ=255)
                                                                                                                    cAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgc                              <  4:85336/149‑1 (MQ=255)
                                                                                                                    cAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgc                              >  4:82907/1‑149 (MQ=255)
                                                                                                                               gcTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGAt                    >  6:33739/1‑148 (MQ=255)
                                                                                                                               gcTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATg                   <  3:233114/149‑1 (MQ=255)
                                                                                                                                cTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                >  3:158438/1‑71 (MQ=255)
                                                                                                                                cTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                <  4:158438/71‑1 (MQ=255)
                                                                                                                                 tGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTc                 >  3:265084/1‑149 (MQ=255)
                                                                                                                                   gCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTg               >  6:93065/1‑149 (MQ=255)
                                                                                                                                    cccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              >  1:7075/1‑149 (MQ=255)
                                                                                                                                     ccTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCCCGGATCAAGCACAAACCTTCCCCGCTGCTGGTCACGCCGGCCGCACAGGCGCCGCAGGATGTCTgg              >  4:280509/1‑148 (MQ=255)
                                                                                                                                     ccTGGCAGCCACCTTGGTGGGCTGCAAGAGCACACGTAAGCTCTACGAATCCGATAGCGTTCGCGCGCCCAACTACCAGGCACGGATCAAGCACAGACCTTACCCGCTGCTGGTCAAGCCGGCCGCACAGGCTCAGCAGGATGTCtag              >  5:440313/1‑146 (MQ=255)
                                                                                                                                      cTGGCAGCCCACCAAGGTGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              >  7:133757/1‑147 (MQ=255)
                                                                                                                                         gCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAg           <  5:33739/147‑1 (MQ=255)
                                                                                                                                              ccACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              >  3:26621/1‑140 (MQ=255)
                                                                                                                                              ccACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              <  4:26621/140‑1 (MQ=255)
                                                                                                                                                aCCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcg   >  7:209819/1‑149 (MQ=255)
                                                                                                                                                  cTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcgc  >  1:416846/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681649‑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: