Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,050,942 T→C 6.8% V583A (GTC→GCC)  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,9420TC 6.8% 223.7 / 4.5 76V583A (GTC→GCC) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (2/3);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAA  >  NC_002947/6050795‑6051090
                                                                                                                                                   |                                                                                                                                                     
gCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCg                                                                                                                                                     >  2:214841/1‑150 (MQ=255)
  aTGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGc                                                                                                                                                   >  1:1781617/1‑150 (MQ=255)
       cTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGCCGGccgaca                                                                                                                                              <  1:1556039/150‑5 (MQ=255)
              gTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggcc                                                                                                                                      >  1:472436/1‑150 (MQ=255)
              gTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggcc                                                                                                                                      <  2:472436/150‑1 (MQ=255)
                   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)
                                          gCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTc                                                                                                                                                      >  2:236286/1‑107 (MQ=255)
                                          gCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTc                                                                                                                                                      <  1:236286/107‑1 (MQ=255)
                                          gCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCCGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGCCGTGCACCtccaga                                                                                                          >  2:344805/1‑145 (MQ=255)
                                              aCAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTg                                                                                                      >  2:2227096/1‑150 (MQ=255)
                                                  cTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTc                                                                                                                                                      >  1:1862144/1‑99 (MQ=255)
                                                  cTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTc                                                                                                                                                      <  2:1862144/99‑1 (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:344805/150‑1 (MQ=255)
                                                           ccAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                         <  1:2227096/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:381890/150‑1 (MQ=255)
                                                              gAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGaa                                                                                      >  2:556521/1‑150 (MQ=255)
                                                               aaGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGaaa                                                                                     <  2:1781617/150‑1 (MQ=255)
                                                                     cACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTc                                                                                                                                                      <  2:2190159/80‑1 (MQ=255)
                                                                     cACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTc                                                                                                                                                      >  1:2190159/1‑80 (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                                                                                                                     >  2:1877394/1‑94 (MQ=255)
                                                                                       ccGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGt                                                                                                                     <  1:1877394/94‑1 (MQ=255)
                                                                                        cGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTc                                                                                                                                          >  2:1998457/1‑72 (MQ=255)
                                                                                        cGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTc                                                                                                                                          <  1:1998457/72‑1 (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                                                                                         <  1:217307/114‑1 (MQ=255)
                                                                                               cTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGt                                                                                         >  2:217307/1‑114 (MQ=255)
                                                                                                      tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCa                                                                                           >  1:1213893/1‑105 (MQ=255)
                                                                                                      tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCa                                                                                           <  2:1213893/105‑1 (MQ=255)
                                                                                                      tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTcc                                                  >  1:358670/1‑146 (MQ=255)
                                                                                                      tCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTcc                                                  <  2:358670/146‑1 (MQ=255)
                                                                                                         gCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                >  1:1606178/1‑113 (MQ=255)
                                                                                                         gCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                <  2:1606178/113‑1 (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                                                              >  2:678436‑M1/11‑102 (MQ=255)
                                                                                                                                      gaccatgcccTTGCCGGCCGACAGGTGGCCGACAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAgtgc                                                              <  1:678436‑M1/92‑1 (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                                                                                <  2:1888364/78‑1 (MQ=255)
                                                                                                                                            ccGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGc                                                                                >  1:1888364/1‑78 (MQ=255)
                                                                                                                                            ccGATTGTCGGCCACCGGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCgt        <  2:2242969/150‑1 (MQ=255)
                                                                                                                                                  GTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGaa  >  1:1848064‑M1/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCGAACTGCGTGTCGAA  >  NC_002947/6050795‑6051090

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: