Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,721,136 C→T L38F (CTT→TTT)  ydhI → DUF1656 domain‑containing protein YdhI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,721,1360CT94.0% 195.6 / ‑4.4 67L38F (CTT→TTT) ydhIDUF1656 domain‑containing protein YdhI
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (35/28);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GCTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGC  >  NC_000913/1720891‑1721371
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
aTTGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGt                                                                                                                                                                                                                                          <  1:642361/250‑1 (MQ=255)
  tGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑Gt                                                                                                                                                                                                                                      >  2:297149/1‑251 (MQ=255)
  tGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGGTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑g                                                                                                                                                                                                                                       <  1:720980/250‑1 (MQ=255)
     tAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGCAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAAAAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTac                                                                                                                                                                                                                                    >  1:744487/1‑250 (MQ=255)
                     ctgcATTTGGTTAGCGTAAGACATATTTTTAAGAAATTGCGCGTAAATTTTTCGAATTTAAAGAATATTTATTAATCTAACGCAATATATTTGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTAAATGCAACAGGATTGTCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTAAAAGCATTCGCGTTTGGATTCGTAATATG‑GTTTGTG‑TTACACTGCCTGATTTGTggt                                                                                                                                                                                                                   >  2:382710/1‑250 (MQ=255)
                       gcATTTGGTTACCGTCAGACATATTTTCCAGAAATTGCGCGTAAATTTTTCGCCCTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑ACGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGGCCACTTTCCTCCGTTTTTCAACGCATTCGCGTGTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGc                                                                                                                                                                                                                   <  1:33996/249‑1 (MQ=255)
                         aTTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑TTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGa                                                                                                                                                                                                               <  1:776959/251‑1 (MQ=255)
                            tGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATc                                                                                                                                                                                                             <  1:667533/250‑1 (MQ=255)
                                  ccGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCGAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTTCTTTCCTCCGTTTTTCAAAGCTTTCGCGTTTGGATTCGTCAGCTG‑GTTTGTC‑GTACACCGCCTGCGTCGTGGCTGGATCTACGcc                                                                                                                                                                                                       <  2:358635/250‑1 (MQ=255)
                                    gTCAGACAGATTTGTCATAAATTGCGTGTAAAGTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGTTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGt                                                                                                                                                                                                    <  2:232324/251‑1 (MQ=255)
                                      tagccATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGAAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTGATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGAGTTCGGCGAGTCCGTCTAATGTCCTCCGTTTTTCAAAGCATTCTCGATTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTAGTGGCTGAATCTACGcagg                                                                                                                                                                                                     >  2:283986/5‑245 (MQ=255)
                                       aGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGAc                                                                                                                                                                                                 >  1:863212/1‑251 (MQ=255)
                                           atatTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTT‑CCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACTCCGGTGACATCTg                                                                                                                                                                                            >  1:264840/1‑251 (MQ=255)
                                                       taTTGCGCGAAAATTTTTCGAACATAAAGAATATTTATTAATATAACAAAATATATTCTGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCACAATGCAACAGGCATGCACTTACAAGATATGGTGTTCGGTGCGTACGACTACTTTCTTCCGTTTTTCAAAGCATTCGCGTTTGGATGCGTCATCTG‑GTTTGTC‑GATCACCGCCTGCTTCGTGTCGGGATATACACCGGTTACATGTGGCATCCATTGt                                                                                                                                                                                 >  2:518866/2‑251 (MQ=255)
                                                       aaTTGCGCGGAAATTTTGCGCACTTAAATAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATTTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGCTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGCC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGt                                                                                                                                                                                 <  1:818631/251‑1 (MQ=255)
                                                                 aaaTTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGAttt                                                                                                                                                                       <  2:332657/251‑1 (MQ=255)
                                                                   aTTTTTCGCACTTAAAGACTATTTATTAATCTACCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGTTCTGGTGTTCGGTTCGTCCGTCTACTTTCCTCCGTTTTTCATATCATTCGCGATTGGATTCGTCCTCTG‑TTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTATGTCGGTGTCATCTGGCATCCCTTGTTAATGTAttttt                                                                                                                                                                     >  2:279333/1‑249 (MQ=255)
                                                                      tttCGCACTTAAAGAATATTTATTAACCGAACTAAATATATTCGGTCGTAAAGGAAACTACTTTGTG‑AAGTTTATGCCCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GCTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCg                                                                                                                                                                   <  2:413201/250‑1 (MQ=255)
                                                                         cgcaatgaaaGAATATTAATTAATCATACGCTATAGACTCGGTCGGAAAGGAATCTTCTGTGTG‑AAGTTTACTCTCAACGAATCAGGATTCCCCTTACAAGATCGGTTGTGCTGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTT‑GTTTGTC‑TAACACCGCCTGCTTCGTGGCTTGATCTACGCCGGTGACATCTGGCATCCCTTTTTAAGGGATTTATCGCTGt                                                                                                                                                               <  2:252453/244‑1 (MQ=255)
                                                                          gCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGAATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGtt                                                                                                                                                              <  1:775938/251‑1 (MQ=255)
                                                                                  ataattttgaatattttaacgacatatatTCGGTCGTAAAGTACTCTACTTTTTG‑AAGTTTCTGCTCAATGCCACAGGATTGCCCTTACAAACTCTGGTGTTCGGTTCGTCCTGCTACTTACCTCCGTTTTTCAAAGCATTCGCTTTTTGATTCGCCAT‑TGCGTTTTTC‑CGACTCCGCCTGCTTCGTCGCTGGATCTACGCCGTTGTCATCTGGCATCCCTTGTTAATGGATTTATCGCTGTGTGCGAtct                                                                                                                                                      <  2:58902/228‑3 (MQ=255)
                                                                                   gAATATTTATTAATCTAACACAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTAAATGCAACAGGATTGCCCTTACAAGTTCTGGTGTTTGTTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGGGACATATGGCATCCCTTGTTAATGGTTTTTTCCCTGTTTgagctttg                                                                                                                                                     >  2:803940/1‑244 (MQ=255)
                                                                                      tatttattAATATAACGCAATATATTCGGACGTAAAGGAATCTACTTTATG‑AAGTTTATGCTCAATGCAAAAGGATTGCCCTTACAAGATCTGGTGTTCGGTGAGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGCTTCGTCATCTG‑TTTTGTC‑GTACACCGCCTGCTTCGTAGCTGACTCTACGCCGGTGCCATCTGGCATCACTTGTTAATGGATTTATCGCTGTTTGCGagtatcgt                                                                                                                                                  >  2:719536/1‑244 (MQ=255)
                                                                                       tattcttaATTTAACGTAATATATTCGGTCGTAAAGGCATCTCCTTT‑TGTAATGTTATGTCCAATGCAACAGGATTGCCCTTACAAGATATGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GCTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATGTATCGCTGTTTGCGATTTGCGt                                                                                                                                                  <  2:393780/245‑1 (MQ=255)
                                                                                         ttattaATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGttt                                                                                                                                                <  1:11621/250‑1 (MQ=255)
                                                                                              aaTCTAACGCAAGATAGTCGGTCGTAAAGGAATCT‑TTTTGTTTAATTTTATGCTCAGTGCAGTCGGATTGCCCTTACAAGATTTGTTGTTTGGTTCGTCCGTTTGCTTTCCTCCGTTTTACAAAGCATCCTCGTTTGTATTCGTCCTCTG‑GTGTGTC‑GTACCCCGCCTGCTTCGTGGCTGGATCTGCGCCGGTGACATCTGGCATCCCTTGTTAATGGAATTATCGCTGTTTGCGATTtgcctttgcctt                                                                                                                                           <  2:767606/250‑1 (MQ=255)
                                                                                                 cTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGAGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGctc                                                                                                                                       >  1:428671/1‑251 (MQ=255)
                                                                                                    aCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTCTGCTCAATGAAACAGGATTGCCCTTCCAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGCTCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGCTTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTgg                                                                                                                                    >  2:852276/1‑251 (MQ=255)
                                                                                                         atatatTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGTTGACATCTGGCATCCCTTGTTAAt                                                                                                                                                                             >  2:22458/1‑205 (MQ=255)
                                                                                                         atatatTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAAt                                                                                                                                                                             <  1:22458/205‑1 (MQ=255)
                                                                                                             atTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCACTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGAGCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTg                                                                                                                            >  2:775791/1‑250 (MQ=255)
                                                                                                                cGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTg                                                                                                                        >  2:734710/1‑251 (MQ=255)
                                                                                                                      tAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGt                                                                                                                  >  1:89195/1‑251 (MQ=255)
                                                                                                                           gAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTa                                                                                                             <  1:478936/251‑1 (MQ=255)
                                                                                                                            aaTCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTAt                                                                                                            <  1:416033/251‑1 (MQ=255)
                                                                                                                                tagtttGTG‑AAGATTATGCTAAATGCAACAGGAGTGCCCTTACAAGATCTGGTGTTCAGTGCGTACGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GCTTGTC‑GTACACAGCCTGCTTAGTGGCTGGGTCTACGCCGGTGACATCTGGCATACCTTGTTAATGGATTTATCGCTGTTGTCGATTTGCGTATGCATTGCTCTGGCAATACTGGTTGAGTAGTAACTAtata                                                                                                         >  2:720832/4‑247 (MQ=255)
                                                                                                                                          aaGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAAc                                                                                                >  2:464586/1‑250 (MQ=255)
                                                                                                                                           agtttagtctcatgtcaacattaTTGCCCTGACAAGATCTGGGGTTCGGTTCGTCCTTCTACTTTCCACCGTTTTTCAAAGCATTCGCTTGTGGAATCGTCATCTG‑GTTTGTC‑GTCCACCGCCTTCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCCGGCAATACTGATTGCGTGCTAACTATGTCAAGTAAAActa                                                                                              <  2:648356/229‑3 (MQ=255)
                                                                                                                                              ttATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATCCTGATTGCGTGGTAACTATGTCAATTAAAACAAtt                                                                                            >  2:170512/1‑250 (MQ=255)
                                                                                                                                              ttATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAAtt                                                                                            <  1:648990/250‑1 (MQ=255)
                                                                                                                                                  gCTCAATGCAACAGGATTGCCCTTACAAGCTCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTaaaa                                                                                                 >  2:49380/1‑241 (MQ=255)
                                                                                                                                                  gCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTaaaa                                                                                                 <  1:49380/241‑1 (MQ=255)
                                                                                                                                                  gCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAACCAATTAAGt                                                                                        >  1:566477/1‑250 (MQ=255)
                                                                                                                                                       aTGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTtct                                                                                  <  1:665546/251‑1 (MQ=255)
                                                                                                                                                         tcaaaagtctggcctttataaGATCTGTGGTTTGGTGCGGCTGTCTTCTTTCCTCGGGTTGTTAAAGCATTCGTGTTTGGATGCGTCATATG‑GTTTTTCTGT‑CACCGCCTTCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTTTTAATGGATTTATCGCTGTTTGCGAGTTTCGTTTGCCATGTTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTtca                                                                                  <  2:236698/230‑2 (MQ=255)
                                                                                                                                                         gCAACAGGATTGCCCTTACAAGATCTGGGGTTCGGGGCGTCCGTCTACTTTCCTCCGGTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTAAACCGCCTGCGTCGTGGCTGGATCTACGCCGGTGACATCTGGCAGCCCTTGTTAAGGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTca                                                                                <  1:442796/251‑1 (MQ=255)
                                                                                                                                                              aGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAAt                                                                           >  1:99911/1‑251 (MQ=255)
                                                                                                                                                                       cTTACAAGATCTGGTGTTCGGAGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGc                                                                  >  1:679187/1‑251 (MQ=255)
                                                                                                                                                                          aCAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCTTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTAGTTCTCAACAATCATTGTAGCGGt                                                               <  2:81437/251‑1 (MQ=255)
                                                                                                                                                                             aGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGt                                                            >  2:523985/1‑251 (MQ=255)
                                                                                                                                                                              gATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGtt                                                           >  2:410447/1‑251 (MQ=255)
                                                                                                                                                                              gATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGtt                                                           >  1:789992/1‑251 (MQ=255)
                                                                                                                                                                               aTCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTg                                                          >  2:392709/1‑251 (MQ=255)
                                                                                                                                                                                             gCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCgg                                             >  2:507651/1‑250 (MQ=255)
                                                                                                                                                                                               gtccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACACTCATTGTAGCGGTAGTTTCGGTTCTTGCCggctg                                          >  2:874758/1‑248 (MQ=255)
                                                                                                                                                                                               gtccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtg                                          >  1:820452/1‑251 (MQ=255)
                                                                                                                                                                                                 ccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GCTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtgg                                         <  1:390944/250‑1 (MQ=255)
                                                                                                                                                                                                    tcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtgtgagc                                     >  2:160087/1‑246 (MQ=255)
                                                                                                                                                                                                       aCTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCt                                    <  1:170512/249‑1 (MQ=255)
                                                                                                                                                                                                                     ttCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATg                    >  1:823195/1‑251 (MQ=255)
                                                                                                                                                                                                                           gCATTCGCTTTTGGATTCTGCATCTG‑GTTGGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCGTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTa                                                                                                    <  2:517076/165‑1 (MQ=255)
                                                                                                                                                                                                                           gCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTa                                                                                                    >  1:517076/1‑165 (MQ=255)
                                                                                                                                                                                                                           gCATTCGCGTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAAtct              <  1:734710/251‑1 (MQ=255)
                                                                                                                                                                                                                            cATTCGCGTTTGGATTCGTAATCTG‑GTTTGTA‑GTACACCGCCTGATTCGTGGCTGGATCTACGCCGGAGAAAACTGTCATACCTTGTGAATGGATTTATCGCTGTTTGAGATTAGCGTTTGACTTGCTCTGGCAATAATGCTTGCGTGGTAACTATGTCAGTTAACACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGAGGTTCTTACCGGATGGCAGCAGTGGAATTTTTACATGAAAtctt             >  1:490941/1‑250 (MQ=255)
                                                                                                                                                                                                                               tCGCGTATGGATTCGTCAGCTG‑GTTTGTC‑GTACACCGCCTGCTTCGAGGATGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGTCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGTTGGTGGCTGTGGAATTATTACGTGCAATATCCgg          >  2:138878/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gTTTGGATTCGTCATCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGac      >  1:174922/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggATTCGTCACCTG‑GTTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACAcgc  <  1:132084/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
ATTGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTG‑AAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTG‑GCTTGTC‑GTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGC  >  NC_000913/1720891‑1721371

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