Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 408,567 A→G 7.7% intergenic (+106/‑42) yaiE → / → ykiA UPF0345 family protein/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913408,5670AG 7.7% 141.1 / 2.5 52intergenic (+106/‑42)yaiE/ykiAUPF0345 family protein/pseudogene
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GAAGGCGAATACACCTTCAGCACCGCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTT  >  NC_000913/408270‑408866
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                           
gAAGGCGAATACACCTTCAGCACCGCTGAGCCGGAAGAGATGACGGTAATCAGTGGCTCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCAGGTTCGGTGTTTAATGTTCCCGGTCACAGTGCGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGACTTCCCCTTTAACGGGAGGGCATTTATCTGTAATATCCTTTCTTTAGCACATAATAAAATTTCCATAGCTCCGAATTTTTACATTTCCAATATATTCaaaa                                                                                                                                                                                                                                                                                                          >  1:474430/1‑301 (MQ=255)
              cTTCAGCACCGTTGAGTCGGAAGAGATGACGGTAATTAGTGGCGCGCTGATGTTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGa                                                                                                                                                                                                                                                                                            <  1:511713/301‑1 (MQ=255)
                  tgtttttccgttctggAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGAtata                                                                                                                                                                                                                                                                                        <  2:20534/287‑1 (MQ=255)
                  aGCACCGCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGAtata                                                                                                                                                                                                                                                                                        >  1:142724/1‑301 (MQ=255)
                   gCACCGCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGACGCCGGTTCGGTTTTTACTGTTCCCGGTCACAGATAGTTTCATCTGCAAGTTACCAAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGTGAGGCCATTTTTCTGCAATATACTTTCTTTAGCCCACAATAATATTTCCTTTGCTTCGCTTTTTTCAATTTCAGATATATTCGTAATTTATCAATGTTGATATaa                                                                                                                                                                                                                                                                                       >  1:633106/1‑301 (MQ=255)
                   gCACCGCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATaa                                                                                                                                                                                                                                                                                       >  2:169226/1‑301 (MQ=255)
                   gCACCGCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGACTTCCCCTTGAACGGGCGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGAAATATTCGTAATTTATAAAGTTTGAAATaa                                                                                                                                                                                                                                                                                       >  1:633086/1‑301 (MQ=255)
                     aCCTCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATTTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAat                                                                                                                                                                                                                                                                                     <  2:142849/301‑1 (MQ=255)
                        gCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGCAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGAAATATTCATAATTTATCAATGTTAATATAAatatt                                                                                                                                                                                                                                                                                  >  1:467249/1‑300 (MQ=255)
                          tGAGCCGGAAGAGATGACGTTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAg                                                                                                                                                                                                                                                                                <  2:89965/301‑1 (MQ=255)
                                ggAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATaaaatcggaagaccg                                                                                                                                                                                                                                                                          >  1:20483/1‑289 (MQ=255)
                                ggAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGa                                                                                                                                                                                                                                                                          >  2:484952/1‑301 (MQ=255)
                                                          tgtgtgctcttccgatctcTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGc                                                                                                                                                                                                                                                <  2:279095/283‑1 (MQ=255)
                                                                          gCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTCTCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATACTACTATTTCCTTTGCTGCGATTTTTTCAATTTCCGATCTCTTCATATTTTATCTAGGTTGATATAAATATCAGTTAAGATCCCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATaaa                                                                                                                                                                                                                                >  2:185658/1‑301 (MQ=255)
                                                                            cTGACGTGACCGACTGGCAGGTGTATTAAGCCGGTTCGGTGTTTAATTTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATc                                                                                                                                                                                                                              <  1:484351/301‑1 (MQ=255)
                                                                            cTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGcagatcggaaaaacgtcgg                                                                                                                                                                                                                              >  1:278792/1‑283 (MQ=255)
                                                                             tGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCa                                                                                                                                                                                                                             >  2:299399/1‑301 (MQ=255)
                                                                              gtcgcgTCCGACTGGCAGGTGTATTAAGCCGGTTCGGTGTTTAATGTTCCCGGGCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAg                                                                                                                                                                                                                            <  1:426879/299‑1 (MQ=255)
                                                                              gACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGAAAAAAGATAAATCaa                                                                                                                                                                                                                            >  2:427443/1‑300 (MQ=255)
                                                                                 gcgACCGACTGGCAGGTGTATGTAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTGTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAgatg                                                                                                                                                                                                                         <  2:475031/301‑1 (MQ=255)
                                                                                        acTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAAtg                                                                                                                                                                                                                  >  2:294982/1‑301 (MQ=255)
                                                                                           ggCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGACATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCACTTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAAATGATGAAtggtg                                                                                                                                                                                                               >  1:305889/1‑301 (MQ=255)
                                                                                            gctggTGTTTGAAGCCGGTTCGGTTGTTAATGTTCCCGGTCACATTGATTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCGGTGCCGCTATCTGTAATTCCTCGCCTGCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAAtggtgg                                                                                                                                                                                                              <  1:285047/298‑1 (MQ=255)
                                                                                                gtgtATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTg                                                                                                                                                                                                          <  2:373771/301‑1 (MQ=255)
                                                                                                        agtcGGTTCGTTGTTTAATGTTCCCGGTCACAGTGAGTTTAATCTGCAAGTTGCCGAACCCAACTCTTATATGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATACTTTCTTTAGCCCATAATAATATTTCCTTTGTTGCGATTTTTTCAATTTCCGATATATTCGTAATGTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATtgctg                                                                                                                                                                                                  <  2:633983/298‑1 (MQ=255)
                                                                                                        aGCCTGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAAGATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCGTAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATtgctg                                                                                                                                                                                                  <  2:633963/301‑1 (MQ=255)
                                                                                                                         aTGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCa                                                                                                                                                                                 >  2:66352/1‑301 (MQ=255)
                                                                                                                           gTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCAtt                                                                                                                                                                               >  2:77614/1‑301 (MQ=255)
                                                                                                                                    aacaGTGAGTTTCATCTGCAAGTTGCCGAACCCACATCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAAGGATCATTTACGCCaaa                                                                                                                                                                      <  1:185416/300‑1 (MQ=255)
                                                                                                                                     tcaGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCGATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATAGATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAAc                                                                                                                                                                     <  2:559638/300‑1 (MQ=255)
                                                                                                                                        gTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGc                                                                                                                                                                  >  2:546860/1‑301 (MQ=255)
                                                                                                                                              ttCGTCTTCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTgcg                                                                                                                                                            <  1:66262/301‑1 (MQ=255)
                                                                                                                                                aaTCTGCAAGTTGCCGAAGCCTCCTCTTATCTGTGCGGCTATCTGGAATTCCGCGCCTTCCCCTTGAACGTGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAAGAATATTTCCTTTGCTGCGATTTTTTAAATTTCCGATATATTCATAATTTAGCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCa                                                                                                                                                          <  1:19422/300‑1 (MQ=255)
                                                                                                                                                    tGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTACACAACGc                                                                                                                                                      >  1:283190/1‑301 (MQ=255)
                                                                                                                                                    tGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGAATTGCGCAACGc                                                                                                                                                      >  2:25500/1‑301 (MQ=255)
                                                                                                                                                        tattcactctttccctacacgacgctctttcTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATAAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACAt                                                                                                                                                  <  1:432589/271‑1 (MQ=255)
                                                                                                                                                            gCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATAccc                                                                                                                                              >  2:173106/1‑301 (MQ=255)
                                                                                                                                                              cGAACCCACCTCTTTTCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCg                                                                                                                                            <  1:25449/301‑1 (MQ=255)
                                                                                                                                                                                   ccGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCAAATACCCAGGCGTTGTTGTTAACGACAttt                                                                                                                       >  1:143901/1‑300 (MQ=255)
                                                                                                                                                                                                 tCCTCGCCTTTCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGa                                                                                                         <  1:546187/301‑1 (MQ=255)
                                                                                                                                                                                                  tcTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGAt                                                                                                        <  2:62062/300‑1 (MQ=255)
                                                                                                                                                                                                     tacttgagttcagacgtgtgctcttccgatCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCgg                                                                                                     <  2:386656/272‑1 (MQ=255)
                                                                                                                                                                                                        cTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACCGATCCTGAGGCTCAGTca                                                                                                  >  2:262607/1‑300 (MQ=255)
                                                                                                                                                                                                          tCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACCGATCCTGAGGCTCAGTcaat                                                                                                <  1:172909/301‑4 (MQ=255)
                                                                                                                                                                                                               ttGAACTGGAGGGCATTTTTCTGAATTTTCCTTTCTTTAGCCCATATTAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACCGATCCTGAGGCTCAGTcaatgtcgt                                                                                           <  1:262339/301‑9 (MQ=255)
                                                                                                                                                                                                                 gAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTACCGACATTGACTGAGCCTCAGGCTCGGTCCCTGGCTAtt                                                                                         >  1:260067/1‑301 (MQ=255)
                                                                                                                                                                                                                                              tttctttAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGaa                                                            <  2:467813/301‑1 (MQ=255)
                                                                                                                                                                                                                                                        ccATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCcagatcggaagagcacacgtctgaactccagt                                                  >  2:135770/1‑271 (MQ=255)
                                                                                                                                                                                                                                                          ataataataTTTCCTTTGTTGAGATTTTTGCAATTTCTGATATATTCATAAGTTATCAAGGTGGATATAAATATCAGTGAAGAGCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATGGCTGCAAAGTTGCAATGATCATTTACGCCAAACGGAATTGCGCAACGCAAATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGtt                                                <  1:348973/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ttGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatc                                 <  2:144027/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  ttttttCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCtt                        >  2:611144/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     ttttAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTAtt                     <  2:260335/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     tttCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTAtt                     <  1:77524/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                        taTAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCtt  <  2:306293/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                           
GAAGGCGAATACACCTTCAGCACCGCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGATATATTCATAATTTATCAAGGTTGATATAAATATCAGTGAAGATCTCCAGATATTGTTGCGGAACTGGCTACGATAAAAGATAAATCAGATGATGAATGGTGGCGTGCATTGCTGCAAAGTCGCAATGATCATTTACGCCAAACGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTT  >  NC_000913/408270‑408866

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: