Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,645,886 (C)5→4 coding (1440/2313 nt) pbpC ← penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,645,8820C.97.1% 451.0 / 5.7 103coding (1444/2313 nt)pbpCpenicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (49/51);  total (50/53)
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.54e-01

Only 100 of 111 total aligned reads displayed.

TGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGTATACGC  >  NC_000913/2645786‑2645975
                                                                                                |                                                                                             
tGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                                                                                            >  1:225443/1‑100 (MQ=255)
tGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                                                                                            >  2:682011/1‑100 (MQ=255)
tGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                                                                                            <  1:2526539/100‑1 (MQ=255)
 gCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccg                                                                                           >  1:1354868/1‑99 (MQ=255)
  cccGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccga                                                                                          >  1:653178/1‑98 (MQ=255)
  cccGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccga                                                                                          >  1:2391924/1‑98 (MQ=255)
   ccGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccgac                                                                                         >  1:1122133/1‑97 (MQ=255)
    cGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccgacg                                                                                        <  1:583788/100‑5 (MQ=255)
         tCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACAtt                                                                                  <  2:574775/100‑1 (MQ=255)
          ccATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGTCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAtt                                                                                  >  2:2011999/1‑100 (MQ=255)
          ccATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACAtt                                                                                  >  1:393725/1‑99 (MQ=255)
          ccATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTa                                                                                 >  1:1309927/1‑100 (MQ=255)
          ccATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTa                                                                                 >  1:2818035/1‑100 (MQ=255)
          ccATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTa                                                                                 <  1:2134186/100‑1 (MQ=255)
             tGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAg                                                                              >  1:1375121/1‑100 (MQ=255)
              gCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCGATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAgg                                                                             >  2:2093109/1‑100 (MQ=255)
              gCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAgg                                                                             <  2:1994649/100‑1 (MQ=255)
                 aGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCg                                                                          <  1:195499/100‑1 (MQ=255)
                    ggggCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTc                                                                       <  1:721684/100‑1 (MQ=255)
                    ggggCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTc                                                                       <  1:2133318/100‑1 (MQ=255)
                      ggCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCaa                                                                     <  2:1445606/100‑1 (MQ=255)
                       gCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAag                                                                    >  1:1647465/1‑100 (MQ=255)
                           cgcgcgGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAagcagc                                                                >  2:37114/1‑100 (MQ=255)
                            gcgtggCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCg                                                               <  1:434163/96‑1 (MQ=255)
                            gcgcgGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCg                                                               <  1:887876/100‑1 (MQ=255)
                              gcgGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGa                                                             >  1:2233818/1‑100 (MQ=255)
                               cgGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGAt                                                            >  1:908028/1‑100 (MQ=255)
                                ggCAAGGCACTATCCGGCAAGGGTTCCGCTTCAGCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATc                                                           <  2:2622792/100‑1 (MQ=255)
                                 gCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCg                                                          <  1:49267/100‑1 (MQ=255)
                                  cAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGt                                                         <  2:125279/100‑1 (MQ=255)
                                   aaGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTc                                                        <  2:225438/100‑1 (MQ=255)
                                     ggCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAg                                                      >  2:3225182/1‑100 (MQ=255)
                                     ggCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAg                                                      >  1:2766890/1‑100 (MQ=255)
                                       cACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGc                                                    <  2:1313551/100‑1 (MQ=255)
                                       cACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGc                                                    >  1:230996/1‑100 (MQ=255)
                                        aCTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCt                                                   <  1:1711629/100‑1 (MQ=255)
                                         cTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTg                                                  <  2:1029781/100‑1 (MQ=255)
                                          tATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGt                                                 >  1:3044291/1‑100 (MQ=255)
                                           aTCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTACAGGCCGTTCAAGCAGCGGATCCTCAGGCTGt                                                 >  2:527956/1‑99 (MQ=255)
                                           aTCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTa                                                >  1:706626/1‑100 (MQ=255)
                                           aTCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTa                                                <  1:140880/100‑1 (MQ=255)
                                            tCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTaa                                               >  1:187232/1‑100 (MQ=255)
                                             ccGGCAAGGTTTGGGCTTAATCAGCCATATTCCGGCGAAGGATCCCCGCCG‑CCCCGACGACATTAAAGGCCGTTAAGGCAGCGGATCGTCAGGCTGTa                                                <  2:1152387/98‑1 (MQ=255)
                                              cGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAgc                                             >  1:2990504/1‑100 (MQ=255)
                                              cGGCAAGGGTTGCGCTCCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGCCAGGCTGTAAg                                              >  2:1328778/1‑99 (MQ=255)
                                               ggCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAgcg                                            >  2:1550045/1‑100 (MQ=255)
                                               ggCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAgcg                                            >  2:237134/1‑100 (MQ=255)
                                                 cAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCa                                          >  1:1056951/1‑100 (MQ=255)
                                                 cAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCa                                          <  2:2166258/100‑1 (MQ=255)
                                                   aGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAAt                                        >  2:682652/1‑100 (MQ=255)
                                                    gggTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAAtt                                       >  1:1812416/1‑100 (MQ=255)
                                                     ggTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTagatcggaaga                                      >  1:1814841/1‑90 (MQ=255)
                                                     ggTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAAtt                                       >  2:2723327/1‑99 (MQ=255)
                                                      gTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAAtt                                       >  2:1980665/1‑98 (MQ=255)
                                                        tGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcc                                   <  1:2990898/100‑1 (MQ=255)
                                                         acgcTTCATCAGCCATAATCCGGCGATTGATCCACGCTG‑CCCCGTCGACATTAAAGGCCGTTTAAGCAGCGGATTGACAGGATGTAAGCGCAATTtgcct                                  >  2:2861741/2‑100 (MQ=255)
                                                          cgcTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcctg                                 <  2:2877456/100‑1 (MQ=255)
                                                          cgcTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcctg                                 <  2:2414331/100‑1 (MQ=255)
                                                          cgcTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcctg                                 <  1:2761765/100‑1 (MQ=255)
                                                           gcTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcctgc                                >  1:1469740/1‑100 (MQ=255)
                                                           gcTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcctgc                                >  2:1680625/1‑100 (MQ=255)
                                                           gcTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcct                                  >  2:1093621/1‑98 (MQ=255)
                                                           gcTTCATAAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcctgc                                >  1:1713117/1‑100 (MQ=255)
                                                            cTTCATCAGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTtgcctgcc                               >  2:2607105/1‑100 (MQ=255)
                                                             ttCATCAGCCATAATCGGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTgcctgcct                              >  2:404078/1‑100 (MQ=255)
                                                               catcaGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTg                            <  2:474659/100‑1 (MQ=255)
                                                                atcaGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGc                           <  1:2395162/100‑1 (MQ=255)
                                                                 tcaGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGcc                          >  2:2914251/1‑100 (MQ=255)
                                                                   aGCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGt                        >  1:2437108/1‑100 (MQ=255)
                                                                    gCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTg                       <  2:613914/100‑1 (MQ=255)
                                                                    gCCATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTg                       >  1:1172376/1‑100 (MQ=255)
                                                                     ccATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTgg                      >  2:2208511/1‑100 (MQ=255)
                                                                      cATAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGc                     <  2:1309913/100‑1 (MQ=255)
                                                                        tAATCCGGCGAATGATCCACGCCG‑CGCCGCCGACATTCACGGCCGTTCAAGCAGCGGCTCGTCAGGCTGTAAGCGCAATTTGCCTGCCTCGCCGTggt                     <  2:1945491/98‑2 (MQ=255)
                                                                        tAATCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGa                   <  2:197491/100‑1 (MQ=255)
                                                                         aaTCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAg                  <  1:3225219/100‑1 (MQ=255)
                                                                           tCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCa                <  1:819845/100‑1 (MQ=255)
                                                                           tCCGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCa                >  1:391635/1‑100 (MQ=255)
                                                                             cGGCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCaaa              <  1:2093135/100‑1 (MQ=255)
                                                                              ggCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGTAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAAc             >  1:1111419/1‑100 (MQ=255)
                                                                              ggCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAAc             >  1:1063339/1‑100 (MQ=255)
                                                                              ggCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAAc             >  2:1276885/1‑100 (MQ=255)
                                                                               gCGAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAAcg            >  1:7247/1‑100 (MQ=255)
                                                                                 gAATGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAAcgcg          <  2:378827/100‑1 (MQ=255)
                                                                                  aaTGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAAcgctg         <  1:3016066/100‑3 (MQ=255)
                                                                                  aaTGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAAcgc           >  2:284308/1‑98 (MQ=255)
                                                                                  aaTGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCgg         <  2:2818000/100‑1 (MQ=255)
                                                                                  aaTGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCgg         <  1:181826/100‑1 (MQ=255)
                                                                                  aaTGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCgg         <  1:21468/100‑1 (MQ=255)
                                                                                   aTGATCCACGCCGCCCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCgg         <  1:2012023/100‑1 (MQ=255)
                                                                                   aTGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGt        <  1:682020/100‑1 (MQ=255)
                                                                                    tGATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGta       <  1:1217140/100‑1 (MQ=255)
                                                                                    tGATCCACGCCG‑CCCCGACGACATTAAAGGCCGGTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGta       <  1:1548815/100‑1 (MQ=255)
                                                                                     gATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGtat      >  1:3185231/1‑100 (MQ=255)
                                                                                     gATCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGtat      <  2:2858364/100‑1 (MQ=255)
                                                                                      aTCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGtata     >  1:2726485/1‑100 (MQ=255)
                                                                                      aTCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGtata     <  2:2948117/100‑1 (MQ=255)
                                                                                      aTCCACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGtata     <  1:682661/100‑1 (MQ=255)
                                                                                       tCCACGCCG‑CCCCGACGACATTAAAGGCCGTCCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGTAAAc    >  1:25425/1‑100 (MQ=255)
                                                                                         cACGCCG‑CCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGTATACGc  >  2:2161187/1‑100 (MQ=255)
                                                                                                |                                                                                             
TGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTTCAAGCAGCGGATCGTCAGGCTGTAAGCGCAATTTGCCTGCCTTGCCGTGGCGAGCAAACGCGGTATACGC  >  NC_000913/2645786‑2645975

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: ‑