Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,764,854 C→T 10.1% A260V (GCC→GTC)  ctpA → carboxy‑terminal‑processing protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,764,8540CT10.1% 183.4 / 10.6 69A260V (GCC→GTC) ctpAcarboxy‑terminal‑processing protease
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base T (4/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

CCGTATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGG  >  NC_002947/5764705‑5764994
                                                                                                                                                     |                                                                                                                                            
ccGTATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGc                                                                                                                                              >  1:1026251/1‑150 (MQ=255)
 cGTATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGGGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGcc                                                                                                                                             <  2:63567/150‑1 (MQ=255)
 cGTATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGcc                                                                                                                                             >  2:166452/1‑150 (MQ=255)
  gTATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCg                                                                                                                                            >  2:1265323/1‑150 (MQ=255)
    aTCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGAc                                                                                                                                          <  1:1256887/150‑1 (MQ=255)
       aCCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCAc                                                                                                                                       >  1:561272/1‑150 (MQ=255)
        cccAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAACGCCCTGGCCGAGCTGCGCAAGGACAACGGCAAGCAGCTTCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACt                                                                                                                                      <  1:455435/150‑1 (MQ=255)
         ccAGTTCCAGGTCAAGACCGGCGACGAAGTGGGCAAGTTCCTGGCCAAGCTGCGCACGGAGAACGGCAGGAAGCTGCGCGGGCTGGTACTGGACCTGCGCAACAACCCAGGCGGCGTGGTGCGATCGGCAGTGGAAGTGGCCGACCACtt                                                                                                                                     >  2:309878/1‑150 (MQ=255)
              tCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGa                                                                                                                                >  1:10380/1‑150 (MQ=255)
              tCCAGGTCAAGACCGGTGACCAAGTGAGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGTTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGa                                                                                                                                <  1:1022781/150‑1 (MQ=255)
                     cAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAgg                                                                                                                         >  2:391082/1‑150 (MQ=255)
                        gACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGCCGTTCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCt                                                                                                                      >  1:1379704/1‑150 (MQ=255)
                            ggTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACCACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATc                                                                                                                  <  2:1026251/150‑1 (MQ=255)
                              tGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGt                                                                                                                <  1:348995/150‑1 (MQ=255)
                                    aGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTacac                                                                                                          >  1:176722/1‑150 (MQ=255)
                                     gTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                             >  1:239282/1‑146 (MQ=255)
                                     gTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                             <  2:239282/146‑1 (MQ=255)
                                         gCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGCGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAg                                                                                                     >  2:925895/1‑150 (MQ=255)
                                          cAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAgg                                                                                                    <  1:575808/150‑1 (MQ=255)
                                               cccTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGt                                                                                                                        >  1:829276/1‑125 (MQ=255)
                                               cccTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGt                                                                                                                        <  2:829276/125‑1 (MQ=255)
                                                ccTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGATGCTCGGGGTGAAGCTGGACCGGCGCAACAACCCAGGCGGCATGCTCAGACCTTCAGCGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTATACACCTAAGGCCGCAt                                                                                              >  1:1391009/1‑150 (MQ=255)
                                                      cAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCaa                                                                                        >  2:922756/1‑150 (MQ=255)
                                                                 aGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGc                                                                             <  2:561272/150‑1 (MQ=255)
                                                                  ggACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCg                                                                            >  1:1274425/1‑150 (MQ=255)
                                                                         ggCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTc                                                                                                                       >  1:1004037/1‑100 (MQ=255)
                                                                         ggCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTc                                                                                                                       <  2:1004037/100‑1 (MQ=255)
                                                                         ggCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCg                                                                     <  1:166452/150‑1 (MQ=255)
                                                                           cAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTc                                                                                                                                    <  1:1083916/85‑1 (MQ=255)
                                                                           cAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTc                                                                                                                                    >  2:1083916/1‑85 (MQ=255)
                                                                            aagaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                 >  2:276130/1‑85 (MQ=255)
                                                                            aagaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                 <  1:276130/87‑3 (MQ=255)
                                                                            aagaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCACCAACTCCGAGCTGCGTTTCTCGGcc                                                                  <  1:315862/150‑1 (MQ=255)
                                                                            aagaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGCCGTGCTGCACTCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGcc                                                                  >  1:683166/1‑150 (MQ=255)
                                                                                   tGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGGCTGATCGTCTACACGAAAGGCCGCATCGCCAACTCCGAGCTGCGGTTCTCGGCCGAcccag                                                           >  2:653739/1‑148 (MQ=255)
                                                                                    gcgcGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGGCGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGc                                                          <  2:550068/150‑1 (MQ=255)
                                                                                        ggggTGGTGCTGGACCTGCGCAACACCCCAGGCGGCGTTCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGAt                                                      <  1:653739/150‑1 (MQ=255)
                                                                                          ggtggtGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGc                                                    <  2:1094379/150‑1 (MQ=255)
                                                                                            tggtgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCa                                                  <  1:391082/150‑1 (MQ=255)
                                                                                               tgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                 <  1:261263/68‑3 (MQ=38)
                                                                                               tgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                 >  2:261263/1‑66 (MQ=38)
                                                                                                 cTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGcgcg                                             >  2:232361/1‑150 (MQ=255)
                                                                                                  tGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCgg                                            <  2:176722/150‑1 (MQ=255)
                                                                                                  acaaCCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCgg                                            >  2:1065534/4‑150 (MQ=255)
                                                                                                     aCCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCg                                         >  1:1239966/1‑150 (MQ=255)
                                                                                                       cTGCGCAACAACCCAGGTGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTg                                       >  1:601864/1‑150 (MQ=255)
                                                                                                        tGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGACGTGGCCCACCACTTCCTGACCCAAGGTCAGACCGTCGACACCAAAGGCAGCCAGGGCAATTCCATGCTGGGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGc                                      >  1:1386758/1‑150 (MQ=255)
                                                                                                            caacaaCCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCt                                  <  2:1233130/150‑1 (MQ=255)
                                                                                                            caacaaCCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCt                                  >  1:1233130/1‑150 (MQ=255)
                                                                                                               caacCCAGGCGGCGTGCTGCACGCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggt                               >  2:172167/1‑150 (MQ=255)
                                                                                                                aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGTCCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                              >  2:1267416/1‑150 (MQ=255)
                                                                                                                 acCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAg                                                                                                                          >  2:997273/1‑57 (MQ=255)
                                                                                                                 acCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAg                                                                                                                          <  1:997273/57‑1 (MQ=255)
                                                                                                                   ccAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTa                           <  2:1239966/150‑1 (MQ=255)
                                                                                                                      ggcggcGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                           >  2:26024/1‑67 (MQ=255)
                                                                                                                      ggcggcGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                           <  1:26024/67‑1 (MQ=255)
                                                                                                                          gcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAg                                                                                 <  2:981255/89‑1 (MQ=255)
                                                                                                                          gcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAg                                                                                 >  1:981255/1‑89 (MQ=255)
                                                                                                                                tgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggt                               >  2:1175548/1‑133 (MQ=255)
                                                                                                                                tgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggt                               <  1:1175548/133‑1 (MQ=255)
                                                                                                                                    aTCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGc                                                          >  1:666438/1‑102 (MQ=255)
                                                                                                                                    aTCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGc                                                          <  2:666438/102‑1 (MQ=255)
                                                                                                                                       ggCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGt                                        >  2:255648/1‑117 (MQ=255)
                                                                                                                                       ggCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGt                                        <  1:255648/117‑1 (MQ=255)
                                                                                                                                         cAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcat     <  2:601864/150‑1 (MQ=255)
                                                                                                                                            tGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtgg                                >  2:1077704/1‑120 (MQ=255)
                                                                                                                                            tGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtgg                                <  1:1077704/120‑1 (MQ=255)
                                                                                                                                            tGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcgg  >  1:554537/1‑150 (MQ=255)
                                                                                                                                            tGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcgg  >  2:1073919/1‑150 (MQ=255)
                                                                                                                                            aggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGAc                                                               >  2:1008024/2‑89 (MQ=255)
                                                                                                                                            aggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGAc                                                               <  1:1008024/88‑1 (MQ=255)
                                                                                                                                            aggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcgg  >  1:514851/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
CCGTATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGG  >  NC_002947/5764705‑5764994

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: