Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,861 T→G 17.1% G27G (GGT→GGG)  rpsU → 30S ribosomal subunit protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,210,8610TG17.1% 152.0 / 4.0 74G27G (GGT→GGG) rpsU30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base T (21/33);  new base G (16/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 2.68e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTT  >  NC_000913/3210637‑3211077
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                          
ctcctcACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTACTCAAAGGCGAGAGGCCCATTCCGGTAATTAAAGGACGTGAAAAC‑AAGCCGTTTGAAGAAGCCCTGCGCTCATTAAACCGTTCCTCGCAAAAAACAGgtggt                                                                                                                                                                                                                         >  2:626745/1‑226 (MQ=255)
   ctcACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAACAGCAGGTGTTATg                                                                                                                                                                                                                      >  1:404101/1‑228 (MQ=255)
   ctcACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGGGCTTTACAAACCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGAAAATGCCGGAAATTAAAGTACGAGAAAAC‑GAGCCGTCCTACGAAGCTCTGCGCGCCTTAAAACGTCACTGCGAAAAAGCAGAGGTTccg                                                                                                                                                                                                                      >  2:376538/1‑226 (MQ=255)
     cACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTCCGAAAATGCAGGGGTTCTgcc                                                                                                                                                                                                                    >  1:98477/1‑226 (MQ=255)
      aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACAGGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTGCCTTCAACCGTTCTTGCGGAAAAGCAGTGGTTCTGCCg                                                                                                                                                                                                                   >  1:47129/1‑228 (MQ=255)
           tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                                                                                                                                                                                              >  1:450293/1‑228 (MQ=255)
           tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑AGCCCGTTCGACGAACCCCTGCGACCCTCCAACGCTTCCTCCGAAAAAGCCGGGGTTTTGCCGGAAGt                                                                                                                                                                                                              >  2:373799/1‑228 (MQ=255)
           tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGAACGTGAAACC‑GAGCCGTTCGACGAAGCCCTGCGTCCGTTCAACCGTTCCTGCGAAAAAGCAGGAGTGCTGGCGGAAGt                                                                                                                                                                                                              >  1:570816/1‑228 (MQ=255)
           tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGACGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGGGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                                                                                                                                                                                              >  1:476799/1‑228 (MQ=255)
             taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGCTGCCAGCGGCAAACCGAATTAATCCAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCCCTTCAAGCGTTCCTTCGAAAAACCAGGGGTTCTGGCGGACGttg                                                                                                                                                                                                            >  2:386165/1‑227 (MQ=255)
             taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGACGGTAATTAAAGCACGTGAAAAC‑GAGCCGTTCGACGAAGCCATGCGTCGCTTCAACCGTTCCAGCAAACAAGCAGGGGTTCTggagagagata                                                                                                                                                                                                            >  2:109317/1‑219 (MQ=255)
                  gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGTCGGAAGTtcgacgt                                                                                                                                                                                                       >  1:624723/1‑228 (MQ=255)
                  gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTtcgtcgt                                                                                                                                                                                                       >  1:124082/1‑228 (MQ=255)
                      cTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATACCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTcggtcgtcggg                                                                                                                                                                                                   >  1:91877/1‑226 (MQ=255)
                       tttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGAACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCCGCGAAAAAGCAGGGGTTCTGGCGAAACTtagtcgccgggg                                                                                                                                                                                                  >  2:98575/1‑225 (MQ=255)
                        ttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGGGAg                                                                                                                                                                                                 >  1:645906/1‑228 (MQ=255)
                        ttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTCCGAAAAAGCAGGCGTTCTGGCGGAAGTTCGCCGTCGTGAg                                                                                                                                                                                                 >  1:47224/1‑228 (MQ=255)
                          cAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATCCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGtt                                                                                                                                                                                               <  1:319847/228‑1 (MQ=255)
                          cAAAAAAGGCCGCGGGGCTTTACAAAGAAGCAGCAATGGCAGTAAAAGTCCGCACAATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGGGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGtt                                                                                                                                                                                               <  1:385699/228‑1 (MQ=255)
                                           cTTTACAAGGCGGAACAAATTGCAGAAAAATCCCGCACCATTGTGAAATAAGCTGCGGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGCGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa                                                                                                                                                                              <  1:329810/228‑1 (MQ=255)
                                             ttAAAAAGCAGTAGCAATTGCAGTAAAATTCCGCACGATTTGGAAATAAGCTGGCGTTGATGCAAGCGGCAAACCGAATGAAGCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTAcc                                                                                                                                                                            <  1:417562/228‑1 (MQ=255)
                                                   agcagcagcAATTGCAGTAAATTTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCCAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGACGTACGCCGTCGGGAGTCCTATGCAAAACAGACTACCTAACgc                                                                                                                                                                      >  1:212638/1‑227 (MQ=255)
                                                   agcagcagcAATTGCAGTAAAATTCCGCACCATTTTGAAATAATCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAACC‑GATCCGTTCGACGTAGCTCTGCGTCGCTTCAAGGGTTCATGCGAAAAAGGAGGGGTTCTGTCGGAAGTTCGTCGCCGTGAGTTCTATGAGAAAACGAATACCGAACGt                                                                                                                                                                      >  1:52834/1‑228 (MQ=255)
                                                   agcagcagcAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGAGGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGGGGAAGTTCGTCGGCGTGAGTTCTATGAAAAACCGCCAACCAAACga                                                                                                                                                                      >  2:544015/1‑227 (MQ=255)
                                                   agcagcagcAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCCCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGGAGTTCGTCCTGGTGGGTTCTTTGAAAAACCGACTACCGAACga                                                                                                                                                                      >  2:329810/1‑227 (MQ=255)
                                                     cagcagcaATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAACCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGGGGGTTTCATGAAAAAACAACAACCCAACGaaa                                                                                                                                                                    >  2:634413/1‑228 (MQ=255)
                                                      ggaagaaaTTGCAGTAAAATTCCGCCCCATTTTGAAATAAGCTGGCGTTGAGGCCAGCGGAAAACCGAATTAATCAAAGGTGAGAGGCCCATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAg                                                                                                                                                                   <  1:649726/222‑1 (MQ=255)
                                                       gaataaaTTGCAGTAAAATTCCGAACAATTTGGAAAGAAGCTGGGGTGGAGCCCCGGGGAAACCCGAAATAATCAAAGGGGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAg                                                                                                                                                                   <  2:557690/222‑1 (MQ=255)
                                                        cagcaATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGGAGCTCTGCGTCGCTTCAAACGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTGGTCGTGGGTTCTATGAAAAAACAACTCCCGAACGaaaac                                                                                                                                                                  >  2:544068/1‑225 (MQ=255)
                                                            aaTTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGGGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GTGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCt                                                                                                                                                             <  1:304124/228‑1 (MQ=255)
                                                            aaTTACAGTAAAATTCCGCACCATTTTGAATTAAGGTGCCGTTGATGCCATCGGCAAGCCGAATGAACCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAACC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCt                                                                                                                                                             <  2:404101/228‑1 (MQ=255)
                                                              tgtaaGTAAAGTTCCACAACATTTTGAAATAATGTGGCTTTGATGCCAGGGGAAAACCGAAGTAATCAAAGGTGAGAGGCACATGCCGTTAATTAAATGACGTGAAGAC‑GAGCCGGTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTaa                                                                                                                                                           <  2:424474/224‑1 (MQ=255)
                                                               tGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTTCTAAAAAGCAGGGGTTCTGGCGGAAGTCCGTCGTCGTGGGTTCTAGGAAAAACCGACTACCGCACGTAAGCGCGCTaaa                                                                                                                                                          >  2:232334/1‑228 (MQ=255)
                                                               tGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTaaa                                                                                                                                                          >  1:112394/1‑228 (MQ=255)
                                                                              cACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAGTCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAACGCGCCAAAGCTTTTGCAGTGaaa                                                                                                                                           >  2:39360/1‑228 (MQ=255)
                                                                                   ttttGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGGCTCCCGAACGTAAGGGCGCTAAAGCTCCGGCAGGGAAACGta                                                                                                                                       >  1:567969/1‑226 (MQ=255)
                                                                                           tAAGCTGGCGTGGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCga                                                                                                                                 <  1:335439/225‑1 (MQ=255)
                                                                                              gCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGt                                                                                                                                                                                                >  1:425520/1‑159 (MQ=255)
                                                                                               tgggggTTGAGGCCAGCGGAAACCCGAATAAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGtt                                                                                                                                                                                               <  2:425520/154‑1 (MQ=255)
                                                                                                                 gCAAACCGAATTACGCAAGGGTAAGAGGCAAATGCCGGTAATTAAAGTACGTGAAAGC‑GAGCCGTTCGACGTAGCTCTTCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAGACTGGCTCGC‑GAAAacgcacg                                                                                                       <  1:387139/228‑1 (MQ=255)
                                                                                                                       ccaaTTAAACAAAGGGGAGAGGCAAATGCCGGTAATTAAAAGTCGTGAAAAC‑GAGCCGTTCGACGTAGCTCGGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCa                                                                                                  <  2:645906/225‑1 (MQ=255)
                                                                                                                       cGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCa                                                                                                  <  1:544068/227‑1 (MQ=255)
                                                                                                                       agAATCATTCAAAGGGGAGGGGACCATGCCGGTAATTAAAGTACGTGAAAAC‑GAACCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCa                                                                                                  <  2:450293/226‑1 (MQ=255)
                                                                                                                       agAAATAATCAAAGGTGAGAGGCCAATGCCGGTAATTAAAGGTCGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCa                                                                                                  <  1:517060/226‑1 (MQ=255)
                                                                                                                              aacAAAGGTAGGAGCAACATGCCGGTAATTAAAGTGCGGTAAAAC‑GAACCGTTCGACGTAGCTCTGCGTCGCTTCAAGCTTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTg                                                                                          <  1:109317/226‑1 (MQ=255)
                                                                                                                                    ggggtgagGCACATTCAGGTAAATAAAGTATGTGAAAAC‑GGGCCGGTGGCCGTAGCTCTGGGTCGCTGCAAGTGTTCCTGCGAAAAATCAGGTGTTCTGGCGGAAGTTCGGCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTaa                                                                                    <  2:567241/223‑1 (MQ=255)
                                                                                                                                               caTGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCCCGCGCAGAAACTGGC‑CGCGGAAAACCCACGCCCCACCACTCTGCACCAATTcacccaca                                                                          >  2:296940/1‑223 (MQ=255)
                                                                                                                                               caTGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACAGGCCCGC‑AAAAAGGCCACCCCCACCACGCTGTACCAATTcccacgaag                                                                         >  2:501338/1‑222 (MQ=255)
                                                                                                                                                  gcgggTAATTAAAGGACGTGAAACC‑GGGCCGTGCGACGTAGCTCTGCCTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCg                                                                       <  2:567969/224‑1 (MQ=255)
                                                                                                                                                    cGGTAATTAAAGGTCGGGAAAAC‑GGGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCTTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGtt                                                                     <  1:477376/227‑1 (MQ=255)
                                                                                                                                                    cGGTAATTAAAGGCAGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGtt                                                                     <  1:152160/227‑1 (MQ=255)
                                                                                                                                                    cGGTAATTAAAGGACGTGAAAGC‑GAGCCGGTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGtt                                                                     <  1:373799/227‑1 (MQ=255)
                                                                                                                                                    aggaaaTTAAAGGTGGTGCAAAC‑GGGCCTTTCGTCGTTGCTCGGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCCGTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGtt                                                                     <  2:475259/223‑1 (MQ=255)
                                                                                                                                                      gTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCCCTAAAGATTCTTCAGGGAAACGTCACCGGAAGAAACTGGCACGC‑GAAAACGCCCCCCCCCCTCGTCTGAACTAATCCCCCgagacagtggat                                                                  >  2:129970/1‑220 (MQ=255)
                                                                                                                                                        aagtAAAGGACGTTCAAAC‑GGGCCGTGCGATTTAGCCTGGCGTCGATTCAAGCGTTCCGGCGAAAAAGCTGGGGTTCTGGCGGAAGGTCGACGCCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCTGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑ctc                                                                <  2:476799/224‑1 (MQ=255)
                                                                                                                                                        aaTTAAAGTACGTTAAAAC‑GAGCCGTTCGACGTAGCTCTGCGGCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTAAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGtt                                                                     <  1:232334/223‑1 (MQ=255)
                                                                                                                                                        aaTTAAAGT‑CACGAAAACGGAGCCGTTCGACGTAACTCTGCGTCGCGTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑ctc                                                                <  1:98575/227‑1 (MQ=255)
                                                                                                                                                         aTTAAAGTTCGTTAAAAC‑GGGCCGTTCGACGTAGCCCTGCGGCGCTTCAAGCGTTCCGGCGAAAAAGCAGGTGTTCTGGCGGAAGTTGGTCGTCGTGAGTTCTATGAAAAACCGACCACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑ctct                                                               <  2:624723/227‑1 (MQ=255)
                                                                                                                                                             aaagtggTGAAAAC‑GAGCCGGTCGGATTTGCTTTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGGTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGAt                                                          <  2:498202/222‑1 (MQ=255)
                                                                                                                                                                   gTGAAAAC‑GTGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGtt                                                                     <  1:30824/212‑1 (MQ=255)
                                                                                                                                                                    tGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGACACCGTTG‑CTCTCCCCTCAGAcca                                                   >  1:458226/1‑227 (MQ=255)
                                                                                                                                                                             ggCCGTTCGACGTAGCTCTGCGTCGGTTCAAGCGTTCCTGGGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGagttgt                                             <  1:162849/225‑1 (MQ=255)
                                                                                                                                                                             ggCCGTTCGACGTAGCTCTGCGTCGCTTAAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTACGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGagttgt                                             <  1:544015/225‑1 (MQ=255)
                                                                                                                                                                                                 ggTCGTTTCAAGGGTTCCTGCGAAAAAGCAGGTGTTCCGGGGGAATTTCGTCGTCGTGAGTTCGATGAAAAACCGACTACCGAACGTAAGCGCGGTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                       <  2:91877/227‑1 (MQ=255)
                                                                                                                                                                                                 cGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGTAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTCAGCGCGCTAAAGCTTCTGCAGTGAAACGCCACGGTACGAACCTGGCTCCC‑CAAAACGACCGCCGCACTCGACCGTACTAATTCCCCGACAGCGCTGACTCT‑CGACCAACCCGAGTTCTAGTATTAAGCCGAGGCTTCCGa                       >  2:432907/1‑228 (MQ=255)
                                                                                                                                                                                                 cGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGCAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCt                                                                                                                                                       >  1:187347/1‑102 (MQ=255)
                                                                                                                                                                                                 cGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGACAGCGTTG‑CTCTCCGATCCAACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                       >  1:410863/1‑228 (MQ=255)
                                                                                                                                                                                                  gTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCtt                                                                                                                                                      <  2:187347/102‑1 (MQ=255)
                                                                                                                                                                                                      cTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAACCTGGCTCGC‑GAAAACGCACGCCGCCCTCGTCTGTACTAATTCCCCCAGAGCGTTG‑CTCTCCGATCAGACCGAGTTTTAGTTGTAAGTCCGt                               >  1:228521/1‑215 (MQ=255)
                                                                                                                                                                                                      cTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGGGTTGTAGTTGTAAGGCCGTGCTTCCGCAAGGa                  >  1:142558/1‑228 (MQ=255)
                                                                                                                                                                                                           aGCGTTCCTGCGAAAAAGCAGGGTTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATgcg             <  1:39360/228‑1 (MQ=255)
                                                                                                                                                                                                                 ccTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGGGTTGGAGTTGGAAGGCCGTGCTTCCCAACGGAAAGCCCGGCtt       >  2:39036/1‑228 (MQ=255)
                                                                                                                                                                                                                     cGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGCCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGCGTTTTAGTTGTAAGCCCCGGCTTCCGCAAGGAATGGCCGGCTTAttt   >  1:628674/1‑228 (MQ=255)
                                                                                                                                                                                                                     aggaaaaGCAGGTGTTCGGGGGGAAGTTCGGCGTCGGTATTTCTAAGGAAAACCGACCACAGAACGTAAGCGCGTGAAAGTTTCTGCAGTGAAACGTCACGCGAAGAAGCTGGCTCGC‑GAAAACCCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTAttt   <  2:419601/225‑1 (MQ=255)
                                                                                                                                                                                                                      gAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCCC‑GAAAACGCACGCCGCACTCGCCTGTCCTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGAGTTGTAGTTGTCAGGCCGGGCTTCCGAAAGGAATGCGCGGCTATtttt  >  2:168459/1‑228 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                          
CTCCTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAAC‑GAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGC‑GAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTG‑CTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTT  >  NC_000913/3210637‑3211077

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