Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,350 Δ1 bp 6.3% intergenic (+27/+31) 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,3500A. 6.3% 167.8 / 7.2 63intergenic (+27/+31)zupT/yqiDheavy metal divalent cation transporter ZupT/protein YqiD
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base . (2/2);  total (33/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 = 9.60e-02

ATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCA  >  NC_000913/3183210‑3183487
                                                                                                                                             |                                                                                                                                         
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)
                                                                                                                                             |                                                                                                                                         
ATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCA  >  NC_000913/3183210‑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: