Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,050,543 C→G 5.5% P450R (CCG→CGG)  spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,050,5430CG 5.5% 343.6 / 8.8 111P450R (CCG→CGG) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base C (44/60);  new base G (2/4);  total (46/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 114 total aligned reads displayed.

GTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTC  >  NC_002947/6050394‑6050688
                                                                                                                                                     |                                                                                                                                                 
gTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGcccc                                                                                                                                                   >  2:2057541/1‑150 (MQ=255)
   gCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCggg                                                                                                                                                <  1:810995/150‑1 (MQ=255)
    cccAAAGGCTCCACGGCCGTCGAATTCGCCTGCGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGgc                                                                                                                                               >  2:1480331/1‑150 (MQ=255)
     ccAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGgcg                                                                                                                                              >  2:1008311/1‑150 (MQ=255)
     ccAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGgcg                                                                                                                                              >  1:290437/1‑150 (MQ=255)
     ccAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGgcg                                                                                                                                              >  1:372193/1‑150 (MQ=255)
      cAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCCCACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGgcgc                                                                                                                                             >  2:1389496/1‑150 (MQ=255)
       aaaGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCt                                                                                                                                            >  1:361769/1‑150 (MQ=255)
        aaGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTc                                                                                                                                           <  1:1894635/150‑1 (MQ=255)
         aGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                          <  2:290437/150‑3 (MQ=255)
         aGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                          <  1:1577030/150‑3 (MQ=255)
         aGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                          <  1:2113470/150‑3 (MQ=255)
         aGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                          <  1:1008311/150‑3 (MQ=255)
         aGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCg                                                                                                                                          <  1:1215360/150‑1 (MQ=255)
           gCTCCACGGCCGTCGACTTCGCCTACGCGGGCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGACCTAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGc                                                                                                                                        <  1:1355195/150‑1 (MQ=255)
               cACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCaaa                                                                                                                                    <  1:517444/150‑1 (MQ=255)
                aCGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAAc                                                                                                                                   >  2:1158244/1‑150 (MQ=255)
                   gCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCg                                                                                                                                >  1:209666/1‑150 (MQ=255)
                     cGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTAGGAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGc                                                                                                                              >  1:1831695/1‑150 (MQ=255)
                       tCGACTTCGCCTACGCGGTCCCCCCCCCCGTCGGTAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCAt                                                                                                                            <  2:1831695/150‑1 (MQ=255)
                       tCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCAt                                                                                                                            <  1:2319953/150‑1 (MQ=255)
                        cGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATg                                                                                                                           <  2:1133117/150‑1 (MQ=255)
                         gACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATgg                                                                                                                          <  2:372193/150‑1 (MQ=255)
                                  tACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGc                                                                                                                         <  2:1650247/142‑1 (MQ=255)
                                  tACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGc                                                                                                                         >  1:1650247/1‑142 (MQ=255)
                                   aCGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtg                                                                                                                <  2:361769/150‑1 (MQ=255)
                                    cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                               <  1:1724508/150‑1 (MQ=255)
                                    cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                               <  1:2121057/150‑1 (MQ=255)
                                    cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                               <  1:1158244/150‑1 (MQ=255)
                                    cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                               <  1:569231/150‑1 (MQ=255)
                                    cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCGGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                               <  2:515465/150‑1 (MQ=255)
                                     gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                              <  2:1016819/150‑1 (MQ=255)
                                     gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                              <  2:2041327/150‑1 (MQ=255)
                                     gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                              <  1:1714777/150‑1 (MQ=255)
                                     gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                              <  1:2323512/150‑1 (MQ=255)
                                     gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                              <  2:2314522/150‑1 (MQ=255)
                                         tCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGtct                                                                                                          >  2:1350316/1‑150 (MQ=255)
                                          ccACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGtctc                                                                                                         <  1:305254/150‑1 (MQ=255)
                                                cGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCaa                                                                                                   <  1:1480331/150‑1 (MQ=255)
                                                cGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCaa                                                                                                   <  1:1788749/150‑1 (MQ=255)
                                                 gACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAg                                                                                                  <  2:2038084/150‑1 (MQ=255)
                                                    gTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgca                                                                                               <  1:76540/150‑1 (MQ=255)
                                                        gCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATg                                                                                                                           >  2:1406115/1‑118 (MQ=255)
                                                        gCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATg                                                                                                                           <  1:1406115/118‑1 (MQ=255)
                                                        gCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgcacgca                                                                                           >  1:117473/1‑150 (MQ=255)
                                                         cAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCCGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCACACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgcacgcac                                                                                          <  2:1747383/150‑1 (MQ=255)
                                                          aaCAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtg                                                                                                                <  1:1698469/127‑1 (MQ=255)
                                                          aaCAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtg                                                                                                                >  2:1698469/1‑127 (MQ=255)
                                                              gTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAata                                                                                     >  1:1254499/1‑150 (MQ=255)
                                                                tGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGGACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTTCTTTCTCGGGCAAGGCACGCACGAATATc                                                                                   <  1:64681/150‑1 (MQ=255)
                                                                 gCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAAc                                                                                                                    <  1:70065/116‑1 (MQ=255)
                                                                 gCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAAc                                                                                                                    >  2:70065/1‑116 (MQ=255)
                                                                   aTCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGc                                                                                <  2:209666/150‑1 (MQ=255)
                                                                   aTCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGc                                                                                <  1:642127/150‑1 (MQ=255)
                                                                        ttGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCg                                                                                                                                          <  1:1060683/87‑1 (MQ=255)
                                                                        ttGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCg                                                                                                                                          >  2:1060683/1‑87 (MQ=255)
                                                                         tGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTGAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCCCGCATATCCGCCAcgcc                                                                          >  2:1939790/1‑149 (MQ=255)
                                                                                  aaCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCaa                                                                 >  2:1331933/1‑150 (MQ=255)
                                                                                   aCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAAc                                                                >  1:1406563/1‑150 (MQ=255)
                                                                                     cgccgcCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCCCAACTTTGTGGTCTCGGGCACGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAg                                                              >  2:410796/1‑150 (MQ=255)
                                                                                       ccgccTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAgcg                                                            <  1:1567370/150‑1 (MQ=255)
                                                                                       ccgccTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAgcg                                                            >  2:1702822/1‑150 (MQ=255)
                                                                                       ccgccTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAgcg                                                            >  1:312078/1‑150 (MQ=255)
                                                                                        cgccTGGCGCCGCTCTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAgcgc                                                           <  1:2263201/150‑1 (MQ=255)
                                                                                         gccTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGcgcc                                                          >  1:503510/1‑150 (MQ=255)
                                                                                                gccgcTGTCCGCACCGATACAAAGCGGCTCGACCGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGa                                                   <  2:127707/150‑1 (MQ=255)
                                                                                                gccgcTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGa                                                   <  2:503510/150‑1 (MQ=255)
                                                                                                gccgaTGTCCGATCCGCTACAAAGCGGCTGGACAGTGGAAATCGTCAGCGCCCTGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAGGGCACGCACGCATATCCGCTACGCGCTCAAGCAACAGCGCCGCTCGGa                                                   <  2:611866/150‑1 (MQ=255)
                                                                                                   gcTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTc                                                <  1:2057541/150‑1 (MQ=255)
                                                                                                    cTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCg                                                                                                                                >  2:242021/1‑69 (MQ=255)
                                                                                                    cTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCg                                                                                                                                <  1:242021/69‑1 (MQ=255)
                                                                                                    cTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTcc                                               <  1:2045873/150‑1 (MQ=255)
                                                                                                       tCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATc                                            >  1:1776795/1‑150 (MQ=255)
                                                                                                          gAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGc                                         >  2:71295/1‑150 (MQ=255)
                                                                                                          gAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGc                                         >  2:808695/1‑150 (MQ=255)
                                                                                                            aCCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCt                                       >  1:2285762/1‑150 (MQ=255)
                                                                                                                cTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCa                                              <  2:747253/139‑1 (MQ=255)
                                                                                                                cTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCa                                              >  1:747253/1‑139 (MQ=255)
                                                                                                                   cAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAg                                <  2:165195/150‑1 (MQ=255)
                                                                                                                          gCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgcac                                                                                              <  1:1679432/81‑1 (MQ=255)
                                                                                                                           cTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAgcgc                                                           >  1:1613911/1‑115 (MQ=255)
                                                                                                                           cTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAgcgc                                                           <  2:1613911/115‑1 (MQ=255)
                                                                                                                           cTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCctgct                        <  1:1331933/150‑1 (MQ=255)
                                                                                                                                cAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACa                   >  2:1556039/1‑150 (MQ=255)
                                                                                                                                 aGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACaa                  <  2:1853472/150‑1 (MQ=255)
                                                                                                                                   tGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAgg                >  2:55725/1‑150 (MQ=255)
                                                                                                                                   tGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAgg                <  1:98676/150‑1 (MQ=255)
                                                                                                                                   cgaaaaTCGCCAGCCCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAgg                <  1:2211651/147‑1 (MQ=255)
                                                                                                                                     gAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGa                                 <  2:99219/131‑1 (MQ=255)
                                                                                                                                     gAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGa                                 >  1:99219/1‑131 (MQ=255)
                                                                                                                                           gTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCctgc                         <  1:962119/133‑1 (MQ=255)
                                                                                                                                           gTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCctgc                         >  2:962119/1‑133 (MQ=255)
                                                                                                                                           gTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACt        >  2:1532974/1‑150 (MQ=255)
                                                                                                                                              aGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGc                                                                                >  2:397285/1‑75 (MQ=255)
                                                                                                                                              aGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGc                                                                                <  1:397285/75‑1 (MQ=255)
                                                                                                                                              aGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGc     >  1:202064/1‑150 (MQ=255)
                                                                                                                                              aGCGCCCCGGGCGCTCGGCCAAACCCGGCACGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGc     >  2:1606302/1‑150 (MQ=255)
                                                                                                                                               gcgcCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCt    >  2:1379567/1‑150 (MQ=255)
                                                                                                                                               gcgcCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCt    >  2:1378710/1‑150 (MQ=255)
                                                                                                                                                cgcCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCtt   <  2:1776795/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTC  >  NC_002947/6050394‑6050688

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