Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,573,517 C→T 5.4% R184C (CGC→TGC)  rsgA → ribosome biogenesis GTPase RsgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,573,5170CT 5.4% 220.4 / 2.1 74R184C (CGC→TGC) rsgAribosome biogenesis GTPase RsgA
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base T (2/2);  total (33/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 = 5.30e-01

CCGGCGCCAGAGCCGCACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAGTCGTCGCTGGTCAACAGCCTGCTGCCA  >  NC_002947/5573370‑5573663
                                                                                                                                                   |                                                                                                                                                  
ccGGCGCCAGAGCCGCACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACcgc                                                                                                                                                  >  1:939648/1‑150 (MQ=255)
     gcCAGAGCCGCACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACt                                                                                                                                             <  2:1032085/150‑1 (MQ=255)
      ccAGAGCCGCACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTg                                                                                                                                            <  1:806006/150‑1 (MQ=255)
       cAGAGCCGCACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTgg                                                                                                                                           >  2:505039/1‑150 (MQ=255)
        agagCCGCACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTggg                                                                                                                                          >  1:441031/1‑150 (MQ=255)
               cACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCg                                                                                                                                   >  1:1207385/1‑150 (MQ=255)
               cACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCg                                                                                                                                   <  1:1535032/150‑1 (MQ=255)
                aCCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCgc                                                                                                                                  >  1:777362/1‑150 (MQ=255)
                 cccGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCgct                                                                                                                                 >  2:316764/1‑150 (MQ=255)
                  ccGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGAGCTACCCgctg                                                                                                                                >  2:1076804/1‑150 (MQ=255)
                        ccgATCGACCGCTCTCTGGTAGCGGCCGAGCACGCCGGCCCGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGACCCGCCATGCATGCGCTGCTGGAGGTCTTCCTCGAACTGGGCTGCCCGCTGCTGGaa                                                                                                                          <  2:73182/148‑1 (MQ=255)
                        cTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGaa                                                                                                                          >  2:920457/1‑150 (MQ=255)
                        cTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGaa                                                                                                                          <  2:939648/150‑1 (MQ=255)
                          gATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGt                                                                                                                        >  2:60380/1‑150 (MQ=255)
                          gATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGt                                                                                                                        >  1:1471840/1‑150 (MQ=255)
                            tCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGttt                                                                                                                      <  1:898419/150‑1 (MQ=255)
                             cGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTc                                                                                                                     <  1:1129823/150‑1 (MQ=255)
                               aCCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCgg                                                                                                                   >  2:1733006/1‑150 (MQ=255)
                               aCCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCgg                                                                                                                   >  1:1472828/1‑150 (MQ=255)
                               aCCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCgg                                                                                                                   <  1:316764/150‑1 (MQ=255)
                               aCCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCgg                                                                                                                   >  2:96397/1‑150 (MQ=255)
                                  gCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGcac                                                                                                                <  1:505039/150‑1 (MQ=255)
                                    tATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCacac                                                                                                              >  1:1193522/1‑150 (MQ=255)
                                    tATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCacac                                                                                                              <  1:60380/150‑1 (MQ=255)
                                       cTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCAcaccac                                                                                                           <  2:1471840/150‑1 (MQ=255)
                                        tGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACg                                                                                                          <  1:742479/150‑1 (MQ=255)
                                                  cGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGGGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCAt                                                                                                >  1:758667/1‑150 (MQ=255)
                                                  cGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCAt                                                                                                <  1:758943/150‑1 (MQ=255)
                                                   gAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACTAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATg                                                                                               >  2:627773/1‑150 (MQ=255)
                                                   gAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGTGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATg                                                                                               <  1:1396937/150‑1 (MQ=255)
                                                       aCGCCGGACTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGc                                                                                              >  2:1302578/1‑147 (MQ=255)
                                                       aCGCCGGACTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGc                                                                                              <  1:1302578/147‑1 (MQ=255)
                                                         gCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAgcgc                                                                                         <  2:1207385/150‑1 (MQ=255)
                                                           cGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCt                                                                                       <  1:920457/150‑1 (MQ=255)
                                                            ggCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTg                                                                                      >  2:1315505/1‑150 (MQ=255)
                                                             gCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTgc                                                                                     <  2:1472828/150‑1 (MQ=255)
                                                                tGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTACAGGTGCTGCGGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTgcagc                                                                                  <  2:60236/150‑1 (MQ=255)
                                                                   gcCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGcaac                                                                               <  1:96397/150‑1 (MQ=255)
                                                                      cGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCaacaac                                                                            <  2:1707929/150‑1 (MQ=255)
                                                                            tgctgctgAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGATCTACCGCGAAGTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGATGGCATGCAGCGCCTGCAGCAACAACTCGACg                                                                      <  1:1689776/150‑1 (MQ=255)
                                                                               tgctgAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAg                                                                                            >  2:833243/1‑125 (MQ=255)
                                                                               tgctgAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAg                                                                                            <  1:833243/125‑1 (MQ=255)
                                                                                     aCAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGTACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCa                                                             >  2:1055164/1‑150 (MQ=255)
                                                                                     aCAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCa                                                             >  1:1240839/1‑150 (MQ=255)
                                                                                       aaGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACCCCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGc                                                           <  1:731891/150‑1 (MQ=255)
                                                                                       aaGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATTAGc                                                           <  2:1193522/150‑1 (MQ=255)
                                                                                          gCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCgtg                                                        <  1:1315505/150‑1 (MQ=255)
                                                                                           ccGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCgtgt                                                       <  2:1929408/150‑1 (MQ=255)
                                                                                            cGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCgt                                                         <  2:1213745/147‑1 (MQ=255)
                                                                                            cGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATCCGCTGCTTGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCgt                                                         >  1:1213745/1‑147 (MQ=255)
                                                                                              aCCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCa                                                                                             >  2:685024/1‑109 (MQ=255)
                                                                                              aCCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCa                                                                                             <  1:685024/109‑1 (MQ=255)
                                                                                                       acgacgAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCgc                                                                                                                                  <  1:949376/63‑1 (MQ=255)
                                                                                                       acgacgAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCgc                                                                                                                                  >  2:949376/1‑63 (MQ=255)
                                                                                                         gacgaGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACTGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGcc                                                                   <  1:1395260/124‑1 (MQ=255)
                                                                                                         gacgaGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACTGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGcc                                                                   <  1:1395293/124‑1 (MQ=255)
                                                                                                         gacgaGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACTGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGcc                                                                   >  2:1395260/1‑124 (MQ=255)
                                                                                                         gacgaGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACTGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGcc                                                                   >  2:1395293/1‑124 (MQ=255)
                                                                                                          acgaGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCgg                                        <  2:777362/150‑1 (MQ=255)
                                                                                                                aTGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTgg                                  <  2:1240839/150‑1 (MQ=255)
                                                                                                                 tGGCCCTGTTCTTCATGCTCTGGTGGAGACCGAGAGCGACCCGGGCTACCCGATGCTGGAAGTTTCGGCACCCCACCGCGACGGAATGCAGGGCCTGAAGCAACGCCTCGACGGCCATATCAGCGTGTTCCTCGGCCAGTCggacgtcgg                                 <  2:356355/150‑8 (MQ=255)
                                                                                                                 tGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTggg                                 <  1:627773/150‑1 (MQ=255)
                                                                                                                         gCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAgtcgt                         >  1:1578559/1‑150 (MQ=255)
                                                                                                                          ccTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCgg                                        >  1:957563/1‑134 (MQ=255)
                                                                                                                          ccTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCgg                                        <  2:957563/134‑1 (MQ=255)
                                                                                                                           cTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACg                                                                      <  1:1906150/103‑1 (MQ=255)
                                                                                                                           cTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACg                                                                      >  2:1906150/1‑103 (MQ=255)
                                                                                                                                    ctgctgGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCCGATCTCCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAGTCGTCGCTGGTCAAc              <  1:1391949/150‑1 (MQ=255)
                                                                                                                                     tgctgGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTg                                   >  1:1405173/1‑128 (MQ=255)
                                                                                                                                     tgctgGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTg                                   <  2:1405173/128‑1 (MQ=255)
                                                                                                                                     tgctgGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAGTCGTCGCTGGTCAACa             >  1:95499/1‑150 (MQ=255)
                                                                                                                                     tgctgGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAGTCGTCGCTGGTCAACa             >  1:541219/1‑150 (MQ=255)
                                                                                                                                       ctggcagTCTACCACGATCTTGCTAACCCGCTCCTTGACGTTTCTCAAAACCACGGCGACGGCATGCAGCGCCTGTAGCAACAACTCGACGTCGATCTTTACCTGTTCGTCGCCCAGTCGGGAGTAGGCAAGTCGTCTCTGGTCAACAGc           <  2:758667/144‑1 (MQ=255)
                                                                                                                                                tACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAGTCGTCGCTGGTCAACAGCCTGCTGCCa  >  1:1461603/1‑150 (MQ=255)
                                                                                                                                                tACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAGTCGTCGCTGGTCAACAGCCTGCTGCCa  >  2:1913634/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGGCGCCAGAGCCGCACCCGAACCTGATCGACCGCTATCTGGTAGCGGCCGAGCACGCCGGCCTGCGCCCGCTGCTGCTGCTGAACAAGGCCGACCTGATCAACGACGAGAATGGCCCCGGCCTGCATGCGCTGCTGGAGGTCTACCGCGAACTGGGCTACCCGCTGCTGGAAGTTTCGGCACACCACGGCGACGGCATGCAGCGCCTGCAGCAACAACTCGACGGCCATATCAGCGTGTTCGTCGGCCAGTCGGGAGTGGGCAAGTCGTCGCTGGTCAACAGCCTGCTGCCA  >  NC_002947/5573370‑5573663

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: