Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,164,545 G→T 6.3% D151Y (GAT→TAT)  nagZ → beta N‑acetyl‑glucosaminidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,164,5450GT 6.3% 648.3 / 7.4 199D151Y (GAT→TAT) nagZbeta N‑acetyl‑glucosaminidase
Reads supporting (aligned to +/- strand):  ref base G (113/67);  new base T (2/15);  total (115/84)
Fisher's exact test for biased strand distribution p-value = 5.02e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

GTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCC  >  NC_000913/1164396‑1164691
                                                                                                                                                      |                                                                                                                                                   
gttggttgATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTg                                                                                                                                                     <  1:1122385/150‑1 (MQ=255)
     ttgATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGt                                                                                                                                                >  1:3231895/1‑150 (MQ=255)
      tgATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTa                                                                                                                                               >  1:552896/1‑150 (MQ=255)
      tgATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTa                                                                                                                                               >  1:1685602/1‑150 (MQ=255)
      tgATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTa                                                                                                                                               >  1:1685678/1‑150 (MQ=255)
       gATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTAt                                                                                                                                              >  2:2871880/1‑150 (MQ=255)
       gATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTAt                                                                                                                                              <  2:336854/150‑1 (MQ=255)
         tGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGc                                                                                                                                           <  2:223751/150‑1 (MQ=255)
           gCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGgcatgcat                                                                                                                                         >  2:3422558/1‑150 (MQ=255)
             cAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGa                                                                                                                                       >  1:767523/1‑150 (MQ=255)
             cAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGa                                                                                                                                       >  2:1206157/1‑150 (MQ=255)
              aGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGaa                                                                                                                                      >  1:1062563/1‑150 (MQ=255)
              aGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGaa                                                                                                                                      <  2:3095205/150‑1 (MQ=255)
               gCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAg                                                                                                                                     >  1:2507925/1‑150 (MQ=255)
                 gAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGcc                                                                                                                                   >  1:517866/1‑150 (MQ=255)
                  aaaTGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCg                                                                                                                                  <  2:408466/150‑1 (MQ=255)
                     tGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGaa                                                                                                                               <  2:4013985/150‑1 (MQ=255)
                     tGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGaa                                                                                                                               <  2:3515619/150‑1 (MQ=255)
                      gATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTCTTGATGGTATGCATGAAGCCGGAAt                                                                                                                              >  1:2320407/1‑150 (MQ=255)
                      gATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAAt                                                                                                                              <  1:634632/150‑1 (MQ=255)
                      gATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAAt                                                                                                                              <  1:592805/150‑1 (MQ=255)
                      gATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAAt                                                                                                                              <  1:769281/150‑1 (MQ=255)
                      gATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAAt                                                                                                                              <  1:3318292/150‑1 (MQ=255)
                       aTCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATg                                                                                                                             >  1:3081035/1‑150 (MQ=255)
                       aTCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATg                                                                                                                             <  1:863737/150‑1 (MQ=255)
                       aTCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATg                                                                                                                             >  2:1686226/1‑150 (MQ=255)
                           cTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGaaaa                                                                                                                         >  2:2063896/1‑150 (MQ=255)
                                 atatTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTTATGGTATGCATGAAGCCGGAATGAAAACGACCg                                                                                                                   >  2:1560985/1‑150 (MQ=255)
                                          tCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACt                                                                                                          >  2:3283703/1‑150 (MQ=255)
                                          tCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACt                                                                                                          >  1:3352812/1‑150 (MQ=255)
                                          tCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACt                                                                                                          >  1:4091878/1‑150 (MQ=255)
                                            aGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTc                                                                                                        >  1:2896853/1‑150 (MQ=255)
                                            aGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTc                                                                                                        >  1:1658911/1‑150 (MQ=255)
                                              cTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTccc                                                                                                      >  2:3155082/1‑150 (MQ=255)
                                               tttGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCa                                                                                                     >  2:1762839/1‑150 (MQ=255)
                                               tttGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCa                                                                                                     >  1:1778221/1‑150 (MQ=255)
                                                tggcgcGTGTGGTGGATTTCGGGCAAATCCGGGCGGCGATGGGCGATCTTTCTTAACATGCCGATCCACAAAAAGCGCGGGAAATTGCCCGCGGGTTTTTTTTTGTTATGCATGAAGCCGCAATCAAAACGACCGGGAAACACTTCCCAg                                                                                                    <  2:3507816/148‑1 (MQ=255)
                                                  gctcgtgttCTGGATGTGGGGATTATCAGCGCGGCGTTTGGCGAGCGTTATTTTCATGCCGTTCCACAAAAAGCCCTGGCATTTGCCTGCCGGTTTTTTTATGGTTTGATTGAAGCGGGAATGAAAACGACCGGGAAACACTTCCCAGGa                                                                                                  <  2:2346668/145‑1 (MQ=255)
                                                  gcgcGTGTGCTGGTTGTGGGGCATATCAGCGGGGCTTTTGGCGTGGGTTCTTATCAGGCCGATCCACAAAAAGCCGTGTAATTTGCCCGCGGTTTTTTTTTTGGTTTGCATTAAGCCGTATTTAAAACGACCGGGAAACACTTCCCAGGa                                                                                                  <  2:3820779/150‑1 (MQ=255)
                                                     ccTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGacac                                                                                               >  2:140945/1‑150 (MQ=255)
                                                       tgtgCTGGATGTCGGGCATATCAGCGCGGCGTTTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCGGGTTTTTTGATGTTATGCTTGATGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACgg                                                                                             <  2:928213/150‑1 (MQ=255)
                                                         tgCTGTATGTCGGCCAAAACAGCGCGCATATTGCCTAGATTAATTAACATCCCAAACAAAAAAAACCCCTGCAAATTGACAGCCTTTTTATTCTTTTTATGCATGAACCCTAAATAAAAACGACCGGGAAACACTTCCCAGGACACGGTg                                                                                           <  2:2806036/150‑1 (MQ=255)
                                                         tgCGGGATGTCGGGCATAGCCGCGCGGCTATTGGCGCGCGTTCTTATCATGCCGACCCAAAAAACGCCCTGGCAATTGCCAGCCGTTTTATTGATGGTTTGCATGAAGCCGCAAGTAAAACGACCGGGAAACACTTCCCAGGACACGGTg                                                                                           <  2:3023199/150‑1 (MQ=255)
                                                           cTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCa                                                                                         >  1:2795206/1‑150 (MQ=255)
                                                             ggATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTGTTATCATGCCGATCCACAAAATGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGTAAACGACCGGGAAACACTTCGCAGGACACGGTGCAGt                                                                                       >  1:1350972/1‑150 (MQ=255)
                                                                  tCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACgg                                                                                  >  2:412483/1‑150 (MQ=255)
                                                                    gggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCa                                                                                >  2:3686602/1‑150 (MQ=255)
                                                                    gggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCa                                                                                >  2:1918048/1‑150 (MQ=255)
                                                                     ggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAg                                                                               >  1:3130662/1‑150 (MQ=255)
                                                                       cATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGAc                                                                             >  1:44689/1‑150 (MQ=255)
                                                                           tCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcac                                                                         >  2:731915/1‑150 (MQ=255)
                                                                           tCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcac                                                                         >  2:1316991/1‑150 (MQ=255)
                                                                              gcgcGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcacaca                                                                      >  2:1949178/1‑150 (MQ=255)
                                                                                gcggctTTTGGCGAGGGTTCTTATCATGCCTATCCAAAAAAAGCCGTGGAATTTGCCAGCGGGTTTTTTGTTGTTTTGATTTATGCGGGATTGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACaaa                                                                    <  2:3231876/150‑1 (MQ=255)
                                                                                    tttttGGCGTGGTTTCTTACCCTGCCCAACAACAAACAGCCGTTTAATTTGCCTGCGGTTTTTTTGTGGGTTGGCCTGAAGCCGTATTCAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGaaa                                                                <  1:1021721/147‑1 (MQ=255)
                                                                                    cgtttGGCGTGGTTTCTTACCTTCCCCATCCAAAAAAAGCGCTGGAATTGGCATGCGTTTTTTTTGTTGTTATGATTGACGCCGCAATCAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGaaa                                                                <  1:1854962/147‑1 (MQ=255)
                                                                                     gATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAac                                                               >  2:2128378/1‑150 (MQ=255)
                                                                                      aTTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAaca                                                              >  2:1279111/1‑150 (MQ=255)
                                                                                      aTTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAaca                                                              >  1:2481316/1‑150 (MQ=255)
                                                                                         ggCGAGCTTTATTATCTTGCCGATCCAAAAAAAGCCCTGGAATTTGCCTGCGGTTTTTTTTTTGTTTTGCTTGACGCGGTATTGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCg                                                           <  1:2414490/150‑1 (MQ=255)
                                                                                            gAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGc                                                        >  2:1593936/1‑150 (MQ=255)
                                                                                             aGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCg                                                       >  1:2375947/1‑150 (MQ=255)
                                                                                              gCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCTGCCGGTTTATTGTTGTTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGa                                                      <  2:692636/150‑1 (MQ=255)
                                                                                                gTTCTTCTCCGGCCTACCAACAAAAGGCCGTGGAATTGGCCAGCGGTTTTTTTGCTGTTATGCTGCCCGCCGCACTTAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATc                                                    <  1:3542005/150‑1 (MQ=255)
                                                                                                   cTTATCTTCCCTACCCACAAAAAGCCGTGCAATTTGCCCGCGGTTTTTTTTTTGTTTGTCATGACGCCGTATTCAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCAc                                                 <  2:2437086/150‑1 (MQ=255)
                                                                                                      aTCATGCCTATCCACAAAAAGCCCTGGCATTTGCCAGCCGGTTTTTTGTTGGTATGATTGAAGCGGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTc                                              <  2:3738883/150‑1 (MQ=255)
                                                                                                        cATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCca                                            >  1:1178180/1‑150 (MQ=255)
                                                                                                           gCCTATCCAAAAAAAGCCCTGGAAATTGCCTGCGGTTTTATTGTTGTTTTGAATGAAGCCGCATTGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACaa                                         <  2:2507925/150‑1 (MQ=255)
                                                                                                           gCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGGTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACaa                                         >  2:562982/1‑150 (MQ=255)
                                                                                                              gATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCg                                      >  1:3228957/1‑150 (MQ=255)
                                                                                                              gATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATGGAGGGTATGCCTGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGCGCCACGGCCACAcgcg                                      >  2:3886394/1‑150 (MQ=255)
                                                                                                                tCCACAAAAAGCCCTGGAAATTGCCAGCCGGTTTATTGTTGGTATGCATGACGCCGGAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGga                                    <  2:887582/150‑1 (MQ=255)
                                                                                                                  cacaAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGgaga                                  >  2:280315/1‑150 (MQ=255)
                                                                                                                   acaAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGAt                                 >  2:3082277/1‑150 (MQ=255)
                                                                                                                     aaacaccccTGGAATTGGCCAGCCGTTTTTTTGTTGTTATGCATGACGCCGTAATCAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTc                               <  1:2063896/144‑1 (MQ=255)
                                                                                                                     aaaacGCCCGGGAATTTGCCCGCGGTTTTATTGTTGTTTTGCATGACGCGGCACTCAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTc                               <  1:2158625/150‑1 (MQ=255)
                                                                                                                      aaaaGCCCTGGCATTTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcg                              <  2:3349656/150‑1 (MQ=255)
                                                                                                                       aacgCCCTGGATATTGCCTGCGGTTTTTTTGATGGTCTGACTGACGCCGTTATCAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcgc                             <  1:460542/147‑1 (MQ=255)
                                                                                                                       aaaGCCGTGGAATTTGCCAGCCGTTTTATTGTTGGTATGCATGACGCCGGATTTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcgc                             <  1:2574122/150‑1 (MQ=255)
                                                                                                                        aaGCCCTGTCATTTGCATGCGGTTTTTTTGTTGTTCTGAATTAAGCGGGAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcgcg                            <  2:1495227/150‑1 (MQ=255)
                                                                                                                             cTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTaaa                       >  1:2092441/1‑150 (MQ=255)
                                                                                                                             cTGGAATTTGCCAGCCGTTTTTTTTTTGTTTTGCTTGAAGCCGTAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTaaa                       <  2:552896/150‑1 (MQ=255)
                                                                                                                              tGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAg                      >  2:2444794/1‑150 (MQ=255)
                                                                                                                              tGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAg                      >  2:1303764/1‑150 (MQ=255)
                                                                                                                               ggaatttGCCAGCCGTTTTATTGTTGGTTTGCATGACGCGGGATTTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGa                     <  2:1944576/145‑1 (MQ=255)
                                                                                                                               ggCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGa                     <  2:3391884/150‑1 (MQ=255)
                                                                                                                                gaatttGCCAGCGGTTTTTTTTATGGTATGCATGACGCGGTAATCAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGat                    <  1:2531830/146‑1 (MQ=255)
                                                                                                                                  attggCCAGCGGTTTTTTTGATGTTCTGCTTTAAGCCGTATTGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGatat                  <  2:84756/146‑1 (MQ=255)
                                                                                                                                  aaTTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGatat                  >  1:1421930/1‑150 (MQ=255)
                                                                                                                                  aaTTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGatat                  >  1:269946/1‑150 (MQ=255)
                                                                                                                                    ttGCCTGCGGTTTTTTTGTTGGTATTCATTAAGCCGGAATCAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGt                <  1:1868221/150‑1 (MQ=255)
                                                                                                                                     tGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCGAAAGATATGTc               >  1:2691511/1‑150 (MQ=255)
                                                                                                                                         ccgCCGTTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGt            <  2:2647164/148‑1 (MQ=255)
                                                                                                                                         aGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTc           >  2:3315569/1‑150 (MQ=255)
                                                                                                                                           gCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCt          >  2:2396188/1‑150 (MQ=255)
                                                                                                                                           gCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCt          >  2:1958883/1‑150 (MQ=255)
                                                                                                                                            ccGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCtt         >  1:1297282/1‑150 (MQ=255)
                                                                                                                                             cGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTc        >  2:2978731/1‑150 (MQ=255)
                                                                                                                                                 tttttgttgTTATGCATGATGCGGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGtt    <  2:517866/147‑1 (MQ=255)
                                                                                                                                                   aTTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTcc  >  2:1093835/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCC  >  NC_000913/1164396‑1164691

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