Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,697,369:1 +A 5.8% intergenic (+23/‑54) hslV → / → hslU ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,697,3691.A 5.8% 275.0 / 10.8 69intergenic (+23/‑54)hslV/hslUATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base . (28/37);  new base A (2/2);  total (30/39)
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.05e-01

CTGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCA  >  NC_002947/5697224‑5697516
                                                                                                                                                  |                                                                                                                                                   
ctgctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCgg                                                                                                                                                 <  2:1751965/150‑1 (MQ=255)
ctgctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCgg                                                                                                                                                 <  1:1310035/150‑1 (MQ=255)
  gctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTg                                                                                                                                               <  2:207284/150‑1 (MQ=255)
        cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgggg                                                                                                                                         <  2:1526977/150‑1 (MQ=255)
        cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgggg                                                                                                                                         >  1:322132/1‑150 (MQ=255)
          aGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                           <  1:1338552/146‑1 (MQ=255)
          aGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                           >  2:1338552/1‑146 (MQ=255)
          aGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCt                                                                                                                                       <  2:187666/150‑1 (MQ=255)
            accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCttt                                                                                                                                     <  1:1464166/150‑1 (MQ=255)
                    gtcGGCCGGCGAAATCGCCGAAACCGCCCTCAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGGAGGGGCTTTTCCACGCt                                                                                                                             <  2:209398/149‑1 (MQ=255)
                            gcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGc                                                                                                                              <  2:4769/141‑1 (MQ=255)
                            gcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGc                                                                                                                              >  1:4769/1‑141 (MQ=255)
                                aaTCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGcc                                                                                                                 >  1:1201803/1‑150 (MQ=255)
                                    gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                           >  2:1632336/1‑120 (MQ=255)
                                    gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                           <  1:1632336/120‑1 (MQ=255)
                                    gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAg                                                                                                             <  2:809028/150‑1 (MQ=255)
                                        agaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGAcc                                                                                                         <  2:632250/150‑1 (MQ=255)
                                           ccGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCATCTGCCCAAGGACCGTa                                                                                                      <  1:1498059/150‑1 (MQ=255)
                                           ccGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTa                                                                                                      <  2:854613/150‑1 (MQ=255)
                                                  gAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGcgcc                                                                                                                                     >  2:1490235/1‑109 (MQ=255)
                                                  gAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGGGAAGGGGcgcc                                                                                                                                     <  1:1490235/112‑4 (MQ=255)
                                                      aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTc                                                                                           >  2:182176/1‑150 (MQ=255)
                                                       tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                           >  1:1051905/1‑101 (MQ=255)
                                                       tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                           <  2:1051905/101‑1 (MQ=255)
                                                         gCCGGCGAGATCTGCGTGTTCAACAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                        <  2:892154/150‑1 (MQ=255)
                                                          ccGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATg                                                                                       >  1:797532/1‑150 (MQ=255)
                                                           cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGa                                                                                      <  1:182176/150‑1 (MQ=255)
                                                           cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGa                                                                                      >  1:1283795/1‑150 (MQ=255)
                                                             gCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTcc                                                                                          >  1:1522133/1‑144 (MQ=255)
                                                             gCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTcc                                                                                          <  2:1522133/144‑1 (MQ=255)
                                                               gACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAcccc                                                                                  >  1:618073/1‑150 (MQ=255)
                                                                aCATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAccccc                                                                                 >  2:174887/1‑150 (MQ=255)
                                                                     tGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCga                                                                            >  1:316916/1‑150 (MQ=255)
                                                                       cGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCgaga                                                                          <  2:797532/150‑1 (MQ=255)
                                                                          gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCAt                                                                                              <  2:1293628/127‑1 (MQ=255)
                                                                          gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCAt                                                                                              >  1:1293628/1‑127 (MQ=255)
                                                                            tCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                     >  1:1680956/1‑150 (MQ=255)
                                                                             cACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTcc                                                                    >  2:1539503/1‑150 (MQ=255)
                                                                               ccaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGCCTCACCTGCCCAACGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCAc                                                                  <  2:322132/150‑1 (MQ=255)
                                                                                 aaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATGGTCCACGa                                                                >  1:331249/1‑150 (MQ=255)
                                                                                           tGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGc                                                      <  2:618073/150‑1 (MQ=255)
                                                                                            gACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGcc                                                     >  1:1846125/1‑150 (MQ=255)
                                                                                               cctCGAGTAGCAGGACCTGGCCGGCTGATTCAGCTGTTCTGGCTCCCCTTCCCGGTGAATGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCTTATTGATCTAATCCAACGCCCCCCGCGAAATCCTACACGAACTCAACCGCcaca                                                  <  2:917899/148‑1 (MQ=255)
                                                                                                aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcat                                                 >  2:1302024/1‑150 (MQ=255)
                                                                                                   gaggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                              <  2:1846125/150‑1 (MQ=255)
                                                                                                     ggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCg                                            <  2:1680956/150‑1 (MQ=255)
                                                                                                        gCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCaa                                                                                                              >  1:1261257/1‑81 (MQ=255)
                                                                                                        gCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCaa                                                                                                              <  2:1261257/81‑1 (MQ=255)
                                                                                                            gACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGga                                     <  2:957605/150‑1 (MQ=255)
                                                                                                              ccTGGACGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGgacg                                   >  2:1536542/1‑150 (MQ=255)
                                                                                                                   ccGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGcc                              <  2:1201803/150‑1 (MQ=255)
                                                                                                                   ccGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTCCACGCTGGCTCACCTGCCCAAGGACCGTATTGa                                                                                                  <  1:8914/81‑1 (MQ=255)
                                                                                                                   ccGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTCCACGCTGGCTCACCTGCCCAAGGACCGTATTGa                                                                                                  >  2:8914/1‑81 (MQ=255)
                                                                                                                    cGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATg                                                                                             >  1:1732744/1‑86 (MQ=255)
                                                                                                                    cGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATg                                                                                             <  2:1732744/86‑1 (MQ=255)
                                                                                                                      aCTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAg                           <  1:1539503/150‑1 (MQ=255)
                                                                                                                              aGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCga                                                                            <  2:249532/93‑1 (MQ=255)
                                                                                                                              aGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCga                                                                            >  1:249532/1‑93 (MQ=255)
                                                                                                                               gCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg                  <  2:1649096/150‑1 (MQ=255)
                                                                                                                                cTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGCCTCACCTGCCCAAGGACTGTATTAATCATCTCCTTGTCCCCCCGCGAGGTCGTCCACGAACTCAACCGCCACATCATCTGCCAGGGCGACCCCAAGCGCTCCGTgg                 >  2:647766/1‑150 (MQ=255)
                                                                                                                                cTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                 <  2:1619634/150‑1 (MQ=255)
                                                                                                                                     cTGGCGCTCCTTCCCAGTGCAGGGGCTTTTCCACGCTCACTTACCTGCCCAAGGACCGCATTGATCATGCCCATGCCCCCCCGCGAGATCGTCCACCAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCAtt            <  1:634315/150‑1 (MQ=255)
                                                                                                                                        gcgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgc                          >  1:1821791/1‑133 (MQ=255)
                                                                                                                                        gcgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgc                          <  2:1821791/133‑1 (MQ=255)
                                                                                                                                           ccccTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                 <  2:723208/139‑1 (MQ=255)
                                                                                                                                           ccccTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                 >  1:723208/1‑139 (MQ=255)
                                                                                                                                            cccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGTCGCCAAGCGCGCCGTGGCCATTGCCCTgc     <  1:453656/150‑1 (MQ=255)
                                                                                                                                              cTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgcgc   >  1:1363237/1‑150 (MQ=255)
                                                                                                                                               ttCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCa  >  1:1062639/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCA  >  NC_002947/5697224‑5697516

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