Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 16.6% intergenic (‑497/‑85) 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,2620AC16.6% 518.5 / 103.7 218intergenic (‑497/‑85)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base A (104/78);  new base C (27/9);  total (131/87)
Fisher's exact test for biased strand distribution p-value = 6.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-02

Only 100 of 225 total aligned reads displayed.

GATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCT  >  NC_000913/4051164‑4051361
                                                                                                  |                                                                                                   
gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTaa                                                                                                     >  2:2112959/1‑99 (MQ=255)
       ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                             <  2:2296467/100‑1 (MQ=255)
       ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                             <  2:3682315/100‑1 (MQ=255)
       ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                             >  1:4415016/1‑100 (MQ=255)
        cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  2:3953062/100‑1 (MQ=255)
        cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                             <  2:3243540/99‑1 (MQ=255)
         ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAATAAAGGCGACAGATTACTCTGTCGCCt                                                                                           <  2:2685026/100‑1 (MQ=39)
         ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                           <  2:848957/100‑1 (MQ=255)
         ggCTAAAGGAAAAACAGGATGTGTTAGCGTACTGCAAGGATTGTAAGACCCGTTAAGGGTAAAAAAAAAGAAAAAAAGGCGACAGAGTAATCTGTCg                                                                                              <  1:1274921/97‑1 (MQ=255)
             aaaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttttt                                                                                       <  1:292278/100‑1 (MQ=255)
               aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                     >  1:1154454/1‑100 (MQ=255)
               aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                     >  1:2242802/1‑100 (MQ=255)
                ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                    >  2:1377337/1‑100 (MQ=255)
                ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                    >  2:1179246/1‑100 (MQ=255)
                ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                    >  1:3295165/1‑100 (MQ=255)
                ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                    >  2:4296225/1‑100 (MQ=255)
                 gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCtt                                                                                   <  1:1878539/100‑1 (MQ=255)
                  aaaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                                  <  1:79965/100‑1 (MQ=255)
                   aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                 >  1:3808170/1‑100 (MQ=255)
                   aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                 >  2:2931076/1‑100 (MQ=255)
                    aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                                <  1:1204379/100‑1 (MQ=255)
                         gggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTtgcttgct                                                                           >  1:290577/1‑100 (MQ=255)
                          ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                          >  1:1386504/1‑100 (MQ=255)
                            tgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                        <  1:919143/100‑1 (MQ=255)
                               gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                     >  2:1773684/1‑100 (MQ=255)
                               gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                     >  1:1798010/1‑100 (MQ=255)
                               gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                     >  2:218096/1‑100 (MQ=255)
                               gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                     <  1:3220363/100‑1 (MQ=255)
                                        cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                            >  1:4113075/1‑100 (MQ=255)
                                        cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                            >  1:1010466/1‑100 (MQ=255)
                                         tGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                           >  1:3992445/1‑100 (MQ=255)
                                         tGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                           <  1:1755889/100‑1 (MQ=255)
                                         tGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                           >  1:1986493/1‑100 (MQ=255)
                                          gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                          >  1:354339/1‑100 (MQ=255)
                                          gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                          >  1:1241883/1‑100 (MQ=255)
                                           cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                         >  2:1070359/1‑100 (MQ=255)
                                               gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTTTGTTAGATTCCGCCTCaaa                                                     >  1:4152329/1‑100 (MQ=255)
                                                aTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAAt                                                    >  2:3278841/1‑100 (MQ=255)
                                                 ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                   <  2:3811110/100‑1 (MQ=255)
                                                 ctGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                   <  1:2003900/99‑1 (MQ=255)
                                                   gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                                 >  1:3650245/1‑100 (MQ=255)
                                                    tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                                <  2:1275954/100‑1 (MQ=255)
                                                    tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                                <  2:1016357/100‑1 (MQ=255)
                                                      aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                              >  2:298490/1‑100 (MQ=255)
                                                       gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                             >  2:2911733/1‑100 (MQ=255)
                                                       gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                             >  1:2066799/1‑100 (MQ=255)
                                                       gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                             <  2:829218/100‑1 (MQ=255)
                                                        aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                            >  1:966828/1‑100 (MQ=255)
                                                        aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                            >  2:1881244/1‑100 (MQ=255)
                                                        aCCCATTAAGGGTTATTAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                            <  2:3462798/100‑1 (MQ=255)
                                                         cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTAGCTTGCTCTCTGTTAGATTCCCCCGCAAATCTATACTGa                                           >  1:3384559/1‑100 (MQ=255)
                                                            gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                        <  2:270173/100‑1 (MQ=255)
                                                              tAAGGGTTATGAGTCAGGAAACACGGCGACCGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaa                                      >  1:472672/1‑100 (MQ=255)
                                                               aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                     >  1:2221616/1‑100 (MQ=255)
                                                               aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                     >  1:3776867/1‑100 (MQ=255)
                                                                aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                    <  2:2078527/100‑1 (MQ=255)
                                                                   gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                                 <  2:3333903/100‑1 (MQ=255)
                                                                   gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCAGCCGCAAATCTATACTGAATAAAACGGc                                 >  2:2337089/1‑100 (MQ=255)
                                                                     tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                                 <  1:2672119/98‑1 (MQ=255)
                                                                     tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                               <  2:1540530/100‑1 (MQ=255)
                                                                          gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                          >  2:1457428/1‑100 (MQ=255)
                                                                           tCAGGAAAAGAGGCGACAGAGTAATCTGTCGCGTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAc                         <  1:1288639/100‑1 (MQ=255)
                                                                            cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                        >  1:4173795/1‑100 (MQ=255)
                                                                              ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGaa                      <  2:3642092/100‑1 (MQ=255)
                                                                              ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                       <  1:1421593/99‑1 (MQ=255)
                                                                              ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTATTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATGTATACTGAATAAAACGGCTAAAAGACGaa                      >  1:2648759/1‑100 (MQ=255)
                                                                               gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                     <  2:96184/100‑1 (MQ=255)
                                                                               gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTGTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                     >  1:2351103/1‑100 (MQ=255)
                                                                                aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                    >  2:1137750/1‑100 (MQ=255)
                                                                                aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                    <  1:2817769/100‑1 (MQ=255)
                                                                                aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                    <  1:388330/100‑1 (MQ=255)
                                                                                   aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                 >  2:357237/1‑100 (MQ=255)
                                                                                   aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                 >  2:2199444/1‑100 (MQ=255)
                                                                                   aaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGCATAAAACGGCTAAAAGACGAACCCtt                 >  1:1224072/1‑100 (MQ=255)
                                                                                   aaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                 >  2:4018646/1‑100 (MQ=255)
                                                                                    aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCTCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                >  1:4054947/1‑100 (MQ=38)
                                                                                    aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                <  2:3666424/100‑1 (MQ=255)
                                                                                     aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTTCGCCCCAAATTTATACTTAATAAACCGTCTAAAACACAACCCAttct               >  1:3700112/1‑98 (MQ=21)
                                                                                     aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  2:3773831/1‑100 (MQ=255)
                                                                                     aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  2:587625/1‑100 (MQ=255)
                                                                                     aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  1:2072969/1‑100 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATCCTGAATAAAACGGCCAAAAGACGAACCATTa                >  1:3246306/1‑98 (MQ=38)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  1:2605883/1‑100 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  1:163365/1‑100 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTTTTTTTTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  1:597169/1‑100 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              <  1:1100069/100‑1 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  2:1877299/1‑100 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              <  2:139996/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)
                                                                                       gCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGa             >  2:2778752/1‑100 (MQ=255)
                                                                                        cGACAGATTACTCTGTCGCTTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaa            >  1:1538889/1‑100 (MQ=255)
                                                                                          aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc          <  1:3051843/100‑1 (MQ=255)
                                                                                          aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc          >  2:534855/1‑100 (MQ=255)
                                                                                          aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc          <  1:155427/100‑1 (MQ=255)
                                                                                          aCAGAGTAATCTGTCGCCTTTTCTCTTTGCTAGCTTTGTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc          <  1:2967828/100‑1 (MQ=255)
                                                                                            agagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca        <  2:1939204/100‑1 (MQ=255)
                                                                                                 aaTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCatcatc   <  1:2279264/100‑1 (MQ=255)
                                                                                                  aTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCt  <  2:4075272/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCT  >  NC_000913/4051164‑4051361

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: