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)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTa                                                                                     <  1:596372/100‑1 (MQ=255)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTa                                                                                     >  1:497900/1‑100 (MQ=255)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTa                                                                                     <  1:1422543/100‑1 (MQ=255)
         gTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTa                                                                                     <  1:785960/100‑1 (MQ=255)
         gTTTTCGCCAACCACACGACTGCCCGCTGGCGGCAGCGGCCAGGCGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTa                                                                                     >  2:236516/1‑100 (MQ=255)
           tttCGCCAACCAGACGACTGCCCGCTGGCGGCACCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGc                                                                                   >  1:508903/1‑100 (MQ=255)
            ttCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGcc                                                                                  <  1:1612879/100‑1 (MQ=255)
             tCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCa                                                                                 >  1:1170211/1‑100 (MQ=255)
              cGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAg                                                                                >  2:2432653/1‑100 (MQ=255)
               gCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGc                                                                               >  1:1119401/1‑100 (MQ=255)
                ccAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCg                                                                              <  1:2064913/100‑1 (MQ=255)
                ccAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCg                                                                              >  1:587619/1‑100 (MQ=255)
                 cAACCAGACGACTACCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCga                                                                             >  1:1743665/1‑100 (MQ=255)
                    ccAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGaa                                                                          >  1:155530/1‑100 (MQ=255)
                     cAGGCGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGaaa                                                                         <  2:1757174/100‑1 (MQ=255)
                       gacgacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCAAGCAAGaaaaa                                                                       <  2:1119393/100‑1 (MQ=255)
                        acgacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAAc                                                                      <  1:2432667/100‑1 (MQ=255)
                         cggcTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTGGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAATCGTTAGCCAGCGAGAAAAAcg                                                                     <  2:414870/97‑1 (MQ=255)
                          gacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAAcgc                                                                    >  2:658389/1‑100 (MQ=255)
                           ccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAAcgc                                                                    >  2:1193222/2‑99 (MQ=38)
                            cTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGt                                                                  >  2:2163459/1‑100 (MQ=255)
                              gCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGGGAACG‑GAAAAACGTGAGCCAGCGAGAAAAACGCGtt                                                                 >  2:1734139/1‑99 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               >  2:2191047/1‑100 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               <  1:2234243/100‑1 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               <  1:1185961/100‑1 (MQ=255)
                               cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGtttt                                                               <  2:638554/100‑1 (MQ=255)
                                  gCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGat                                                            >  1:121783/1‑100 (MQ=255)
                                   cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAACCGTTAGCCAGCGAGAAATACGCGTTTTGagc                                                           <  1:1727412/100‑3 (MQ=255)
                                   cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatc                                                           <  1:1916134/100‑1 (MQ=255)
                                    tGGCGGCAGCGGCCAGGTGTTGGCTTTTGCTGGTAGCGTCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatca                                                          <  2:2408611/100‑1 (MQ=255)
                                    tGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCAAGAAAAACGCGTTTTGatc                                                           <  2:809695/99‑1 (MQ=255)
                                     ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcat                                                         >  2:569474/1‑100 (MQ=255)
                                     ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcat                                                         >  2:1914550/1‑100 (MQ=255)
                                      gcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTCGatc                                                           <  2:201015/97‑1 (MQ=255)
                                       cggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGAtcatca                                                       >  2:5003/1‑100 (MQ=255)
                                       cggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGAtcatca                                                       >  1:2485668/1‑100 (MQ=255)
                                        ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCAc                                                      <  1:997801/100‑1 (MQ=255)
                                        ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCAc                                                      <  2:2405705/100‑1 (MQ=255)
                                          cAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGAAAGAAAGAAAAACGCGTTTTGATCATCACCa                                                    <  2:512635/100‑1 (MQ=255)
                                            gCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaa                                                  <  2:928529/100‑1 (MQ=255)
                                            gCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAAAAAGAGAGAAAAACGCGTTTTGATCATCACCaaa                                                  <  2:1347055/100‑1 (MQ=255)
                                             cGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaaa                                                 >  2:2142481/1‑100 (MQ=255)
                                                 cAGGTGTTGGCTTTTGCCGGGAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAAtgc                                             >  2:1560453/1‑98 (MQ=255)
                                                   ggTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCat                                           <  1:872254/100‑1 (MQ=255)
                                                   ggTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCat                                           >  1:2066824/1‑100 (MQ=255)
                                                   ggTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTCTTGATCATCACCAAAAATCCat                                           >  2:1961547/1‑100 (MQ=255)
                                                      gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataa                                         >  1:2340777/1‑100 (MQ=255)
                                                      gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                        >  1:1519497/1‑100 (MQ=255)
                                                      gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                        >  1:639537/1‑100 (MQ=255)
                                                       ttGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataata                                       >  1:793592/1‑100 (MQ=255)
                                                        tGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCAtaataa                                      <  1:2101397/100‑1 (MQ=255)
                                                        tGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCAtaataa                                      >  2:1767598/1‑100 (MQ=255)
                                                         ggCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaa                                     >  1:1713036/1‑100 (MQ=255)
                                                          gCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaat                                    <  2:386104/100‑1 (MQ=255)
                                                          gCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaat                                    <  2:1947876/100‑1 (MQ=255)
                                                           cTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaata                                   >  1:1708242/1‑100 (MQ=255)
                                                            ttttGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAATCGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATACATaataaata                                   >  2:1378438/1‑99 (MQ=255)
                                                               tGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAg                               >  1:1452129/1‑100 (MQ=255)
                                                                gCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAgg                              <  2:242834/100‑1 (MQ=255)
                                                                  cGGTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGtt                            >  1:1983218/1‑100 (MQ=255)
                                                                   ggTAGCGCCAGCGCCACGTCAGCGGCGAACGT‑TAAAATGTTAGCCAGCGGGACAGACGGGGATTGATCATCAGCAAACATCCATAATCAATATAGGGTTa                           >  2:542804/1‑100 (MQ=255)
                                                                   ggTAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTTTTTTTTTTTTTAAATCCATAATAAATATAAGGTTa                           >  2:162345/1‑100 (MQ=255)
                                                                     tAGCGCCAGCGCCACGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTAtt                         >  2:1302320/1‑100 (MQ=255)
                                                                       gcgccagcgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                       >  1:1500931/1‑100 (MQ=255)
                                                                       gcgccagcgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAACATAAGGTTAtttt                       <  1:1736446/100‑3 (MQ=255)
                                                                       gcgccagcgccaCGGAAACGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAAGCCATAATAAATATAAGGTTATTGt                       <  2:2341570/100‑1 (MQ=255)
                                                                         gccagcgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaa                     >  1:1544924/1‑100 (MQ=255)
                                                                           cagcgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaata                   <  1:152751/100‑1 (MQ=255)
                                                                           cagcgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaata                   >  1:1928184/1‑100 (MQ=255)
                                                                            agcgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                  >  1:396883/1‑100 (MQ=255)
                                                                              cgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAAc                >  2:1046487/1‑100 (MQ=255)
                                                                              cgccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAAc                >  1:363426/1‑100 (MQ=255)
                                                                                ccaCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAg              <  1:2149613/100‑1 (MQ=255)
                                                                                  aCGGCAGCTGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGa             <  2:797286/99‑1 (MQ=255)
                                                                                  aCGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGAt            >  2:1507865/1‑100 (MQ=255)
                                                                                   cGGCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATa           <  2:2496735/100‑1 (MQ=255)
                                                                                     gCAGCGGCGAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAg         >  1:1817550/1‑100 (MQ=255)
                                                                                       aGCGGCGAACGT‑AAAAACGTTAGCAAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGcc       <  2:971877/100‑1 (MQ=255)
                                                                                         cggcgAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATTATAAATATAAGGTTATTGTAATAAACAGATAAGCCTg     <  2:1372733/100‑1 (MQ=255)
                                                                                            cgAACGT‑AAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAAt  >  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: ‑