Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,350 Δ1 bp 15.2% intergenic (+27/+463) zupT → / ← ribB zinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,3500A.15.2% 600.5 / 107.5 204intergenic (+27/+463)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (99/74);  new base . (19/12);  total (118/86)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 207 total aligned reads displayed.

CTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTG  >  NC_000913/3183251‑3183448
                                                                                                    |                                                                                                  
cTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCa                                                                                                    <  1:1713197/100‑1 (MQ=255)
cTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCa                                                                                                    >  2:3390473/1‑100 (MQ=255)
  aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                                  <  1:960722/100‑1 (MQ=255)
  aTGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                                  <  1:2181478/100‑1 (MQ=255)
    ggCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATgg                                                                                                 <  1:3048724/99‑1 (MQ=255)
      cGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACa                                                                                              >  1:487049/1‑100 (MQ=255)
       gTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACaa                                                                                             >  2:2889516/1‑100 (MQ=255)
       gTACTATGTGGAACGTCAGTGACGGGATTCATTTTAGTGCTGCTACAAAGGGCGAAAAAAAAAAAAAAACAACATCAGAAGTTGTCCAGCGCATGGAc                                                                                               <  1:1692025/98‑1 (MQ=255)
         aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGTTGCTACAAACGGCGGAAATCAGAAAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAAc                                                                                            <  1:2102583/99‑1 (MQ=25)
         aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                           <  1:713794/100‑1 (MQ=39)
         aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACt                                                                                           <  1:2845132/100‑1 (MQ=39)
         aCTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACt                                                                                           >  2:1927681/1‑100 (MQ=255)
           tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTc                                                                                         >  1:458071/1‑100 (MQ=255)
            aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGTGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCt                                                                                        <  2:131924/100‑1 (MQ=38)
            aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                        >  1:3007460/1‑100 (MQ=255)
             tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCt                                                                                        <  2:3286605/99‑1 (MQ=38)
              gtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGATTAGAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCt                                                                                        >  1:3569499/1‑98 (MQ=255)
              gtgGAATGTCAGTGATGGGATTCAGGTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAGGTTGTCCAGCGCATGGACAACTTCTGc                                                                                      <  1:1132041/100‑1 (MQ=255)
              gtgCAATGTAAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTg                                                                                       <  1:2943773/99‑1 (MQ=37)
               tgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                     <  1:2512615/100‑1 (MQ=255)
               tgGAATGTCAGTGATGGGATTCAGTGTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                     <  1:84920/100‑1 (MQ=255)
                ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgatg                                                                                    <  2:2987449/100‑4 (MQ=38)
                ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCa                                                                                    >  2:3006330/1‑100 (MQ=255)
                ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCa                                                                                    <  2:2052747/100‑1 (MQ=255)
                  aaTGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGc                                                                                  <  2:260432/100‑1 (MQ=255)
                    tataaatGACGGGATTCAGTATTATACTGCTACATAAGGATGGAATTAAAGACAAACAAAAATATAATTTGTCCGGCGCATGGACAACTGCTGCCAg                                                                                   >  1:1634616/7‑97 (MQ=255)
                       cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCAc                                                                             >  1:1267218/1‑100 (MQ=255)
                       cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCAc                                                                             <  2:686172/100‑1 (MQ=255)
                       cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCAc                                                                             >  1:195580/1‑100 (MQ=255)
                         gTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGa                                                                           >  1:230910/1‑100 (MQ=255)
                            aTGGGAGTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAg                                                                        <  2:2399810/100‑1 (MQ=255)
                             tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGc                                                                       >  2:3143430/1‑100 (MQ=255)
                              gggATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                      <  2:3816713/100‑1 (MQ=255)
                               ngATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCt                                                                     >  2:3310371/2‑100 (MQ=255)
                                gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCCGCGCATGGACAACTGCTGCCAGCTTCACGAAAGCCtt                                                                    >  1:2595743/1‑100 (MQ=255)
                                gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                    >  2:482586/1‑100 (MQ=255)
                                gATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                    >  2:305108/1‑100 (MQ=255)
                                 aTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTa                                                                   <  1:2734192/100‑1 (MQ=255)
                                   tcagggtAGTGCGGCGATAAACGGAGGTGTTTGGTTAAAAACAGTATCAGAAGTAGTCCACGGCATGGGCAACTTCTGCCGGCTTCACGAAAGCCTTAtt                                                                 <  1:3349851/94‑1 (MQ=255)
                                    cAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACATCTGCCAGCTTAACGAAAGCCTTATTc                                                                >  1:220909/1‑100 (MQ=255)
                                     aGTTTAGTGCTGCTACAAACGGCGGGCATTGGTTAAACACAACATCAGCAGTTGTCCAGCGCATGGACAGCTTCTGGCAGCTTCACGATAGCCTTAtt                                                                 <  1:3744436/98‑1 (MQ=255)
                                      gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTa                                                              >  2:731371/1‑100 (MQ=255)
                                        ttAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                            <  2:3108944/100‑1 (MQ=255)
                                         tAGTGCTGCTACAAACGGCTGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGTACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAcc                                                            <  2:3642697/99‑1 (MQ=255)
                                             gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                       <  2:1257082/100‑1 (MQ=255)
                                             gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGAACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                       >  2:2119693/1‑100 (MQ=255)
                                               tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                     <  1:2226696/100‑1 (MQ=255)
                                               tgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCgg                                                     >  1:3332997/1‑100 (MQ=255)
                                                gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                    >  2:409469/1‑100 (MQ=255)
                                                gctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                    <  2:2312662/100‑1 (MQ=255)
                                                gctACAAACGGCGGGAATTGGTCAAAAACAACATCAGAAGTTGTACAGCGCATGGACACCTTCTGCCCGCTACACGAAAGCCTTATTCTACCTTTTCGGc                                                    <  1:3658552/100‑1 (MQ=255)
                                                  tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                                  >  1:1464853/1‑100 (MQ=255)
                                                    cAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtc                                                >  2:3518404/1‑100 (MQ=255)
                                                     aaaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtct                                               >  1:2765134/1‑100 (MQ=255)
                                                          gcggGAATTGGTTAAAAACAACATCAGAAGTTGTCGAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                           <  1:2130905/99‑1 (MQ=255)
                                                          gcggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttct                                          >  1:3538317/1‑100 (MQ=255)
                                                            gggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTATACCTTTTCGGCTGtctcttctct                                        >  1:1902922/1‑100 (MQ=255)
                                                              gAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                      <  1:1061325/100‑1 (MQ=255)
                                                              gAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                      >  2:1061333/1‑100 (MQ=255)
                                                                   ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                                 <  1:3536714/100‑1 (MQ=255)
                                                                   ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTg                                 <  2:1963602/100‑1 (MQ=255)
                                                                    gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGt                                >  2:1837096/1‑100 (MQ=255)
                                                                     ttAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGtt                               >  2:1648970/1‑100 (MQ=255)
                                                                     ttAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGtt                               >  2:3243691/1‑100 (MQ=255)
                                                                      tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                              >  2:84449/1‑100 (MQ=255)
                                                                      tAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGttt                              >  1:1022466/1‑100 (MQ=255)
                                                                       aaaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                             >  2:1209781/1‑100 (MQ=255)
                                                                         aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                           <  2:1899393/100‑1 (MQ=255)
                                                                         aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                           <  2:3556376/100‑1 (MQ=255)
                                                                              ctggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAttt                      >  2:449055/5‑100 (MQ=37)
                                                                              acaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                        <  1:1920742/98‑1 (MQ=255)
                                                                               caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                     >  2:792266/1‑100 (MQ=255)
                                                                               caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                     <  1:404459/100‑1 (MQ=255)
                                                                                ggcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                    >  2:2480245/3‑100 (MQ=38)
                                                                                aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                    >  2:658505/1‑100 (MQ=255)
                                                                                aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGt                    >  2:1178913/1‑100 (MQ=255)
                                                                                aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTCCTCTCGTACTGTTTAAGTATTTGt                    >  1:342649/1‑100 (MQ=255)
                                                                                 tCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGtt                   >  2:1467555/1‑100 (MQ=255)
                                                                                  cAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc                  >  1:3298902/1‑100 (MQ=38)
                                                                                  cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACGTTTTCGTCTGTCTCTTCTCTCGTAATGTTTAAGTATTTGTTc                  >  1:3193471/1‑100 (MQ=255)
                                                                                  cAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTc                  >  1:2876207/1‑100 (MQ=255)
                                                                                   aGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTcc                 >  2:1357124/1‑100 (MQ=38)
                                                                                    gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg                >  2:891348/1‑100 (MQ=39)
                                                                                    gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg                >  2:26469/1‑100 (MQ=39)
                                                                                    gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg                >  1:1410709/1‑100 (MQ=39)
                                                                                    gAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCACTCATACTGTTTAAGTATTTCtt                   >  1:1458425/1‑97 (MQ=37)
                                                                                    gAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg                >  1:1768226/1‑100 (MQ=255)
                                                                                    gAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg                >  2:2223076/1‑100 (MQ=255)
                                                                                    gAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCGTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg                >  1:2468301/1‑100 (MQ=255)
                                                                                      aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcgcg              >  2:2583710/1‑100 (MQ=39)
                                                                                       gTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGACAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc             >  1:2374296/1‑100 (MQ=39)
                                                                                        ttGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCt            <  2:220913/100‑1 (MQ=255)
                                                                                          gTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTc          >  2:1036083/1‑100 (MQ=255)
                                                                                           tCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCa         >  2:2011727/1‑100 (MQ=255)
                                                                                           tCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAGTTGTTCCGCGCTTCa         <  1:2448821/100‑1 (MQ=255)
                                                                                             cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg       >  2:592042/1‑100 (MQ=255)
                                                                                             cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAg       >  1:241676/1‑100 (MQ=255)
                                                                                              aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc      >  1:3220949/1‑100 (MQ=255)
                                                                                                   cATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTg  >  1:3483003/1‑100 (MQ=255)
                                                                                                   cATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTg  >  2:2170552/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
CTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTG  >  NC_000913/3183251‑3183448

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