Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,644,603 (G)6→5 coding (743/2313 nt) pbpC ← putative peptidoglycan enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,644,5980G.99.6% 1062.0 / ‑3.7 236coding (748/2313 nt)pbpCputative peptidoglycan enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (110/125);  total (110/126)
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.07e-01

Only 100 of 269 total aligned reads displayed.

AAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTA  >  NC_000913/2644506‑2644686
                                                                                            |                                                                                        
aaGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCg                                                                                          >  2:2866962/1‑93 (MQ=255)
aaGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCg                                                                                          <  2:1032983/93‑1 (MQ=255)
 aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  2:8256465/93‑1 (MQ=255)
 aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  2:7784055/93‑1 (MQ=255)
  gACCGCCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                        <  1:1499108/93‑1 (MQ=255)
    ccGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggggg                                                                                      <  1:3279425/93‑1 (MQ=255)
     cGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggc                                                                                     >  1:4458949/1‑92 (MQ=255)
      ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    >  1:6095332/1‑91 (MQ=255)
      ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    <  1:7815265/93‑3 (MQ=255)
       gCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcca                                                                                   >  1:6933374/1‑90 (MQ=255)
        cATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggccag                                                                                  >  2:2164000/1‑89 (MQ=255)
         aTCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGgccag                                                                                 >  1:4731224/1‑93 (MQ=255)
         aTCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagc                                                                                <  1:130497/93‑1 (MQ=255)
         aTCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagc                                                                                >  1:2009732/1‑93 (MQ=255)
          tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagcc                                                                               <  2:9281683/93‑1 (MQ=255)
          tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGCTGCGGCATTTGTCGGGGgccagcc                                                                               >  1:2416919/1‑93 (MQ=255)
           ccAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagcca                                                                              >  1:7414196/1‑93 (MQ=255)
           ccAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGgccagcca                                                                              <  2:4503436/93‑1 (MQ=255)
            cAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGccagccag                                                                             >  1:7605684/1‑93 (MQ=255)
             aaCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGTTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGa                                                                            <  1:3754220/93‑1 (MQ=255)
              aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCGGCCAGAt                                                                           >  1:2730123/1‑93 (MQ=255)
              aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGAt                                                                           <  1:70249/93‑1 (MQ=255)
              aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGAt                                                                           >  1:937434/1‑93 (MQ=255)
                 tagtagTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATggg                                                                        <  2:6656215/93‑1 (MQ=255)
                 tagtagTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATggg                                                                        >  2:4900630/1‑93 (MQ=255)
                 tagtagTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATggg                                                                        <  1:9524150/93‑1 (MQ=255)
                 tagtagTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATggg                                                                        >  1:9826639/1‑93 (MQ=255)
                  agtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGt                                                                       <  2:544631/93‑1 (MQ=255)
                   gtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGtt                                                                      >  1:4313723/1‑93 (MQ=255)
                    tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttc                                                                     >  2:1855675/1‑93 (MQ=255)
                     agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttct                                                                    <  2:7876261/93‑1 (MQ=255)
                      gtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttctt                                                                   >  2:5588340/1‑93 (MQ=255)
                       tGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttcttc                                                                  >  2:2527701/1‑93 (MQ=255)
                       tGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGttcttc                                                                  >  2:7671355/1‑93 (MQ=255)
                        gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTcc                                                                 >  2:2311706/1‑93 (MQ=255)
                         aTTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTccc                                                                >  2:4148328/1‑93 (MQ=255)
                            ttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTg                                                             >  1:4989848/1‑93 (MQ=255)
                             tGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGa                                                            <  1:2233463/93‑1 (MQ=255)
                             tGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGa                                                            <  2:9383381/93‑1 (MQ=255)
                              gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGAc                                                           <  1:6057782/93‑1 (MQ=255)
                              gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGAc                                                           <  1:5410041/93‑1 (MQ=255)
                                cGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActc                                                         <  2:9176178/93‑1 (MQ=255)
                                cGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActc                                                         >  2:5406825/1‑93 (MQ=255)
                                 gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCTAGCCAGATGGGTTCTTCCCTTGActct                                                        <  1:1323237/93‑1 (MQ=255)
                                 gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                        >  1:2509832/1‑93 (MQ=255)
                                 gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                        >  2:6986533/1‑93 (MQ=255)
                                  cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                       <  1:6335695/93‑1 (MQ=255)
                                   ttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCttt                                                      <  1:7292328/93‑1 (MQ=255)
                                    tttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtttt                                                     >  2:8076473/1‑93 (MQ=255)
                                    tttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtttt                                                     >  1:9113964/1‑93 (MQ=255)
                                     ttgcttttACCGAGCATCATGCGCGAAAGCAGCGGTGCCAGTTGCGGCATTTGTCGGGGTCCAGCCAGATGGGTTCTTCCCTTGACCCTTTTa                                                    <  2:1485421/93‑1 (MQ=255)
                                     ttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTa                                                    <  2:6043936/93‑1 (MQ=255)
                                     ttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTa                                                    >  1:6731902/1‑93 (MQ=255)
                                      tgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc                                                   <  2:3963369/93‑1 (MQ=255)
                                      tgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTCAc                                                   >  1:5090150/1‑93 (MQ=255)
                                        cttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCt                                                 >  1:9055315/1‑93 (MQ=255)
                                            tACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTc                                             >  2:1593723/1‑93 (MQ=255)
                                                gAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgg                                         >  2:6823615/1‑93 (MQ=255)
                                                gAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACgg                                         <  2:5932165/93‑1 (MQ=255)
                                                 aGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACggg                                        >  1:5695897/1‑93 (MQ=255)
                                                  gCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGa                                       >  2:612333/1‑93 (MQ=255)
                                                   catcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGa                                       <  2:9723640/93‑1 (MQ=255)
                                                    atcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGAcc                                     <  1:9681067/93‑1 (MQ=255)
                                                    atcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGGTCTTCCCTTGACTCTTTTACCTGCTCACGGGAcc                                     >  2:9222917/1‑93 (MQ=255)
                                                     tcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTTCTCACGGAACca                                    >  1:1935297/1‑93 (MQ=255)
                                                     tcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACca                                    >  2:119642/1‑93 (MQ=255)
                                                     tcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACca                                    >  2:7143251/1‑93 (MQ=255)
                                                     tcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACca                                    >  2:5904803/1‑93 (MQ=255)
                                                       atGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcaca                                  <  1:5893612/93‑1 (MQ=255)
                                                        tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                 <  1:4015239/93‑1 (MQ=255)
                                                         gcgcGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacaca                                <  2:3360435/93‑1 (MQ=255)
                                                         gcgcGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacaca                                <  2:3492124/93‑1 (MQ=255)
                                                         gcgcGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAATTGGGATCTTCCCTTTACTCTTTTACCTGCTCACGGGACcacaca                                >  2:8653321/1‑93 (MQ=255)
                                                           gcgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACAcc                              <  1:2970433/93‑1 (MQ=255)
                                                           gcgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACAcc                              <  1:171712/93‑1 (MQ=255)
                                                             gAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt                            <  1:4518518/93‑1 (MQ=255)
                                                             gAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt                            >  2:4326564/1‑93 (MQ=255)
                                                               acaCAGCGGTGCGAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGc                          <  1:8826640/91‑1 (MQ=255)
                                                                aaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCa                         <  1:4044773/93‑1 (MQ=255)
                                                                  cAGCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACg                       <  1:3578638/93‑1 (MQ=255)
                                                                    gCGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGc                     <  1:212090/93‑1 (MQ=255)
                                                                     cGGTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGcc                    <  1:5794392/93‑1 (MQ=255)
                                                                      ggTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGATTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCa                   <  1:5334322/93‑1 (MQ=255)
                                                                       gTGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCAt                  >  1:8283315/1‑93 (MQ=255)
                                                                        tGCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATc                 >  1:6254601/1‑93 (MQ=255)
                                                                         gCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATcc                <  1:8530202/93‑1 (MQ=255)
                                                                         gCCAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATcc                <  2:4610149/93‑1 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt              >  2:1328746/1‑93 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt              <  1:5632704/93‑1 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt              <  2:7306058/93‑1 (MQ=255)
                                                                           cAGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt              <  2:6607333/93‑1 (MQ=255)
                                                                            aGTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGtt             <  2:1555275/93‑1 (MQ=255)
                                                                            aGTTGCGGCATGTGGCGGGGGGAAGCAAGATGGGTTCTTCCCTTGACTCTTTTGCCTGCTCACGGGACCACACACCTTGCGCGGCCATCCGtt             <  1:3470388/93‑1 (MQ=255)
                                                                             gTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCCTCCGTTc            >  1:9001857/1‑93 (MQ=255)
                                                                             gTTGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc            <  1:3836727/93‑1 (MQ=255)
                                                                              ttGCGGCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCg           >  2:6517709/1‑93 (MQ=255)
                                                                                   gCATTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGAACCACACACCTTGCACGGCCATCCGTTCGAGCAc      <  2:6372062/93‑1 (MQ=255)
                                                                                     aTTTGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACtt    >  2:5026801/1‑93 (MQ=255)
                                                                                      tttGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACttt   >  1:10026986/1‑93 (MQ=255)
                                                                                      tttGTCGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACttt   <  1:4390618/93‑1 (MQ=255)
                                                                                            |                                                                                        
AAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTA  >  NC_000913/2644506‑2644686

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