Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,764,857 A→G 10.3% D261G (GAC→GGC)  ctpA → carboxy‑terminal‑processing protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,764,8570AG10.3% 182.7 / 10.2 68D261G (GAC→GGC) ctpAcarboxy‑terminal‑processing protease
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base G (4/3);  total (37/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 = 2.59e-01

ATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGG  >  NC_002947/5764709‑5764994
                                                                                                                                                    |                                                                                                                                         
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                                                                                                                                    >  2:1083916/1‑85 (MQ=255)
                                                                       cAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTc                                                                                                                                    <  1:1083916/85‑1 (MQ=255)
                                                                        aagaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                 <  1:276130/87‑3 (MQ=255)
                                                                        aagaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                 >  2:276130/1‑85 (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                                                                                                                                 >  2:261263/1‑66 (MQ=38)
                                                                                           tgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                 <  1:261263/68‑3 (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                                  >  1:1233130/1‑150 (MQ=255)
                                                                                                        caacaaCCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCt                                  <  2:1233130/150‑1 (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                                                                                                           <  1:26024/67‑1 (MQ=255)
                                                                                                                  ggcggcGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                           >  2:26024/1‑67 (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  >  2:1073919/1‑150 (MQ=255)
                                                                                                                                        tGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcgg  >  1:554537/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)
                                                                                                                                                    |                                                                                                                                         
ATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGG  >  NC_002947/5764709‑5764994

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