Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 470,566 G→C 5.6% D566H (GAT→CAT)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913470,5660GC 5.6% 254.2 / 6.0 71D566H (GAT→CAT) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base C (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

AGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGG  >  NC_000913/470420‑470711
                                                                                                                                                  |                                                                                                                                                 
aGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGat                                                                                                                                                 <  2:35592/149‑1 (MQ=255)
aGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATCTGCTGGCACAACAAAGCGGCTGGTATCGCGat                                                                                                                                                 <  4:686450/149‑1 (MQ=255)
 gTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGata                                                                                                                                                <  1:704770/149‑1 (MQ=255)
 gTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGata                                                                                                                                                >  2:604079/1‑149 (MQ=255)
 gTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCCata                                                                                                                                                >  3:606332/1‑149 (MQ=255)
 gTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCCata                                                                                                                                                >  3:510357/1‑149 (MQ=255)
      ggCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTAt                                                                                                                                           <  2:85525/149‑1 (MQ=255)
         aGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGc                                                                                                                                        >  4:622637/1‑149 (MQ=255)
          gggAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCt                                                                                                                                       >  2:200296/1‑149 (MQ=255)
             aagaaCGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATc                                                                                                                                    >  4:89663/1‑149 (MQ=255)
                   ggTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCaacaac                                                                                                                              <  3:286303/149‑1 (MQ=255)
                   ggTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCaacaac                                                                                                                              <  1:791005/149‑1 (MQ=255)
                    gTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACt                                                                                                                             >  4:132919/1‑149 (MQ=255)
                     tAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTa                                                                                                                                      >  1:3893/1‑139 (MQ=255)
                     tAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTa                                                                                                                                      <  2:3893/139‑1 (MQ=255)
                     tAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTg                                                                                                                            <  4:478041/149‑1 (MQ=255)
                     tAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTg                                                                                                                            <  2:541420/149‑1 (MQ=255)
                     tAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTg                                                                                                                            <  2:254306/149‑1 (MQ=255)
                       atcatcAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGa                                                                                                                          >  2:542870/1‑149 (MQ=255)
                       atcatcAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGa                                                                                                                          <  1:54972/149‑1 (MQ=255)
                          atcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAtgc                                                                                                                       >  3:402256/1‑149 (MQ=255)
                          atcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAggc                                                                                                                       >  2:11059/1‑149 (MQ=255)
                          atcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAggc                                                                                                                       >  3:666921/1‑149 (MQ=255)
                          atcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAggc                                                                                                                       >  1:743933/1‑149 (MQ=255)
                              gTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGTTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGgcggcg                                                                                                                   >  3:489734/1‑149 (MQ=255)
                                             cTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCg                                                                                                    >  4:680034/1‑149 (MQ=255)
                                                cACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGaa                                                                                                 <  2:392590/149‑1 (MQ=255)
                                                  cTGACGCAAGGCAGTGAAATTATTGTGTAGCAGCACGGTCATATCGCCCATCGTGGCAATCATGATGTGCTATCACAACAAAGCGGCTGTGATCGCGATATGTATCGCTAGCAACAACTGGAGGCGGCGCTCGACGACGTTCCGGataa                                                                                               <  4:408712/149‑1 (MQ=255)
                                                    gACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATc                                                                                             <  3:622637/149‑1 (MQ=255)
                                                     aCGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcg                                                                                            >  2:111947/1‑149 (MQ=255)
                                                     aCGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcg                                                                                            >  3:304219/1‑149 (MQ=255)
                                                         aaGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCgag                                                                                        <  1:542870/149‑1 (MQ=255)
                                                            ccAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGTTATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCgaggag                                                                                     >  3:461/1‑149 (MQ=255)
                                                             cAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGaggagg                                                                                    >  4:64936/1‑149 (MQ=255)
                                                               gTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGcc                                                                                  >  2:719728/1‑149 (MQ=255)
                                                                tGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCg                                                                                 >  4:590506/1‑149 (MQ=255)
                                                                tGAAATTATTGTGATGCAGCACGCACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCg                                                                                 >  4:49383/1‑149 (MQ=255)
                                                                    cttattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaacaa                                                                                                                               <  2:489131/98‑1 (MQ=255)
                                                                    attattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTAt                                                                                                                                     >  4:557589/1‑93 (MQ=255)
                                                                    attattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTAt                                                                                                                                     <  3:557589/93‑1 (MQ=255)
                                                                    attattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaacaa                                                                                                                               >  1:489131/1‑99 (MQ=255)
                                                                       attGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGTTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGc                                                                          <  4:489734/149‑1 (MQ=255)
                                                                           tgATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAg                                                                      >  3:181956/1‑149 (MQ=255)
                                                                             aTGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCCATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGtt                                                                    <  4:510357/149‑1 (MQ=255)
                                                                             aTGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCCATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGtt                                                                    <  4:606332/149‑1 (MQ=255)
                                                                                    aCGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGAt                                                                            >  1:182540/1‑134 (MQ=255)
                                                                                    aCGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGAt                                                                            <  2:182540/134‑1 (MQ=255)
                                                                                          atatCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACtgtg                                                       <  3:89663/149‑1 (MQ=255)
                                                                                             tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGTTATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGcc                                                    <  4:461/149‑1 (MQ=255)
                                                                                             tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGcc                                                    >  4:684964/1‑149 (MQ=255)
                                                                                               gCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGa                                                  <  3:590506/149‑1 (MQ=255)
                                                                                                cccAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTCAGCCAACTGTGGCCGAc                                                 >  1:180678/1‑149 (MQ=255)
                                                                                                 ccAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTcgacga                                                                                                           >  3:537337/1‑90 (MQ=255)
                                                                                                 ccAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGATGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTcgacga                                                                                                           <  4:537337/90‑1 (MQ=255)
                                                                                                      gTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCaa                                           <  1:462857/149‑1 (MQ=255)
                                                                                                      gTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCaa                                           <  2:734793/149‑1 (MQ=255)
                                                                                                         gCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAgcg                                        >  1:30152/1‑149 (MQ=255)
                                                                                                           aaTCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAAt                                                                                              >  1:459518/1‑93 (MQ=255)
                                                                                                           aaTCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAAt                                                                                              <  2:459518/93‑1 (MQ=255)
                                                                                                                tgatgTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTa                                 >  3:2944/1‑149 (MQ=255)
                                                                                                                tgatgTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTa                                 >  4:520199/1‑149 (MQ=255)
                                                                                                                     tgCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTCAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                            <  2:180678/149‑1 (MQ=255)
                                                                                                                        tGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACgg                         >  1:224988/1‑149 (MQ=255)
                                                                                                                          gCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGtt                       >  2:709988/1‑149 (MQ=255)
                                                                                                                            acaacaAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCg                     >  3:711363/1‑149 (MQ=255)
                                                                                                                                caaAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGt                 <  1:11059/149‑1 (MQ=255)
                                                                                                                                caaAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGt                 >  4:327654/1‑149 (MQ=255)
                                                                                                                                 aaaGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTg                <  2:743933/149‑1 (MQ=255)
                                                                                                                                 aaaGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTg                <  4:666921/149‑1 (MQ=255)
                                                                                                                                  aaGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTgg               <  4:402256/149‑1 (MQ=255)
                                                                                                                                         tGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAAc        >  4:284046/1‑149 (MQ=255)
                                                                                                                                               cgcgATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTgg  <  2:6191/149‑1 (MQ=255)
                                                                                                                                               cgcgATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTgg  <  4:181956/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGG  >  NC_000913/470420‑470711

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