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

ATTGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGC  >  NC_000913/1720891‑1721371
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
aTTGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGt                                                                                                                                                                                                                                          <  1:642361/250‑1 (MQ=255)
  tGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGt                                                                                                                                                                                                                                      >  2:297149/1‑251 (MQ=255)
  tGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGGTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCg                                                                                                                                                                                                                                       <  1:720980/250‑1 (MQ=255)
     tAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGCAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAAAAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTac                                                                                                                                                                                                                                    >  1:744487/1‑250 (MQ=255)
                     ctgcATTTGGTTAGCGTAAGACATATTTTTAAGAAATTGCGCGTAAATTTTTCGAATTTAAAGAATATTTATTAATCTAACGCAATATATTTGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTAAATGCAACAGGATTGTCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTAAAAGCATTCGCGTTTGGATTCGTAATATGGTTTGTGTTACACTGCCTGATTTGTggt                                                                                                                                                                                                                   >  2:382710/1‑250 (MQ=255)
                       gcATTTGGTTACCGTCAGACATATTTTCCAGAAATTGCGCGTAAATTTTTCGCCCTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGACGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGGCCACTTTCCTCCGTTTTTCAACGCATTCGCGTGTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGc                                                                                                                                                                                                                   <  1:33996/249‑1 (MQ=255)
                         aTTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGTTTTGTCGTACACCGCCTGCTTCGTGGCTGGa                                                                                                                                                                                                               <  1:776959/251‑1 (MQ=255)
                            tGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATc                                                                                                                                                                                                             <  1:667533/250‑1 (MQ=255)
                                  ccGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCGAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTTCTTTCCTCCGTTTTTCAAAGCTTTCGCGTTTGGATTCGTCAGCTGGTTTGTCGTACACCGCCTGCGTCGTGGCTGGATCTACGcc                                                                                                                                                                                                       <  2:358635/250‑1 (MQ=255)
                                    gTCAGACAGATTTGTCATAAATTGCGTGTAAAGTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGTTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGt                                                                                                                                                                                                    <  2:232324/251‑1 (MQ=255)
                                      tagccATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGAAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTGATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGAGTTCGGCGAGTCCGTCTAATGTCCTCCGTTTTTCAAAGCATTCTCGATTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTAGTGGCTGAATCTACGcagg                                                                                                                                                                                                     >  2:283986/5‑245 (MQ=255)
                                       aGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGAc                                                                                                                                                                                                 >  1:863212/1‑251 (MQ=255)
                                           atatTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACTCCGGTGACATCTg                                                                                                                                                                                            >  1:264840/1‑251 (MQ=255)
                                                       taTTGCGCGAAAATTTTTCGAACATAAAGAATATTTATTAATATAACAAAATATATTCTGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCACAATGCAACAGGCATGCACTTACAAGATATGGTGTTCGGTGCGTACGACTACTTTCTTCCGTTTTTCAAAGCATTCGCGTTTGGATGCGTCATCTGGTTTGTCGATCACCGCCTGCTTCGTGTCGGGATATACACCGGTTACATGTGGCATCCATTGt                                                                                                                                                                                 >  2:518866/2‑251 (MQ=255)
                                                       aaTTGCGCGGAAATTTTGCGCACTTAAATAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATTTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGCTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGCCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGt                                                                                                                                                                                 <  1:818631/251‑1 (MQ=255)
                                                                 aaaTTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGAttt                                                                                                                                                                       <  2:332657/251‑1 (MQ=255)
                                                                   aTTTTTCGCACTTAAAGACTATTTATTAATCTACCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGTTCTGGTGTTCGGTTCGTCCGTCTACTTTCCTCCGTTTTTCATATCATTCGCGATTGGATTCGTCCTCTGTTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTATGTCGGTGTCATCTGGCATCCCTTGTTAATGTAttttt                                                                                                                                                                     >  2:279333/1‑249 (MQ=255)
                                                                      tttCGCACTTAAAGAATATTTATTAACCGAACTAAATATATTCGGTCGTAAAGGAAACTACTTTGTGAAGTTTATGCCCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCg                                                                                                                                                                   <  2:413201/250‑1 (MQ=255)
                                                                         cgcaatgaaaGAATATTAATTAATCATACGCTATAGACTCGGTCGGAAAGGAATCTTCTGTGTGAAGTTTACTCTCAACGAATCAGGATTCCCCTTACAAGATCGGTTGTGCTGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTTGTTTGTCTAACACCGCCTGCTTCGTGGCTTGATCTACGCCGGTGACATCTGGCATCCCTTTTTAAGGGATTTATCGCTGt                                                                                                                                                               <  2:252453/244‑1 (MQ=255)
                                                                          gCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGAATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGtt                                                                                                                                                              <  1:775938/251‑1 (MQ=255)
                                                                                  ataattttgaatattttaacgacatatatTCGGTCGTAAAGTACTCTACTTTTTGAAGTTTCTGCTCAATGCCACAGGATTGCCCTTACAAACTCTGGTGTTCGGTTCGTCCTGCTACTTACCTCCGTTTTTCAAAGCATTCGCTTTTTGATTCGCCATTGCGTTTTTCCGACTCCGCCTGCTTCGTCGCTGGATCTACGCCGTTGTCATCTGGCATCCCTTGTTAATGGATTTATCGCTGTGTGCGAtct                                                                                                                                                      <  2:58902/228‑3 (MQ=255)
                                                                                   gAATATTTATTAATCTAACACAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTAAATGCAACAGGATTGCCCTTACAAGTTCTGGTGTTTGTTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGGGACATATGGCATCCCTTGTTAATGGTTTTTTCCCTGTTTgagctttg                                                                                                                                                     >  2:803940/1‑244 (MQ=255)
                                                                                      tatttattAATATAACGCAATATATTCGGACGTAAAGGAATCTACTTTATGAAGTTTATGCTCAATGCAAAAGGATTGCCCTTACAAGATCTGGTGTTCGGTGAGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGCTTCGTCATCTGTTTTGTCGTACACCGCCTGCTTCGTAGCTGACTCTACGCCGGTGCCATCTGGCATCACTTGTTAATGGATTTATCGCTGTTTGCGagtatcgt                                                                                                                                                  >  2:719536/1‑244 (MQ=255)
                                                                                       tattcttaATTTAACGTAATATATTCGGTCGTAAAGGCATCTCCTTTTGTAATGTTATGTCCAATGCAACAGGATTGCCCTTACAAGATATGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATGTATCGCTGTTTGCGATTTGCGt                                                                                                                                                  <  2:393780/245‑1 (MQ=255)
                                                                                         ttattaATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGttt                                                                                                                                                <  1:11621/250‑1 (MQ=255)
                                                                                              aaTCTAACGCAAGATAGTCGGTCGTAAAGGAATCTTTTTGTTTAATTTTATGCTCAGTGCAGTCGGATTGCCCTTACAAGATTTGTTGTTTGGTTCGTCCGTTTGCTTTCCTCCGTTTTACAAAGCATCCTCGTTTGTATTCGTCCTCTGGTGTGTCGTACCCCGCCTGCTTCGTGGCTGGATCTGCGCCGGTGACATCTGGCATCCCTTGTTAATGGAATTATCGCTGTTTGCGATTtgcctttgcctt                                                                                                                                           <  2:767606/250‑1 (MQ=255)
                                                                                                 cTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGAGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGctc                                                                                                                                       >  1:428671/1‑251 (MQ=255)
                                                                                                    aCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTCTGCTCAATGAAACAGGATTGCCCTTCCAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGCTCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGCTTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTgg                                                                                                                                    >  2:852276/1‑251 (MQ=255)
                                                                                                         atatatTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGTTGACATCTGGCATCCCTTGTTAAt                                                                                                                                                                             >  2:22458/1‑205 (MQ=255)
                                                                                                         atatatTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAAt                                                                                                                                                                             <  1:22458/205‑1 (MQ=255)
                                                                                                             atTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCACTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGAGCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTg                                                                                                                            >  2:775791/1‑250 (MQ=255)
                                                                                                                cGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTg                                                                                                                        >  2:734710/1‑251 (MQ=255)
                                                                                                                      tAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGt                                                                                                                  >  1:89195/1‑251 (MQ=255)
                                                                                                                           gAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTa                                                                                                             <  1:478936/251‑1 (MQ=255)
                                                                                                                            aaTCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTAt                                                                                                            <  1:416033/251‑1 (MQ=255)
                                                                                                                                tagtttGTGAAGATTATGCTAAATGCAACAGGAGTGCCCTTACAAGATCTGGTGTTCAGTGCGTACGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACAGCCTGCTTAGTGGCTGGGTCTACGCCGGTGACATCTGGCATACCTTGTTAATGGATTTATCGCTGTTGTCGATTTGCGTATGCATTGCTCTGGCAATACTGGTTGAGTAGTAACTAtata                                                                                                         >  2:720832/4‑247 (MQ=255)
                                                                                                                                          aaGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAAc                                                                                                >  2:464586/1‑250 (MQ=255)
                                                                                                                                           agtttagtctcatgtcaacattaTTGCCCTGACAAGATCTGGGGTTCGGTTCGTCCTTCTACTTTCCACCGTTTTTCAAAGCATTCGCTTGTGGAATCGTCATCTGGTTTGTCGTCCACCGCCTTCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCCGGCAATACTGATTGCGTGCTAACTATGTCAAGTAAAActa                                                                                              <  2:648356/229‑3 (MQ=255)
                                                                                                                                              ttATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATCCTGATTGCGTGGTAACTATGTCAATTAAAACAAtt                                                                                            >  2:170512/1‑250 (MQ=255)
                                                                                                                                              ttATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAAtt                                                                                            <  1:648990/250‑1 (MQ=255)
                                                                                                                                                  gCTCAATGCAACAGGATTGCCCTTACAAGCTCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTaaaa                                                                                                 >  2:49380/1‑241 (MQ=255)
                                                                                                                                                  gCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTaaaa                                                                                                 <  1:49380/241‑1 (MQ=255)
                                                                                                                                                  gCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAACCAATTAAGt                                                                                        >  1:566477/1‑250 (MQ=255)
                                                                                                                                                       aTGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTtct                                                                                  <  1:665546/251‑1 (MQ=255)
                                                                                                                                                         tcaaaagtctggcctttataaGATCTGTGGTTTGGTGCGGCTGTCTTCTTTCCTCGGGTTGTTAAAGCATTCGTGTTTGGATGCGTCATATGGTTTTTCTGTCACCGCCTTCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTTTTAATGGATTTATCGCTGTTTGCGAGTTTCGTTTGCCATGTTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTtca                                                                                  <  2:236698/230‑2 (MQ=255)
                                                                                                                                                         gCAACAGGATTGCCCTTACAAGATCTGGGGTTCGGGGCGTCCGTCTACTTTCCTCCGGTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTAAACCGCCTGCGTCGTGGCTGGATCTACGCCGGTGACATCTGGCAGCCCTTGTTAAGGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTca                                                                                <  1:442796/251‑1 (MQ=255)
                                                                                                                                                              aGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAAt                                                                           >  1:99911/1‑251 (MQ=255)
                                                                                                                                                                       cTTACAAGATCTGGTGTTCGGAGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGc                                                                  >  1:679187/1‑251 (MQ=255)
                                                                                                                                                                          aCAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCTTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTAGTTCTCAACAATCATTGTAGCGGt                                                               <  2:81437/251‑1 (MQ=255)
                                                                                                                                                                             aGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGt                                                            >  2:523985/1‑251 (MQ=255)
                                                                                                                                                                              gATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGtt                                                           >  2:410447/1‑251 (MQ=255)
                                                                                                                                                                              gATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGtt                                                           >  1:789992/1‑251 (MQ=255)
                                                                                                                                                                               aTCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTg                                                          >  2:392709/1‑251 (MQ=255)
                                                                                                                                                                                             gCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCgg                                             >  2:507651/1‑250 (MQ=255)
                                                                                                                                                                                               gtccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACACTCATTGTAGCGGTAGTTTCGGTTCTTGCCggctg                                          >  2:874758/1‑248 (MQ=255)
                                                                                                                                                                                               gtccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtg                                          >  1:820452/1‑251 (MQ=255)
                                                                                                                                                                                                 ccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtgg                                         <  1:390944/250‑1 (MQ=255)
                                                                                                                                                                                                    tcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtgtgagc                                     >  2:160087/1‑246 (MQ=255)
                                                                                                                                                                                                       aCTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCt                                    <  1:170512/249‑1 (MQ=255)
                                                                                                                                                                                                                     ttCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATg                    >  1:823195/1‑251 (MQ=255)
                                                                                                                                                                                                                           gCATTCGCTTTTGGATTCTGCATCTGGTTGGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCGTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTa                                                                                                    <  2:517076/165‑1 (MQ=255)
                                                                                                                                                                                                                           gCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTa                                                                                                    >  1:517076/1‑165 (MQ=255)
                                                                                                                                                                                                                           gCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAAtct              <  1:734710/251‑1 (MQ=255)
                                                                                                                                                                                                                            cATTCGCGTTTGGATTCGTAATCTGGTTTGTAGTACACCGCCTGATTCGTGGCTGGATCTACGCCGGAGAAAACTGTCATACCTTGTGAATGGATTTATCGCTGTTTGAGATTAGCGTTTGACTTGCTCTGGCAATAATGCTTGCGTGGTAACTATGTCAGTTAACACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGAGGTTCTTACCGGATGGCAGCAGTGGAATTTTTACATGAAAtctt             >  1:490941/1‑250 (MQ=255)
                                                                                                                                                                                                                               tCGCGTATGGATTCGTCAGCTGGTTTGTCGTACACCGCCTGCTTCGAGGATGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGTCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGTTGGTGGCTGTGGAATTATTACGTGCAATATCCgg          >  2:138878/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGac      >  1:174922/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggATTCGTCACCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACAcgc  <  1:132084/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
ATTGCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGC  >  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: