Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,050,960 G→A 5.7% G589D (GGC→GAC)  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,9600GA 5.7% 206.1 / 2.5 70G589D (GGC→GAC) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base A (2/2);  total (32/38)
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.46e-01

GCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAACTGGAACACCAGCGCGGGC  >  NC_002947/6050814‑6051109
                                                                                                                                                   |                                                                                                                                                     
gcgcCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCaa                                                                                                                                                    <  2:1501679/150‑1 (MQ=255)
 cgccgcCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAg                                                                                                                                                   <  2:715939/150‑1 (MQ=255)
  gccgccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAgg                                                                                                                                                  >  1:1103183/1‑150 (MQ=255)
    cgccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGc                                                                                                                                                <  2:1103183/150‑1 (MQ=255)
       ctgctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATg                                                                                                                                             <  1:214841/150‑1 (MQ=255)
          ctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTc                                                                                                                                          <  2:1878958/150‑1 (MQ=255)
           tgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCg                                                                                                                                         >  1:406990/1‑150 (MQ=255)
             tCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCaa                                                                                                                                                    <  1:441476/137‑1 (MQ=255)
             tCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCaa                                                                                                                                                    >  2:441476/1‑137 (MQ=255)
                  cGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCt                                                                                                                                  <  1:1124764/150‑1 (MQ=255)
                       gCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCCGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGCCGTGCACCtccaga                                                                                                                             >  2:344805/1‑145 (MQ=255)
                           aCAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTg                                                                                                                         >  2:2227096/1‑150 (MQ=255)
                                 gCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCaa                                                                                                                   <  1:236976/150‑1 (MQ=255)
                                      cgcCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCa                                                                                                              <  1:39050/150‑1 (MQ=255)
                                      cgcCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCAAGCCCAAGCCTCAGGACCCGATTCAGGTCCTCCTGTCGTCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCa                                                                                                              <  1:1498344/150‑1 (MQ=255)
                                        ccAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                                            <  1:2227096/150‑1 (MQ=255)
                                        ccAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                                            <  1:344805/150‑1 (MQ=255)
                                        ccAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                                            >  2:1848064/1‑150 (MQ=255)
                                         cAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGgaga                                                                                                                             >  1:827379/1‑132 (MQ=255)
                                         cAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGgaga                                                                                                                             <  2:827379/132‑1 (MQ=255)
                                           gAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGaa                                                                                                         >  2:556521/1‑150 (MQ=255)
                                           gAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGaa                                                                                                         <  2:381890/150‑1 (MQ=255)
                                            aaGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGaaa                                                                                                        <  2:1781617/150‑1 (MQ=255)
                                                       gCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAAc                                                                                             >  2:1910701/1‑150 (MQ=255)
                                                         gATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAAccc                                                                                           <  2:1632081/150‑1 (MQ=255)
                                                                 gCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAgt                                                                                   <  2:1161380/150‑1 (MQ=255)
                                                                  cACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAgtg                                                                                  >  2:840745/1‑150 (MQ=255)
                                                                    ccGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGt                                                                                                                                        <  1:1877394/94‑1 (MQ=255)
                                                                    ccGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGt                                                                                                                                        >  2:1877394/1‑94 (MQ=255)
                                                                     cGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGc                                                                                                                     <  1:2058474/112‑1 (MQ=255)
                                                                     cGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGc                                                                                                                     >  2:2058474/1‑112 (MQ=255)
                                                                       aaGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAAc                                                                                                                  <  1:574519/113‑1 (MQ=255)
                                                                       aaGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAAc                                                                                                                  >  2:574519/1‑113 (MQ=255)
                                                                       aaGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAgtgcgtgc                                                                             >  1:2078212/1‑150 (MQ=255)
                                                                          gCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGCAACGTGCGTGCacc                                                                          >  2:1532965/1‑148 (MQ=255)
                                                                            cTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                                            >  2:217307/1‑114 (MQ=255)
                                                                            cTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                                            <  1:217307/114‑1 (MQ=255)
                                                                                   tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCa                                                                                                              <  2:1213893/105‑1 (MQ=255)
                                                                                   tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCa                                                                                                              >  1:1213893/1‑105 (MQ=255)
                                                                                   tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTcc                                                                     <  2:358670/146‑1 (MQ=255)
                                                                                   tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTcc                                                                     >  1:358670/1‑146 (MQ=255)
                                                                                      gCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                                   <  2:1606178/113‑1 (MQ=255)
                                                                                      gCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                                   >  1:1606178/1‑113 (MQ=255)
                                                                                        tACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAg                                                            <  1:913763/150‑1 (MQ=255)
                                                                                           gCCAGGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGAc                                                         >  1:2252554/1‑150 (MQ=255)
                                                                                            ccAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACa                                                        >  1:2012271/1‑150 (MQ=255)
                                                                                                    tgcACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGc                                                >  1:861736/1‑150 (MQ=255)
                                                                                                    tgcACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGc                                                >  1:2324456/1‑150 (MQ=255)
                                                                                                      cACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGa                                              <  2:861736/150‑1 (MQ=255)
                                                                                                       ctccaggtgcacgaccatgcccTTGCCGGCCGACAGGTGGCCGACAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                                   <  1:1738537‑M1/74‑1 (MQ=255)
                                                                                                       ctccaggtgcacgaccatgcccTTGCCGGCCGACAGGTGGCCGACAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                                   >  2:1738537‑M1/23‑96 (MQ=255)
                                                                                                             aTCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAAt                                       <  2:1437171/150‑1 (MQ=255)
                                                                                                             aTCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAAAATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGCGCCAAGGACATCACTGGCGAGttaact                                       <  2:406990/150‑5 (MQ=255)
                                                                                                                   gaccatgcccTTGCCGGCCGACAGGTGGCCGACAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAgtgc                                                                                 <  1:678436‑M1/92‑1 (MQ=255)
                                                                                                                   gaccatgcccTTGCCGGCCGACAGGTGGCCGACAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAgtgc                                                                                 >  2:678436‑M1/11‑102 (MQ=255)
                                                                                                                     tGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACt                               <  2:2252554/150‑1 (MQ=255)
                                                                                                                       aCCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGc                             <  1:840745/150‑1 (MQ=255)
                                                                                                                         ccGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                                            <  1:671843/69‑1 (MQ=255)
                                                                                                                         ccGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                                            >  2:671843/1‑69 (MQ=255)
                                                                                                                         ccGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                                   >  1:1888364/1‑78 (MQ=255)
                                                                                                                         ccGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                                   <  2:1888364/78‑1 (MQ=255)
                                                                                                                         ccGATTGTCGGCCACCGGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCgt                           <  2:2242969/150‑1 (MQ=255)
                                                                                                                               GTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGaa                     >  1:1848064‑M1/1‑150 (MQ=255)
                                                                                                                                  ggCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAACTg                  >  2:2013494/1‑150 (MQ=255)
                                                                                                                                     cACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGaaaa                                                                                     <  1:1840540/80‑1 (MQ=255)
                                                                                                                                     cACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGaaaa                                                                                     >  2:1840540/1‑80 (MQ=255)
                                                                                                                                       aCCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGtt                                           <  2:2303013/121‑1 (MQ=255)
                                                                                                                                       aCCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTAGGCCAAGGACATCACTGGCGAGtt                                           >  1:2303013/1‑121 (MQ=255)
                                                                                                                                           gTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAACTGGa                <  1:1681543/144‑1 (MQ=255)
                                                                                                                                           gTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAACTGGa                >  2:1681543/1‑144 (MQ=255)
                                                                                                                                                  ggcAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAACTGGAACACCAGCGCggg   <  2:2078212/150‑1 (MQ=255)
                                                                                                                                                   gcAAGGGAACAGTCGTGCACCTGGAGAACTGCCGCAACCTCAGTGAAAGCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAACTGGAACACCAGCGCGGGc  <  1:1910701/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAACTGGAACACCAGCGCGGGC  >  NC_002947/6050814‑6051109

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