Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,057,759 +C coding (1773/1776 nt) yihK → putative GTP‑binding factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,057,7581.C99.7% 1113.1 / ‑3.8 296coding (1772/1776 nt)yihKputative GTP‑binding factor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (146/149);  total (146/150)
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.93e-01

Only 100 of 322 total aligned reads displayed.

TATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATG  >  NC_000913/4057666‑4057845
                                                                                              |                                                                                       
tATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCa                                                                                          >  2:2617111/1‑93 (MQ=255)
tATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCa                                                                                          >  1:4102785/1‑93 (MQ=255)
  tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                                       >  1:3933020/1‑92 (MQ=255)
  tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                                       >  1:3298034/1‑92 (MQ=255)
   ccGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcct                                                                                      <  2:832929/93‑3 (MQ=255)
   ccGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcct                                                                                      <  1:2211240/93‑3 (MQ=255)
    cGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctg                                                                                     >  2:5882874/1‑90 (MQ=255)
    cGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctg                                                                                     >  1:2880828/1‑90 (MQ=255)
     gCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctga                                                                                    <  2:232452/93‑5 (MQ=255)
       aTGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACg                                                                                   >  2:2997027/1‑93 (MQ=255)
         gACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGa                                                                                 >  1:1293358/1‑93 (MQ=255)
          aCTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaa                                                                                >  2:1800860/1‑93 (MQ=255)
          aCTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaa                                                                                >  2:4199817/1‑93 (MQ=255)
           ctctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                                               >  2:6219199/1‑93 (MQ=255)
            tctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaaa                                                                              <  1:1040307/93‑1 (MQ=255)
               ggAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGa                                                                           <  1:6648542/93‑1 (MQ=255)
               ggAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGa                                                                           <  2:465943/93‑1 (MQ=255)
                gAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAt                                                                          <  1:6178345/93‑1 (MQ=255)
                  aCAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcg                                                                        <  1:3965955/93‑1 (MQ=255)
                   cAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgt                                                                       >  2:6345559/1‑93 (MQ=255)
                   cAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgt                                                                       <  1:213346/93‑1 (MQ=255)
                    aatctctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                                      <  1:224108/90‑1 (MQ=255)
                    aaGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                                      <  2:1521636/93‑1 (MQ=255)
                    aaGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                                      <  2:3550074/93‑1 (MQ=255)
                     aGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtcg                                                                     >  2:1889426/1‑93 (MQ=255)
                     aGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtcg                                                                     >  2:4337438/1‑93 (MQ=255)
                      gCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgc                                                                    >  2:2681337/1‑93 (MQ=255)
                       ctctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                                   <  2:1358294/93‑1 (MQ=255)
                          tGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCcgcg                                                                >  2:2119427/1‑93 (MQ=255)
                          tGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCcgcg                                                                >  2:5391913/1‑93 (MQ=255)
                           tgAGTTCAACGATGATGAAGAACTGGTAGAAGTGACTCCGACCTCTATACGTATTCGTAACCGTCACCTGACGGAAAACGATCGTCGCCgcgc                                                               <  2:3233480/92‑1 (MQ=255)
                           ggAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCgcgc                                                               <  1:3023418/93‑1 (MQ=255)
                           ggAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCgcgc                                                               >  1:1447527/1‑93 (MQ=255)
                            gAGTTCATCGAGGATGACGAGCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGcc                                                              <  2:4585181/93‑1 (MQ=255)
                              gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCaa                                                            >  1:2547563/1‑93 (MQ=255)
                              gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCaa                                                            >  2:1461129/1‑93 (MQ=255)
                               ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCCCCTGACGGAAAACGATCGTCGCCGCGCCAAc                                                           <  1:5866231/93‑1 (MQ=255)
                               ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAc                                                           >  2:3192726/1‑93 (MQ=255)
                               ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAc                                                           <  1:5823371/93‑1 (MQ=255)
                                 cATCGATGATGACGAACTTGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcg                                                         <  1:3263382/93‑1 (MQ=255)
                                    cGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCgcgc                                                      <  2:103620/93‑1 (MQ=255)
                                      atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAc                                                    >  2:1611997/1‑93 (MQ=255)
                                      atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAc                                                    <  2:574085/93‑1 (MQ=255)
                                       tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                                   <  2:4962223/93‑1 (MQ=255)
                                       agatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                                   <  2:5225953/92‑1 (MQ=255)
                                         atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGa                                                 >  1:8086/1‑93 (MQ=255)
                                         atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGa                                                 >  1:1588629/1‑93 (MQ=255)
                                          tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaa                                                >  1:6512239/1‑93 (MQ=255)
                                           gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaaa                                               <  1:2009750/93‑1 (MQ=255)
                                            aCGAACTGGTAGAAGTGCCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAg                                              <  1:5696043/93‑1 (MQ=255)
                                             cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGa                                             <  2:3583222/93‑1 (MQ=255)
                                             cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGa                                             <  2:6449155/93‑1 (MQ=255)
                                              gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGAc                                            >  1:6055033/1‑93 (MQ=255)
                                               aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                                           <  2:2241236/93‑1 (MQ=255)
                                               aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                                           >  1:1100408/1‑93 (MQ=255)
                                                aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                                          <  1:3097839/93‑1 (MQ=255)
                                                 cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                                         <  2:2477051/93‑1 (MQ=255)
                                                      aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAtt                                    <  2:3432404/93‑1 (MQ=255)
                                                       gAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAttt                                   >  2:5478151/1‑93 (MQ=255)
                                                        aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTc                                  <  1:98195/93‑1 (MQ=255)
                                                           tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  2:1273716/93‑1 (MQ=255)
                                                           tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  2:4021717/93‑1 (MQ=255)
                                                           tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  2:1157901/93‑1 (MQ=255)
                                                             aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttt                             <  1:3130340/93‑1 (MQ=255)
                                                              cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt                            <  2:3720655/93‑1 (MQ=255)
                                                              cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt                            >  2:1383887/1‑93 (MQ=255)
                                                              cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt                            <  2:5650517/93‑1 (MQ=255)
                                                               tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttttt                           >  1:4026364/1‑93 (MQ=255)
                                                                ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTa                          >  2:3404617/1‑93 (MQ=255)
                                                                ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTa                          >  2:5135330/1‑93 (MQ=255)
                                                                  gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGt                        >  2:6619606/1‑93 (MQ=255)
                                                                   aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTc                       >  2:4443508/1‑93 (MQ=255)
                                                                   aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACAGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTc                       >  1:3215261/1‑93 (MQ=255)
                                                                    ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa                      >  2:6023992/1‑93 (MQ=255)
                                                                    ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa                      >  2:6153301/1‑93 (MQ=255)
                                                                      tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATa                    >  1:4889157/1‑93 (MQ=255)
                                                                      tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATa                    <  1:6237011/93‑1 (MQ=255)
                                                                       ttATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaa                   <  2:4108684/92‑1 (MQ=255)
                                                                         aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaa                 >  1:3653845/1‑93 (MQ=255)
                                                                         aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaa                 >  2:3839641/1‑93 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                >  1:5147271/1‑93 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                <  1:5404143/93‑1 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGACATaaaaa                >  1:4472142/1‑93 (MQ=255)
                                                                           ccGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaaa               <  1:2726891/93‑1 (MQ=255)
                                                                            cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAc              <  1:870193/93‑1 (MQ=255)
                                                                            cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAc              >  2:1232915/1‑93 (MQ=255)
                                                                             gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAcc             >  2:5861870/1‑93 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAACACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCt            <  2:4868026/93‑1 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCt            >  2:5526960/1‑93 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCt            <  1:5350181/93‑1 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCt            <  1:3543126/93‑1 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCt            <  1:4963696/93‑1 (MQ=255)
                                                                                  cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCa        <  1:2465645/93‑1 (MQ=255)
                                                                                   gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg       >  1:4707642/1‑93 (MQ=255)
                                                                                    tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGc      >  1:5222720/1‑93 (MQ=255)
                                                                                     aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCTTTTTTTAGTCATAAAAAACCTGCCAGCg     >  2:6565131/1‑93 (MQ=255)
                                                                                       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt   >  1:3734187/1‑93 (MQ=255)
                                                                                       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt   <  1:1156663/93‑1 (MQ=255)
                                                                                        cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg  >  1:3854817/1‑93 (MQ=255)
                                                                                        cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg  <  2:685290/93‑1 (MQ=255)
                                                                                              |                                                                                       
TATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATG  >  NC_000913/4057666‑4057845

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