Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 6.6% intergenic (+22/+36) zupT → / ← yqiD heavy metal divalent cation transporter ZupT/protein YqiD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,3451.T 6.6% 225.7 / 11.0 61intergenic (+22/+36)zupT/yqiDheavy metal divalent cation transporter ZupT/protein YqiD
Reads supporting (aligned to +/- strand):  ref base . (29/28);  new base T (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-02

TGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCA  >  NC_000913/3183203‑3183487
                                                                                                                                               |                                                                                                                                              
tgaGTTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGCTGTCCa                                                                                                                                                 >  2:651436/1‑143 (MQ=255)
tgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                                                                 >  2:934694/1‑143 (MQ=255)
 gaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATACTAACCCCAGCTATGGCTTATTTTGTGGTATTACAGAAAAGGGATTAAGTTTAGTGATGCTACAAACGGCGGGAATTGGTTAAAAACAACATAAGAAGTTGTCAAg                                                                                                                                               >  1:149758/1‑143 (MQ=255)
 gaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAg                                                                                                                                               >  2:125455/1‑143 (MQ=255)
   aTTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcg                                                                                                                                             >  2:907324/1‑143 (MQ=255)
    ttAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                                                                            >  2:286794/1‑143 (MQ=255)
    ttAATGCCGCACGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTACTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                                                                            >  1:802489/1‑143 (MQ=255)
       atgtctcTCGCCAAAGAGAGTGACGCTAATACTATCCCCAGCTATGGCGTAGTGTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCGGGTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                                                                         <  2:367895/137‑1 (MQ=255)
       aTGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAAACACATCAGACGTTGTCCAGCGCATg                                                                                                                                         >  2:331147/1‑143 (MQ=255)
             cTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAAc                                                                                                                                   >  2:1028131/1‑143 (MQ=255)
              tCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                                                                  <  1:794488/143‑1 (MQ=255)
                  cAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTGCTATGTGGAATGTCACTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGCCAACTTCTg                                                                                                                              <  1:1006562/143‑1 (MQ=255)
                          ttGCCCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTc                                                                                                                      <  2:527760/143‑1 (MQ=255)
                            gtcccTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGGCATGGGATTCAGTTTAGTGCTGCGACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAATTTCTGCCAGCTTCAc                                                                                                                    <  2:896514/141‑1 (MQ=255)
                             aCCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCCGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACg                                                                                                                   >  1:337975/1‑143 (MQ=255)
                             aCCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACg                                                                                                                   >  1:891669/1‑143 (MQ=255)
                             aCCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAAACAGATAAGAAAATCACCAGCGCAGGGCCAAATCCTGCAAGCTCCACg                                                                                                                   >  2:630545/1‑143 (MQ=255)
                               ccTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGaa                                                                                                                 >  1:311977/1‑143 (MQ=255)
                                cTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGCATTCAGTTTAGTGCTGCTACACACGGCGGGAATTGGTTAAACACAACATCAGAAGTTGCCCCGCGCATGGACCACTTCTGCCAACTTccacaac                                                                                                                >  2:349397/1‑137 (MQ=255)
                                 taataataaCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAg                                                                                                               >  1:259554/1‑143 (MQ=255)
                                   ataataaCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                                                             <  1:934694/143‑1 (MQ=255)
                                          ccccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCCGCTTCACGAAAGCCTTATACt                                                                                                      >  1:219585/1‑143 (MQ=255)
                                                 tATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCtttt                                                                                               <  1:284400/143‑1 (MQ=255)
                                                   tGGCGTACTAGGTGGCAGGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAACCAACATCAGACGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                                                                                             <  1:670821/143‑1 (MQ=255)
                                                     gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTGAAAACAACATCAGCAGTTGTCCAGCGCAGGGACAACTTCTGCGAGCTTCACGAAAGCCTTATTGTACCTTTTCGGc                                                                                           >  1:532309/1‑143 (MQ=255)
                                                        gaCTATGTGGAATGTCAGGGATGGGATTCAGTTTAGTGCTGCGACAAACGGCGGGAATTGGTTAAAAACAACAGCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGt                                                                                        <  2:883254/142‑1 (MQ=255)
                                                         aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtc                                                                                       <  1:286794/143‑1 (MQ=255)
                                                            gggtgGAACTTCAGTGATGGGTTTCAGTTTGGTGCTGCTACCAACTGCGGGAATTGGTTAAAAACAGCATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATGCTACCTTTTCGGCTGtttct                                                                                    <  2:274399/141‑1 (MQ=255)
                                                            aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGCTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctct                                                                                    <  2:802489/143‑1 (MQ=255)
                                                             tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctctt                                                                                   <  1:1028131/143‑1 (MQ=255)
                                                               tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                                                                 >  1:129688/1‑143 (MQ=255)
                                                                ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAACCCGCCGGAATTGGTTACAAACAACATCAGCAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATCCTACCTTTTCGCCTGtctcttctc                                                                                >  2:1039470/1‑143 (MQ=255)
                                                                 gAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctct                                                                               >  1:106873/1‑143 (MQ=255)
                                                                    tGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGACGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCGTATTCTACCTTGTCGGCTGTCCCGTCTCTAGt                                                                            >  2:232114/1‑143 (MQ=255)
                                                                      tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTATCTCGTAc                                                                          >  1:360030/1‑143 (MQ=255)
                                                                       cggtgtTGGTATTCAGTTTAGTGCTGCTACAAACGGCGGTAATTGGTTAAAAACTACATCAGAAGTTGTCCCGCGCATGGACATCTTCTGCCAGCTTCACGATAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACt                                                                         <  2:1006454/141‑1 (MQ=255)
                                                                        aGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                                                                        >  1:136238/1‑143 (MQ=255)
                                                                          ggATGGCATCCAATTTAGTGCGGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGCACTGtt                                                                      <  1:655895/142‑1 (MQ=255)
                                                                             taggAGTCAGTTTAGTGCTGCTACACACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                                                                   <  2:409583/141‑1 (MQ=255)
                                                                               ggATTCAGTTGAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGt                                                                 <  1:895555/143‑1 (MQ=255)
                                                                                 aTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                                                               >  1:465388/1‑143 (MQ=255)
                                                                                   tCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                                                             <  2:953348/143‑1 (MQ=255)
                                                                                    cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCCCCTCCCGTTTAAGTATTTg                                                            >  2:650723/1‑143 (MQ=255)
                                                                                          aGTGCTGCTACAAACGGCGGGAAGTGGTTAAAAACGTCTACGAAACGTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgc                                                      <  2:219585/143‑1 (MQ=255)
                                                                                            tgctgcTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc                                                    >  2:378062/1‑143 (MQ=255)
                                                                                             gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTCAGTATTTGTTCCGCGCt                                                   >  2:185293/1‑143 (MQ=255)
                                                                                              ctgctACAAAGCGCGGGAATTTGTTAAAAACAGCATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCtt                                                  <  2:1005038/143‑1 (MQ=255)
                                                                                                gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCa                                                >  2:976423/1‑143 (MQ=255)
                                                                                                      aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTCTAAGTATTTGTTCCGCGCTTGcagcagc                                          >  2:623447/1‑142 (MQ=255)
                                                                                                      aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTACTCGACCTTTTCGGCTGTCTCTGCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGt                                          <  2:682079/143‑1 (MQ=255)
                                                                                                             ggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTGCACGAAAGTCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGtctc                                   >  2:529648/1‑143 (MQ=255)
                                                                                                              gAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTctct                                  <  1:960075/143‑1 (MQ=255)
                                                                                                                 ttGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTGCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTGCGCGCTTCAAGCAGTGCGTTTCTGtg                               >  1:548319/1‑142 (MQ=255)
                                                                                                                     ttAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGACTCTCTGTCCTgg                           >  2:113605/1‑143 (MQ=255)
                                                                                                                      tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAACTTCACGAAAGACTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTATTGTTTAAGTATTTGTTCGGGCCTTGAAGCAGTGCTTCTCCGtacgtga                          >  2:247154/1‑137 (MQ=255)
                                                                                                                        aaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGa                        >  2:325386/1‑143 (MQ=255)
                                                                                                                         aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTGCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGaa                       <  1:850042/143‑1 (MQ=255)
                                                                                                                           acaacaTCACAAGTTGTCCAGCGCATGGACAACTTCCGCCAGCTTCACGAAAGCCTTATTCTACCTTGTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAAc                     <  2:149758/143‑1 (MQ=255)
                                                                                                                            caacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCACTGCGTCTCTGTCCTTGCGAAACg                    >  2:257524/1‑143 (MQ=255)
                                                                                                                              acaTCAGAAGTTGTCTAGAGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCGTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGcc                  <  1:498500/143‑1 (MQ=255)
                                                                                                                              acaTCAGAAGTTGCCCAGCGCATGGACATCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGcc                  <  1:219642/143‑1 (MQ=255)
                                                                                                                                agcAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCg                <  1:53460/141‑1 (MQ=255)
                                                                                                                                    gatgTTGTCCTCCGCATGGTCATCTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACaa            <  2:789594/140‑1 (MQ=255)
                                                                                                                                     aaGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAAt           >  2:1003794/1‑143 (MQ=255)
                                                                                                                                     aaGTTGTCCAGCGCATCGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAAt           <  1:190301/143‑1 (MQ=255)
                                                                                                                                       gTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGAGAATaa         >  1:252601/1‑143 (MQ=255)
                                                                                                                                           tCCAGCCCATGGACAACTTCTGCCAGCCTCACGAAAGCCTTATTCCACCTTTGCGGCTGTCTCTTCTCGCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGtt     <  1:2616/143‑1 (MQ=255)
                                                                                                                                             cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCGCGACAATAACGTTTc   >  1:344820/1‑143 (MQ=255)
                                                                                                                                              aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAGCGTTTCa  >  1:586463/1‑143 (MQ=255)
                                                                                                                                              aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCa  >  1:755761/1‑143 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCA  >  NC_000913/3183203‑3183487

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