Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,338 (A)6→5 100% coding (37/963 nt) ldtC ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,3330A.97.5% 352.3 / ‑0.5 80coding (42/963 nt)ldtCL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (43/35);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

ATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAAT  >  NC_000913/1170234‑1170426
                                                                                                   |                                                                                             
aTGAAATTTGTTTTCGCCAACCAGACGACTGCGCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTa                                                                                               <  2:2537214/100‑1 (MQ=255)
    aaTTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaaa                                                                                           <  1:724713/100‑1 (MQ=255)
    aaTTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaaa                                                                                           >  1:146274/1‑100 (MQ=255)
    aaTTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaaa                                                                                           <  2:800521/100‑1 (MQ=255)
    aaTTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaaa                                                                                           <  1:706152/100‑1 (MQ=255)
    aaTTTGTTTGCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaa                                                                                            <  2:586726/99‑1 (MQ=255)
       ttGATTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGGAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaaacg                                                                                         >  1:583927/1‑97 (MQ=255)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTa                                                                                     <  1:596372/100‑1 (MQ=255)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTa                                                                                     >  1:497900/1‑100 (MQ=255)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTa                                                                                     <  1:1422543/100‑1 (MQ=255)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTa                                                                                     <  1:785960/100‑1 (MQ=255)
         gTTTTCGCCAACCACACGACTGCCCGCTGGCGGCAGCGGCCAGGCGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTa                                                                                     >  2:236516/1‑100 (MQ=255)
           tttCGCCAACCAGACGACTGCCCGCTGGCGGCACCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGc                                                                                   >  1:508903/1‑100 (MQ=255)
            ttCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGcc                                                                                  <  1:1612879/100‑1 (MQ=255)
             tCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCa                                                                                 >  1:1170211/1‑100 (MQ=255)
              cGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAg                                                                                >  2:2432653/1‑100 (MQ=255)
               gCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGc                                                                               >  1:1119401/1‑100 (MQ=255)
                ccAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCg                                                                              <  1:2064913/100‑1 (MQ=255)
                ccAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCg                                                                              >  1:587619/1‑100 (MQ=255)
                 cAACCAGACGACTACCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCga                                                                             >  1:1743665/1‑100 (MQ=255)
                    ccAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGaa                                                                          >  1:155530/1‑100 (MQ=255)
                     cAGGCGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGaaa                                                                         <  2:1757174/100‑1 (MQ=255)
                       gacgacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCAAGCAAGaaaaa                                                                       <  2:1119393/100‑1 (MQ=255)
                        acgacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAAc                                                                      <  1:2432667/100‑1 (MQ=255)
                         cggcTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTGGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAATCGTTAGCCAGCGAGAAAAAcg                                                                     <  2:414870/97‑1 (MQ=255)
                          gacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAAcgc                                                                    >  2:658389/1‑100 (MQ=255)
                           ccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAAcgc                                                                    >  2:1193222/2‑99 (MQ=38)
                            cTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGt                                                                  >  2:2163459/1‑100 (MQ=255)
                              gCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGGGAACGGAAAAACGTGAGCCAGCGAGAAAAACGCGtt                                                                 >  2:1734139/1‑99 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               >  2:2191047/1‑100 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               <  1:2234243/100‑1 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               <  1:1185961/100‑1 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               <  2:638554/100‑1 (MQ=255)
                                  gCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGat                                                            >  1:121783/1‑100 (MQ=255)
                                   cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAACCGTTAGCCAGCGAGAAATACGCGTTTTGagc                                                           <  1:1727412/100‑3 (MQ=255)
                                   cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatc                                                           <  1:1916134/100‑1 (MQ=255)
                                    tGGCGGCAGCGGCCAGGTGTTGGCTTTTGCTGGTAGCGTCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatca                                                          <  2:2408611/100‑1 (MQ=255)
                                    tGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCAAGAAAAACGCGTTTTGatc                                                           <  2:809695/99‑1 (MQ=255)
                                     ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcat                                                         >  2:569474/1‑100 (MQ=255)
                                     ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcat                                                         >  2:1914550/1‑100 (MQ=255)
                                      gcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTCGatc                                                           <  2:201015/97‑1 (MQ=255)
                                       cggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGAtcatca                                                       >  2:5003/1‑100 (MQ=255)
                                       cggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGAtcatca                                                       >  1:2485668/1‑100 (MQ=255)
                                        ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCAc                                                      <  1:997801/100‑1 (MQ=255)
                                        ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCAc                                                      <  2:2405705/100‑1 (MQ=255)
                                          cAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGAAAGAAAGAAAAACGCGTTTTGATCATCACCa                                                    <  2:512635/100‑1 (MQ=255)
                                            gCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaa                                                  <  2:928529/100‑1 (MQ=255)
                                            gCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAAAAAGAGAGAAAAACGCGTTTTGATCATCACCaaa                                                  <  2:1347055/100‑1 (MQ=255)
                                             cGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaaa                                                 >  2:2142481/1‑100 (MQ=255)
                                                 cAGGTGTTGGCTTTTGCCGGGAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAAtgc                                             >  2:1560453/1‑98 (MQ=255)
                                                   ggTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCat                                           <  1:872254/100‑1 (MQ=255)
                                                   ggTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCat                                           >  1:2066824/1‑100 (MQ=255)
                                                   ggTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTCTTGATCATCACCAAAAATCCat                                           >  2:1961547/1‑100 (MQ=255)
                                                      gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataa                                         >  1:2340777/1‑100 (MQ=255)
                                                      gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                        >  1:1519497/1‑100 (MQ=255)
                                                      gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                        >  1:639537/1‑100 (MQ=255)
                                                       ttGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataata                                       >  1:793592/1‑100 (MQ=255)
                                                        tGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCAtaataa                                      <  1:2101397/100‑1 (MQ=255)
                                                        tGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCAtaataa                                      >  2:1767598/1‑100 (MQ=255)
                                                         ggCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaa                                     >  1:1713036/1‑100 (MQ=255)
                                                          gCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaat                                    <  2:386104/100‑1 (MQ=255)
                                                          gCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaat                                    <  2:1947876/100‑1 (MQ=255)
                                                           cTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaata                                   >  1:1708242/1‑100 (MQ=255)
                                                            ttttGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAATCGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATACATaataaata                                   >  2:1378438/1‑99 (MQ=255)
                                                               tGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAg                               >  1:1452129/1‑100 (MQ=255)
                                                                gCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAgg                              <  2:242834/100‑1 (MQ=255)
                                                                  cGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGtt                            >  1:1983218/1‑100 (MQ=255)
                                                                   ggTAGCGCCAGCGCCACGTCAGCGGCGAACGTTAAAATGTTAGCCAGCGGGACAGACGGGGATTGATCATCAGCAAACATCCATAATCAATATAGGGTTa                           >  2:542804/1‑100 (MQ=255)
                                                                   ggTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTTTTTTTTTTTTTAAATCCATAATAAATATAAGGTTa                           >  2:162345/1‑100 (MQ=255)
                                                                     tAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTAtt                         >  2:1302320/1‑100 (MQ=255)
                                                                       gcgccagcgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                       >  1:1500931/1‑100 (MQ=255)
                                                                       gcgccagcgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAACATAAGGTTAtttt                       <  1:1736446/100‑3 (MQ=255)
                                                                       gcgccagcgccaCGGAAACGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAAGCCATAATAAATATAAGGTTATTGt                       <  2:2341570/100‑1 (MQ=255)
                                                                         gccagcgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaa                     >  1:1544924/1‑100 (MQ=255)
                                                                           cagcgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaata                   <  1:152751/100‑1 (MQ=255)
                                                                           cagcgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaata                   >  1:1928184/1‑100 (MQ=255)
                                                                            agcgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                  >  1:396883/1‑100 (MQ=255)
                                                                              cgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAAc                >  2:1046487/1‑100 (MQ=255)
                                                                              cgccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAAc                >  1:363426/1‑100 (MQ=255)
                                                                                ccaCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAg              <  1:2149613/100‑1 (MQ=255)
                                                                                  aCGGCAGCTGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGa             <  2:797286/99‑1 (MQ=255)
                                                                                  aCGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGAt            >  2:1507865/1‑100 (MQ=255)
                                                                                   cGGCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATa           <  2:2496735/100‑1 (MQ=255)
                                                                                     gCAGCGGCGAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAg         >  1:1817550/1‑100 (MQ=255)
                                                                                       aGCGGCGAACGTAAAAACGTTAGCAAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGcc       <  2:971877/100‑1 (MQ=255)
                                                                                         cggcgAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATTATAAATATAAGGTTATTGTAATAAACAGATAAGCCTg     <  2:1372733/100‑1 (MQ=255)
                                                                                            cgAACGTAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAAt  >  2:2352498/1‑100 (MQ=255)
                                                                                                   |                                                                                             
ATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAAT  >  NC_000913/1170234‑1170426

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: