Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 15.4% intergenic (‑494/‑88) yihA ← / → yihI cell division GTP‑binding protein/activator of Der GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,2590GT15.4% 442.5 / 59.2 161intergenic (‑494/‑88)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base G (77/59);  new base T (20/5);  total (97/64)
Fisher's exact test for biased strand distribution p-value = 4.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 168 total aligned reads displayed.

CACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051160‑4051358
                                                                                                   |                                                                                                   
cacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTATTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACagag                                                                                                     >  1:1857868/1‑100 (MQ=255)
  cGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTa                                                                                                   >  2:1664216/1‑100 (MQ=255)
    gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                 >  2:1286882/1‑100 (MQ=255)
    gATTGTTCCGGAGAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTAAAAAAAAAAAAAAAAGGCGACAGAGt                                                                                                    <  1:1260487/97‑1 (MQ=255)
     aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                <  2:664866/100‑1 (MQ=255)
        gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                             >  1:790586/1‑100 (MQ=255)
        gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                             >  1:1691516/1‑100 (MQ=255)
          tCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCg                                                                                           >  2:2916194/1‑100 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  2:2134805/100‑1 (MQ=255)
             ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                        >  1:2608243/1‑100 (MQ=255)
              gatAAAGGAAAAACAGGGTGTGTTGGCGGAATGCAAGGATTGTAGGGCCCGTTAAGGGTAAAAAAAAAAAAAAAAAGGCGACAGAGTAATCTGTCGCCtt                                                                                       <  1:1800583/98‑1 (MQ=255)
               cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttt                                                                                      <  2:2608256/100‑1 (MQ=255)
               cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttt                                                                                      >  2:680301/1‑100 (MQ=255)
               cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttt                                                                                      <  1:435640/100‑1 (MQ=255)
                 aaaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttttt                                                                                    <  2:840170/100‑1 (MQ=255)
                 aaaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttttt                                                                                    >  2:3051863/1‑100 (MQ=255)
                 aaaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttttt                                                                                    <  2:3275672/100‑1 (MQ=255)
                  aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAGGACCGGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  1:85316/94‑1 (MQ=38)
                   aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttttt                                                                                   >  1:3399386/1‑99 (MQ=255)
                    ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGAACCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  1:2100499/1‑100 (MQ=255)
                      aaaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                               <  2:3552871/100‑1 (MQ=255)
                       aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  2:1907312/1‑100 (MQ=255)
                        aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                             <  1:1302092/100‑1 (MQ=255)
                           cAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttg                                                                          <  2:238106/100‑1 (MQ=255)
                              ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                       >  2:1458082/1‑100 (MQ=255)
                                tgtgtTGGCGGCCTGCAAGGATTGTAAGGCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                     <  1:3724566/100‑1 (MQ=255)
                                   gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                  >  1:3144303/1‑100 (MQ=255)
                                   gtTGGCGGCCTGCAAGGATTGTAAGACCCCTTAAGCGTTATGAGACAGGAAAAAAGGCGACAGAGTCATCTGTCGCCTTTTTTCTTTGCTTCCTTTCTg                                                                   >  1:3324258/1‑99 (MQ=38)
                                    ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAGGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                 >  1:2739202/1‑100 (MQ=255)
                                      ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                >  1:2248669/1‑99 (MQ=255)
                                       gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                              <  1:3139630/100‑1 (MQ=255)
                                          gcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTc                                                           >  2:318061/1‑100 (MQ=255)
                                           ccTGCAAGGATTGTAAGCCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTcc                                                          <  1:3123786/100‑1 (MQ=255)
                                             tGCAAGGATTGTAAGACCCGTTAAGGTTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                        <  1:638127/100‑1 (MQ=255)
                                              gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       <  2:1857877/100‑1 (MQ=255)
                                              gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       >  1:3510057/1‑100 (MQ=255)
                                              gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       <  2:2284882/100‑1 (MQ=255)
                                               cAAGGATTGTCAGACCCGTTAAGGGTTATGCGTCAGGAAAAAAGGCGACAGAGTAATCTCTCGCCTTTTTTATTTGCTTGCTTCCTGTTAGATTccgccg                                                      <  1:1123754/100‑1 (MQ=38)
                                               cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                      >  2:1525545/1‑100 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                    <  2:3458447/100‑1 (MQ=255)
                                                    aTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAAt                                                 >  2:835897/1‑100 (MQ=255)
                                                     ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  2:2112176/100‑1 (MQ=255)
                                                     ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  2:3513854/100‑1 (MQ=255)
                                                     ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  1:558951/100‑1 (MQ=255)
                                                     ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  1:1202797/100‑1 (MQ=255)
                                                      tGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCt                                               >  1:3746158/1‑100 (MQ=255)
                                                       gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAGAt                                                 >  1:2484120/1‑97 (MQ=255)
                                                       gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCt                                               >  1:1622528/1‑99 (MQ=255)
                                                       gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCACAAATCta                                              <  1:2964246/100‑1 (MQ=255)
                                                        tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                             <  2:1623508/100‑1 (MQ=255)
                                                        tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                             >  2:1543702/1‑100 (MQ=255)
                                                           gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          >  1:147003/1‑100 (MQ=255)
                                                            aCCCGTTAAGGGTTATGAGTCAGGAAAAAGGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                         <  2:2739216/100‑1 (MQ=255)
                                                              ccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtagatcgga                                       >  1:1966998/1‑93 (MQ=255)
                                                                gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                     <  2:3333539/100‑1 (MQ=255)
                                                                gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                     >  2:2226116/1‑100 (MQ=255)
                                                                gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                     >  2:989984/1‑100 (MQ=255)
                                                                    aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                 <  1:3051846/100‑1 (MQ=255)
                                                                      ggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                               >  1:2765128/1‑100 (MQ=255)
                                                                       gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                              <  2:39448/100‑1 (MQ=255)
                                                                         tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                            >  2:1665265/1‑100 (MQ=255)
                                                                           tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCTAATCTATATTGAATAAAACGGCTaaa                          >  1:3706415/1‑100 (MQ=255)
                                                                           tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                          <  2:3405711/100‑1 (MQ=255)
                                                                            gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                         >  1:1780606/1‑100 (MQ=255)
                                                                              gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTGTCTTTGCTTGCTTTCTGTTAGATTCTGACGCAAATCTATACTGAATAAAACGGCTaaaa                         >  1:1040527/1‑98 (MQ=255)
                                                                               caaagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAc                      >  1:635120/5‑100 (MQ=37)
                                                                               caaagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCAGCAAATATATACTGAATAAAACGGCTAAAAGAc                      >  1:3538647/5‑100 (MQ=25)
                                                                                aaagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCCCCGCAAATCTATACTGAATAAAACGGCTAAAAGAc                      >  1:1840591/4‑99 (MQ=37)
                                                                                 aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                    >  1:2691247/1‑100 (MQ=255)
                                                                                 aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGATAGATTCCGCCGCAAATTTATACTGAATAAAACGGCTAAAAGACGa                    >  1:1386041/1‑100 (MQ=255)
                                                                                  ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                    >  2:3202094/1‑99 (MQ=255)
                                                                                  agAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCATTCTGTTAGATTCCCCCGCAAATCTATACTCCATACAACGGCTaca                          >  1:3677837/2‑91 (MQ=21)
                                                                                   gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                  <  2:2888992/100‑1 (MQ=255)
                                                                                    aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGTAAATCTATACTGAATAAAACGGCTAAAAGACGaa                   >  1:905593/1‑98 (MQ=255)
                                                                                    aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                 >  1:3358414/1‑100 (MQ=255)
                                                                                     agaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  2:3534539/3‑100 (MQ=255)
                                                                                     aaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAACTCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  1:2003906/1‑100 (MQ=255)
                                                                                        aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             <  1:2226106/100‑1 (MQ=255)
                                                                                        aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             >  1:3579055/1‑100 (MQ=255)
                                                                                         aGGCGGCAGATTACTCTGTCGCCTTTTTTCTTTGATTGCTTTCTGTGAGATTCCCCCTGTATTTTATACTTAATGAGACGGCTAAGAGACGCACCATTa             >  1:754429/1‑99 (MQ=16)
                                                                                         aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  2:728576/1‑100 (MQ=255)
                                                                                         aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTTCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGTATAGAACGGCTAAAAGACGAACCATTAt            >  1:2870460/1‑100 (MQ=255)
                                                                                         aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  2:890702/1‑100 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTGTGTCAGATTCTGCTGCAAATCTCTACTGAATTAAACGGCTAAAAGCCGAACCATCATg           >  1:2541063/1‑100 (MQ=25)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  2:2654761/1‑99 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:381499/1‑100 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:1002325/1‑100 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:2795486/1‑100 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:1157597/1‑100 (MQ=255)
                                                                                          ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:464347/1‑100 (MQ=255)
                                                                                           gCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGa          >  2:1205402/1‑100 (MQ=255)
                                                                                            cGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaa         >  2:1897402/1‑100 (MQ=255)
                                                                                             gACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa        >  2:1226714/1‑100 (MQ=255)
                                                                                              aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc       >  2:2056807/1‑100 (MQ=255)
                                                                                               cAGAGTAATCTGTCGCCTTTTTTTTTTTTTTTTTTTCTGATAGAATCCGCCGCAAATCTATACTGAATAAAATGGCTAAAAGACGAACCATTATGAAAcc      >  1:3207880/1‑100 (MQ=255)
                                                                                               cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAcc      >  1:986707/1‑100 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    <  2:2763670/100‑1 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    >  2:1172792/1‑100 (MQ=255)
                                                                                                   gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  <  2:2775768/100‑1 (MQ=255)
                                                                                                   gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  <  1:1907303/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051160‑4051358

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