Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,893 C→T 75.1% R46R (CGC→CGT ynfO → uncharacterized protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,8930CT75.1% 49.8 / 77.3 177R46R (CGC→CGTynfOuncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (19/24);  new base T (61/67);  total (83/94)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

Only 100 of 178 total aligned reads displayed.

GGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCCCCTCATATCACATGG  >  NC_000913/1636747‑1637035
                                                                                                                                                  |                                                                                                                                              
ggCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTACCCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCCTTGAAGCCTGccggcg                                                                                                                                              >  1:54863/1‑149 (MQ=12)
  catccTCATGTCAACGAAGTACAGAACCCGCAGAACATCAAGCCGCAACTTCCGATTTCCTAACCAAATGATTTAACTATTTAACATCGCTGTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGga                                                                                                                                            <  2:98314/146‑1 (MQ=1)
    aCCTCATGTCAACGACGAACAGAACCCGCAGAACCACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  6:209721/149‑1 (MQ=2)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACACCCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          >  2:73132/1‑149 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACACACCGCAACCTCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  6:296289/149‑1 (MQ=2)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGGAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  7:191474/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAACAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  6:312604/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  6:27979/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  1:384695/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  3:353705/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  7:107926/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  7:310604/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  4:284378/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  7:395150/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          >  2:219223/1‑149 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  7:78261/149‑1 (MQ=11)
    aCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  8:245271/149‑1 (MQ=11)
     ccTCATGTCAACGAAGAACAGACCCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  3:350162/1‑149 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          >  7:213185/1‑148 (MQ=11)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  8:213185/148‑1 (MQ=11)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         <  7:255149/149‑1 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         <  2:48803/149‑1 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  4:271583/1‑149 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         <  8:219807/149‑1 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         <  3:334392/149‑1 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         <  3:313900/149‑1 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  5:292939/1‑149 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  5:312481/1‑149 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         <  4:91595/149‑1 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  1:30012/1‑149 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         <  8:395408/149‑1 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGAc                                                                                                                                         <  5:170756/149‑1 (MQ=35)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGACGTCGGAGAc                                                                                                                                         >  1:259013/1‑149 (MQ=14)
     ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGAGTGAACAAATTAACATCGCACCTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGACGCCTGCCGACGGAGAc                                                                                                                                         >  2:361125/1‑149 (MQ=11)
      ctaaTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAACTTAACATCGCTCTTTAGCAAAAAGGGGCCGTGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACt                                                                                                                                        <  5:90057/146‑1 (MQ=11)
      cTCATGTCAACGAAGAACAGAACCCGCAGAACACCAACCCGCCACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  5:2475/1‑148 (MQ=9)
      cTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAACAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  5:351300/1‑148 (MQ=14)
      cTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  4:197339/1‑148 (MQ=14)
      cTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  6:135527/1‑148 (MQ=14)
      cTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  4:67500/1‑148 (MQ=14)
      cTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  5:225209/1‑148 (MQ=14)
      cTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCACTTGAGCAAAAAGGTTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGAc                                                                                                                                         >  7:148928/1‑148 (MQ=9)
       tCATGTCAACGGAGAACAGATCCCGCAGAACAACCCCCGGCACCATCCGCTTTCCTAACCAAATGTATGAACAATTTAACATCGCGCTTGAGCTAAAAGGGGGAGGGACAGGCGCAGCCTCGCACATTGAAGCCTGCCGTCGGAGAcg                                                                                                                                        >  4:13311/1‑147 (MQ=2)
       tCATGTCAACGAAGAACAGACCCCGCAGAACAACAACCCGCGACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACt                                                                                                                                        >  8:369858/1‑148 (MQ=12)
       tCATGTCAACGAAGAACAGAACCCGCAGAACAACACCCCGCAACATCCGCTTTCCTAACCAAATGAATGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGATTTTGTCAGCCTGGGTCATTGAAGCCTGCCGTCGCAGACt                                                                                                                                        >  4:225690/1‑148 (MQ=11)
       tCATGTCAACGAAGAACAGAACCCGCAGAACAACACCCCGCAACATCCCATTTCCTAACCAAATGATTGAACAAATTAACATCCATCTTCAGCAAAAAGGGTCACGCAATTTCTCAGACTGGGTCATTCAAGCCTGACGTCGgagaat                                                                                                                                        >  1:242153/1‑146 (MQ=11)
       tCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          >  5:306894/1‑146 (MQ=11)
       tCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGgaga                                                                                                                                          <  6:306894/146‑1 (MQ=11)
       tCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACt                                                                                                                                        >  3:200223/1‑148 (MQ=17)
       tCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACt                                                                                                                                        >  6:204011/1‑148 (MQ=17)
       tCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACt                                                                                                                                        >  1:84337/1‑148 (MQ=17)
        cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      >  3:298723/1‑147 (MQ=11)
        cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      >  3:264677/1‑147 (MQ=11)
        cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      >  2:298854/1‑147 (MQ=11)
        cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      >  1:148160/1‑147 (MQ=11)
        cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      >  4:360830/1‑147 (MQ=2)
        cATGTCAACGAACAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGGCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      <  1:117018/149‑3 (MQ=2)
         aTGTCACCGAAGAACAGAACCCGCAGCACAACAACCCGCAACATCCGCTTTCCTAACCAATTGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCGGGGAATTTGTCAGCCTGGGTCATTGAAGCCTGCCGACGGAGACtaac                                                                                                                                     >  8:257140/1‑146 (MQ=11)
         aTGTCACCGAAGAACAGAACCCGCAGAACAACCCCCCGCAACATCCGCTTTCCTAACCAACTGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCCTTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  5:278913/149‑4 (MQ=1)
         aTGTCACCGAAGAACAGAACCAGCAGAAAAACAACCCGAAACATCCGATTTACTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCAGGGAATTTCTCAGCCTGTGTCATTGAAGCATGCCGTCTGAGACtaa                                                                                                                                      <  8:352365/148‑3 (MQ=11)
         aTGTCAACGAAGAACAGCACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACAACGATATAGAGCAAAAAGGGACCGGGAATTTCTCAGCATGGGTAATTGAAGCCGGCCGTCGGAGACtaac                                                                                                                                     >  3:118964/1‑146 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGACCAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGACTTTCTCAGCCTGGGTCCTTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  2:228722/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAACTTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      <  7:284859/148‑3 (MQ=2)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGATCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  2:126615/149‑4 (MQ=0)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     >  8:93891/1‑146 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  8:140013/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     >  1:115946/1‑146 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  5:169781/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  2:281294/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  1:345145/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  4:64715/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  3:74416/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  4:159443/149‑4 (MQ=1)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaa                                                                                                                                      <  3:164705/148‑3 (MQ=11)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACta                                                                                                                                       <  2:327377/147‑2 (MQ=14)
         aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACta                                                                                                                                       <  5:106542/147‑2 (MQ=14)
         aTGTCAACGAACAACAGAACCCGCAGAACAACAACCCGCACCATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACtaac                                                                                                                                     <  2:253106/149‑4 (MQ=255)
           gTCAACGAAGAACAGAACCCGCAGCACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTAATTGAAGCCTGCCGGCGGAGACtaacg                                                                                                                                    >  8:313084/1‑144 (MQ=255)
            tCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  4:287394/146‑4 (MQ=1)
            tCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     >  3:144763/1‑143 (MQ=1)
            tCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGCTCATTGAAGCCTGCCGTCGGAGAcaaacgt                                                                                                                                   >  3:144893/1‑142 (MQ=0)
            tCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCGATTGAGCAAAAAGGGACCGGGAATTTCTCAGCATGGGAAATTGAAGCATGCCGTCGGTGACtaactt                                                                                                                                   >  6:176974/1‑143 (MQ=1)
            tCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCATAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaac                                                                                                                                     <  4:144763/146‑4 (MQ=0)
              aaCGAAGAACAGACCCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAAATTCTCAGCCTGGGTCATTGAAGCCTGCCGTCGGAGACtaactt                                                                                                                                   >  4:5558/1‑141 (MQ=0)
                            cccGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGGCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTttct                                                                                                                  <  5:126514/149‑1 (MQ=255)
                                  gAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGaaaaaa                                                                                                                          <  2:148963/135‑1 (MQ=255)
                                       acaacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaa                                                                                                       <  7:221916/149‑1 (MQ=255)
                                            ccGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCa                                                                                                                                 >  1:234739/1‑118 (MQ=255)
                                            ccGCAACATCCGCTTTCCTAACCAAATGATTGAACAAAGTAACATCGCTCTTGAGCAAAACGGGTCCGGGAATTTCTCAGCCTGGTTCATTGACGCCTGCCGCCGGAGCCTGTGCTCAGAAAAAAGCGTTTCTTCTGAAGcaaaca                                                                                                     >  3:275522/1‑146 (MQ=18)
                                            ccGCAACATCCGCTTTCCTAAACAAATGATAGAACAAAATAACATCGCTCATGAGCAAAAAGGGTCCGGGAATTTCTCAGCATGGGTCATTGAATCATCCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaaca                                                                                                     >  2:196583/1‑146 (MQ=18)
                                                aaCATCCGCTTTCCTAACCAAATGATTGAACAAATTCCCATCGCTCTTGAGCAAAAAGGGGCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaaa                                                                                              <  8:169173/149‑1 (MQ=255)
                                                aaCATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGTGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaaa                                                                                              <  7:34713/149‑1 (MQ=255)
                                                  cATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaaaga                                                                                            <  3:100232/149‑1 (MQ=255)
                                                    tCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGTAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGt                                                                                          <  1:174174/149‑1 (MQ=255)
                                                        cTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAaga                                                                                                                        >  8:352435/1‑115 (MQ=255)
                                                        cTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAaga                                                                                                                        <  7:352435/115‑1 (MQ=255)
                                                                                            tctTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTTTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGAtt                                                  >  5:125567/1‑149 (MQ=255)
                                                                                                              cGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATc                                >  8:222054/1‑149 (MQ=255)
                                                                                                                     ttCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAAtt                                         >  8:129190/1‑133 (MQ=255)
                                                                                                                                     ttGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACa                                                                                      >  4:343556/1‑72 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCCCCTCATATCACATGG  >  NC_000913/1636747‑1637035

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: