Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 16.7% 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,2590GT16.7% 578.8 / 93.6 215intergenic (‑494/‑88)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base G (101/78);  new base T (27/9);  total (128/87)
Fisher's exact test for biased strand distribution p-value = 4.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 221 total aligned reads displayed.

ACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051161‑4051358
                                                                                                  |                                                                                                   
   gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                 >  2:3924328/1‑100 (MQ=255)
       gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGGTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTg                                                                                              >  2:1310521/1‑99 (MQ=255)
        ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTc                                                                                            <  1:3385615/100‑1 (MQ=255)
          ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                          <  2:3682315/100‑1 (MQ=255)
          ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                          >  1:4415016/1‑100 (MQ=255)
           cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                         <  1:1847909/100‑1 (MQ=255)
            ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                        <  2:848957/100‑1 (MQ=255)
             gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtt                                                                                       >  2:1878546/1‑100 (MQ=255)
                  aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                  >  1:1154454/1‑100 (MQ=255)
                  aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                  >  1:2242802/1‑100 (MQ=255)
                   ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  2:4296225/1‑100 (MQ=255)
                   ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  1:2130909/1‑100 (MQ=255)
                   ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  2:1179246/1‑100 (MQ=255)
                      aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  1:3808170/1‑100 (MQ=255)
                      aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  2:2931076/1‑100 (MQ=255)
                       aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                             <  1:1204379/100‑1 (MQ=255)
                          cAGCGTGTTTTGGCGGCCTGGAAGGATTGTAAGACCTTTTAAGGGTTATGAGTCAGTTAAAAAGGCTACAGCGTAATCTGTCGCCTTTTTTCTttgc                                                                             <  1:1568675/97‑1 (MQ=255)
                           aGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttg                                                                          <  2:3767019/99‑1 (MQ=255)
                             ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                       >  1:1386504/1‑100 (MQ=255)
                               tgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                     <  1:919143/100‑1 (MQ=255)
                                 tgtTGGCGGTCTGCAAGGAATGTAAGACCCTGTGAGGATTGTGCGTTAGAAAAGATGAAGACAGAGTAATCTGTGGACTTTCATCTCTgattcctt                                                                       >  1:3098463/1‑96 (MQ=255)
                                  gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                  >  2:219008/1‑100 (MQ=255)
                                  gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                  <  1:3220363/100‑1 (MQ=255)
                                    tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                <  2:1831236/100‑1 (MQ=255)
                                      gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                              <  1:1498957/100‑1 (MQ=255)
                                        ggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                <  1:3725601/96‑1 (MQ=255)
                                           cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                         >  1:4113075/1‑100 (MQ=255)
                                           cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                         >  1:1010466/1‑100 (MQ=255)
                                            tGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                        <  1:1755889/100‑1 (MQ=255)
                                            tGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                        >  1:1986493/1‑100 (MQ=255)
                                             gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       >  1:1241883/1‑100 (MQ=255)
                                             gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       >  1:354339/1‑100 (MQ=255)
                                                 ggATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaa                                                   >  1:728140/1‑100 (MQ=255)
                                                  gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTTTGTTAGATTCCGCCTCaaa                                                  >  1:4152329/1‑100 (MQ=255)
                                                  gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                  <  1:1215843/100‑1 (MQ=255)
                                                   aTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAAt                                                 >  2:3278841/1‑100 (MQ=255)
                                                    ctGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  1:2003900/99‑1 (MQ=255)
                                                     tGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCt                                               <  2:1052041/100‑1 (MQ=255)
                                                      gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                              >  1:3650245/1‑100 (MQ=255)
                                                       tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                             <  2:1275954/100‑1 (MQ=255)
                                                       tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                             <  2:1016357/100‑1 (MQ=255)
                                                         aGACCCGTTAAGGGTTATGAGTCAGGAAAAACGGCGACAGGGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCTCATCTacac                                           >  1:132112/1‑100 (MQ=255)
                                                         aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                           >  2:298490/1‑100 (MQ=255)
                                                          gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          >  1:2066799/1‑100 (MQ=255)
                                                          gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          <  2:829218/100‑1 (MQ=255)
                                                          gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          >  2:2911733/1‑100 (MQ=255)
                                                          gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          >  2:2109350/1‑100 (MQ=255)
                                                          gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          <  2:4048166/100‑1 (MQ=255)
                                                           aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                         <  2:2297331/100‑1 (MQ=255)
                                                           aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                         >  1:966828/1‑100 (MQ=255)
                                                           aCCCATTAAGGGTTATTAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                         <  2:3462798/100‑1 (MQ=255)
                                                             ccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGaa                                       >  2:3784967/1‑100 (MQ=255)
                                                              cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAt                                      <  2:212148/100‑1 (MQ=255)
                                                               gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                     >  1:3360470/1‑100 (MQ=255)
                                                                 tAAGGGTTATGAGTCAGGAAACACGGCGACCGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaa                                   >  1:472672/1‑100 (MQ=255)
                                                                  aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                  >  1:3776867/1‑100 (MQ=255)
                                                                  aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                  >  1:2221616/1‑100 (MQ=255)
                                                                   aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                 <  2:2078527/100‑1 (MQ=255)
                                                                      gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                              <  2:3333903/100‑1 (MQ=255)
                                                                      gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                              <  2:3523230/100‑1 (MQ=255)
                                                                        tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                              <  1:2672119/98‑1 (MQ=255)
                                                                        tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                            <  2:1540530/100‑1 (MQ=255)
                                                                             gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                       >  2:1457428/1‑100 (MQ=255)
                                                                               cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCGGCCGCCAATCTATAGTGAATAAAACGGCTAAAAGACg                     >  1:2979983/1‑100 (MQ=255)
                                                                               cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                     >  1:4173795/1‑100 (MQ=255)
                                                                                aagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                    >  2:3240906/3‑100 (MQ=255)
                                                                                aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                    >  2:1676662/1‑100 (MQ=255)
                                                                                 ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAACGACGaa                   >  1:1288397/1‑100 (MQ=255)
                                                                                 ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCAGCAAATCTATACTGAATAAAACGGCTAAAAGACGaa                   >  1:1549892/1‑100 (MQ=255)
                                                                                 ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTATTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATGTATACTGAATAAAACGGCTAAAAGACGaa                   >  1:2648759/1‑100 (MQ=255)
                                                                                  gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                  <  2:96184/100‑1 (MQ=255)
                                                                                  gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                  >  2:1100072/1‑100 (MQ=255)
                                                                                  gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTGTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                  >  1:2351103/1‑100 (MQ=255)
                                                                                   aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                 <  1:388330/100‑1 (MQ=255)
                                                                                   aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                 <  1:2817769/100‑1 (MQ=255)
                                                                                     aaaaGGCGACAGAGTAATCTGTCGCCTTTTTTTTTTTTTTTTTTTTTGTTAGATTCCGCCGCAAATCTTTTCTGTATAAAACGGCTAAAAGTCGAAcctt               >  1:2836274/1‑98 (MQ=255)
                                                                                      aaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt              <  2:266839/100‑1 (MQ=255)
                                                                                      aaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt              >  2:4018646/1‑100 (MQ=255)
                                                                                       aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCTCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             >  1:4054947/1‑100 (MQ=38)
                                                                                       aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             <  2:3666424/100‑1 (MQ=255)
                                                                                        aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  2:587625/1‑100 (MQ=255)
                                                                                         ggtgACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:1254690/4‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATCCTGAATAAAACGGCCAAAAGACGAACCATTa             >  1:3246306/1‑98 (MQ=38)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:2605883/1‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:2961154/1‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:163365/1‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:2274478/1‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:2726998/1‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:2699896/1‑100 (MQ=255)
                                                                                         ggCGACAGAGTAATCTGTCGCCTTTTTTTTTTTTTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:597169/1‑100 (MQ=255)
                                                                                         ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           <  1:1100069/100‑1 (MQ=255)
                                                                                          gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:3475370/1‑99 (MQ=255)
                                                                                          gCGACAGATTACTCTGTCGCCTTTTTTATTTGCTTGCTTTCTGTTAGATTCCGCCGCAGATCTATACTGAATAAAACGGCTAAAAGACGAACAGTTATg           >  1:552042/1‑99 (MQ=38)
                                                                                          gCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTCTGa          >  2:564937/1‑100 (MQ=255)
                                                                                             aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc       <  2:1380224/100‑1 (MQ=255)
                                                                                             aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc       >  2:534855/1‑100 (MQ=255)
                                                                                              cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAAGAACCATTATGAAAc       >  1:4031916/1‑99 (MQ=255)
                                                                                               agagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca     <  2:1939204/100‑1 (MQ=255)
                                                                                                gagTAATCTGTCGCCTTTTTTCTTTGTTTGCTTTCTGTTAGATTCCGCCGCAAATCTATTCTGAATAAAACGGCTAAAAGACGAACCCTTCTGAAACcat    >  1:1962403/1‑100 (MQ=255)
                                                                                                  gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  >  2:360277/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051161‑4051358

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: