Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 540,546 C→T 12.1% intergenic (+38/‑19) allB → / → ybbY allantoinase/putative uracil/xanthine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913540,5460CT12.1% 168.8 / 3.3 73intergenic (+38/‑19)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base C (17/45);  new base T (6/3);  total (24/49)
Fisher's exact test for biased strand distribution p-value = 5.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGG  >  NC_000913/540308‑540778
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
ttCGTTTTCCTTTAGCCGAATAGCAGGGTTGTTCTTACCAATGCCTATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCAGACCAGTGGCGTGCGTATCACGAAAACCATCTTGCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑Ttc                                                                                                                                                                                                                                               <  1:320190/239‑1 (MQ=255)
ttCGTCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑Ttc                                                                                                                                                                                                                                               <  2:287940/239‑1 (MQ=255)
ttCGGCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAAGTACGATCTGGAATAGCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAATCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑Ttc                                                                                                                                                                                                                                               <  2:483997/239‑1 (MQ=255)
        tattcaGCCGAATAGCAGTTATGTTCTTACCAATTACGATCTGGAATATCGCCACAAAGACAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTa                                                                                                                                                                                                                                      <  1:353385/238‑1 (MQ=255)
        cattccCCCGAATAGCAGCTTTGTTCTTACCAAGGCCGTTCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTa                                                                                                                                                                                                                                      <  1:121293/239‑1 (MQ=255)
        cattcaGCGGAATGGCAGCTTTTTTCTTACCAATGGCGATCGGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCCTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTa                                                                                                                                                                                                                                      <  2:1025519/239‑1 (MQ=255)
        cattcaGCCGAATAGCCCCTATGTCTGTCCGAATGTCTGTCTGGAATATCGCCACAAAGTCATGCCTTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTa                                                                                                                                                                                                                                      <  1:676556/239‑1 (MQ=255)
        cattcaGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTa                                                                                                                                                                                                                                      <  2:173022/239‑1 (MQ=255)
        cattcaGCCGAAATACACTTATGTTCTTACCAATGCCGATCTGGAATATCGCCACAAAGTCAGCCCGGAGGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTa                                                                                                                                                                                                                                      <  1:877885/239‑1 (MQ=255)
          ttcaGCCGAATAGCAGCTATGTTCTTCCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAg                                                                                                                                                                                                                                    <  2:634566/239‑1 (MQ=255)
           tcaGCCGAATAGCAGCTAGGTTCTTTCCAAGGACGATCTGGAATATCGCCAAAAAGGCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAgg                                                                                                                                                                                                                                   <  2:854905/239‑1 (MQ=255)
            caGCCGAATAGCAGCTATGTTCCTACCAATGACGATCTGGAATATCCCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGgt                                                                                                                                                                                                                                  <  2:572210/239‑1 (MQ=255)
                                         tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGCAAGGCCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCC‑TTTTCC‑GGTTTAAGTGG‑GTTTTTGTTCCAtttttgcgttcgcg                                                                                                                                                                                                     >  2:834427/1‑228 (MQ=255)
                                         tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGGGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGACGTAATCTGGGCCCTGCCATGCCCGTCCTTGCGGCGGGGATTTTTCC‑GATAAAGGGGG‑TTTTTTTTTCATTTTTGCAggagccg                                                                                                                                                                                                     >  1:809357/1‑233 (MQ=255)
                                           gcctatgggAATATCGGCCATAATTCAGGCCGTATGTTTGGCGTACCATTGGGGAGCGTATCACGAAAGCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCATTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCcgc‑g                                                                                                                                                                                                <  1:128007/232‑1 (MQ=255)
                                             gATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑Gaa                                                                                                                                                                                              <  2:448170/239‑1 (MQ=255)
                                              aTCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGC‑TTTTCC‑GGGTTAGGGGG‑GTTTATTTTGAATTTTTCAGCTA‑GCCCC‑Gaaa                                                                                                                                                                                             >  2:31975/1‑239 (MQ=255)
                                                cTGGAATATCGCCAAAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGc                                                                                                                                                                                           <  2:967691/239‑1 (MQ=255)
                                                       agcGCCACAAATTCATCCCGTATGTTTGCCGTACCATTGGCGCGCGTAACACGAAAACCATCTTACGGGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAAGCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTAt                                                                                                                                                                                    <  2:809357/237‑1 (MQ=255)
                                                         cGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCCGGCCCCTGCAATGCCCGTCCTTGCGGCGGACA‑ATATCC‑CGGTCAGGGGG‑GTTTATGTTAAATTTTTTAAGCA‑TCCGC‑GGAAGccctgtttca                                                                                                                                                                                  >  1:562041/1‑231 (MQ=255)
                                                                aaGTCAGCCCGTTTGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTc                                                                                                                                                                          <  2:281732/239‑1 (MQ=255)
                                                                      gCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGAGG‑GTTTGTGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTTTCA‑GGGTTTCatgggg                                                                                                                                                                    >  2:397691/1‑234 (MQ=255)
                                                                       cccGTATGTTGGCCGTACCATTGGCGCGGGTTTCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGCCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGTCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGtt                                                                                                                                                                   <  2:809444/239‑1 (MQ=255)
                                                                           tATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATC‑TACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCTGACA‑TTCTCC‑GGTTAAGGGGTGGTTTAT‑TTCAATTTTTCAGTCA‑GCCGC‑GAAAACCCGTTATTA‑GGATTTCATGGGTTTTTCt                                                                                                                                                              >  2:504444/1‑239 (MQ=255)
                                                                            aTGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTTCGGCGGGCA‑TTCTCC‑GGTTTAGGGGG‑GTTTTTGTTTCATTTTTCAGTCA‑GCCCC‑CCAAACCTGTTTTTA‑GGAGTTTAGTGGgtttttt                                                                                                                                                              >  1:591378/1‑237 (MQ=255)
                                                                                ggcctttccctttGGCGGGCG‑ACCACGAAAACCCGCTTACGGGGTGATGTTATTTACGACATTTAACCGGGCTTCCATGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTAtt                                                                                                                                                         <  1:808662/228‑1 (MQ=255)
                                                                                     gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTAttttt                                                                                                                                                      <  2:228955/238‑1 (MQ=255)
                                                                                           ttGGCGCGCGTTGCACGAAAACCATCTTACGTGGTGATGTGATTTACGGCATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                               <  2:452637/239‑1 (MQ=255)
                                                                                           ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCC‑TTTTCCGGGTTAAGG‑GG‑GCTTTTGTTTCATTTTTCGGTCA‑GCCCG‑GAAAACAGGTTATTA‑GGGGTTTCGAGGGTTTTTTTTTTTTTTTGctaaa                                                                                                                                               >  1:228955/1‑235 (MQ=255)
                                                                                           ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCC‑ATCTTC‑CGGTAAGGTGT‑GTTTTTGTTCCATTTTTCCGTGC‑GCCGC‑GAAAGCCCGGTTTCAGGGTTTTGAGGGGTttttttttttttttgacaaa                                                                                                                                                >  1:377637/1‑233 (MQ=255)
                                                                                           ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                               <  2:562041/239‑1 (MQ=255)
                                                                                            tGGCGCGCGTTTCACTAAAACCATCTTACGTGGGGAGTTGATTTACCACAGTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGGTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac                                                                                                                                              <  1:407283/239‑1 (MQ=255)
                                                                                                  gcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCaa                                                                                                                                                 <  2:752731/230‑1 (MQ=255)
                                                                                                      ttcgcGAAAGCCATCTTACGTGGTTATTTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                                                                                                     <  2:717188/234‑1 (MQ=255)
                                                                                                            aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGGTAAGGTGT‑GTTTATGTTTCATTTTTCTGTCA‑GCCGC‑CAAAACCTGTTTTCA‑GGGTTTTAGTGGGTTTTTTTTTTTTTTTGCAACAAGGGTGTgggtggttct                                                                                                                              >  1:752731/1‑231 (MQ=255)
                                                                                                                  aTCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGG‑GGTTTTGTTTAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGGTTTCAGTGGTTTTTTTTTTTTTTTTTTAACACCGTTTTGGTTTCCCCTAcccca                                                                                                                        >  2:500929/1‑235 (MQ=255)
                                                                                                                         ggtgtgatgttatTTCACACATTGAACGGGGCTTCCCTGTTGCGCCGAAAGTTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTg                                                                                                                 <  1:31975/236‑1 (MQ=255)
                                                                                                                               atgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGGGT‑GTTTATGTTCAATTTTGCGGTCA‑GCCGC‑GAAAGCCTGTTTTCA‑GGATTTCACGGGGTTTTCTTTATTTTTTTGAACACGGTTGTGGTTCTTCCTACCCGACTTTTTGCttttt                                                                                                           >  2:108902/1‑238 (MQ=255)
                                                                                                                                  gtatttacactgtggaaacgggCTTCCCTGTTGCGCCGAAAGGTCAATTTAGCCTTAAACACAAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGt                                                                                                        <  1:39206/220‑1 (MQ=255)
                                                                                                                                       tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGGTTTTCTTTATTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGcc                                                                                                    >  2:999791/1‑238 (MQ=255)
                                                                                                                                              ttGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                           <  2:346120/208‑1 (MQ=255)
                                                                                                                                                taccaggctttCCCGGTTGCCCCGAAAGGTCAATTTTTCCTTAAACCTCCGCAGTAATCTGGCCCATGCAATGCCCGTCCTTGCGGCGGGCA‑TTGTCC‑GGTTAATGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                          <  1:1023516/230‑1 (MQ=255)
                                                                                                                                                taacctgggttCCCTTGTGGGCCGAAAGGTCAATTTTTCCTTAAACATCAGCAGTAATCTGGCCTCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                          <  1:500929/230‑1 (MQ=255)
                                                                                                                                                       gCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGGGT‑GTTTTTGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACCCTACTTTTTGCTTTTCAGTTGCCGCAAAGTAGCCTGCtt                                                                                    >  2:471421/1‑238 (MQ=255)
                                                                                                                                                       gCTTCCCTGTTGCGCAGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCt                                                                                     <  2:899792/237‑1 (MQ=255)
                                                                                                                                                            ccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACCTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTc                                                                              <  2:194568/239‑1 (MQ=255)
                                                                                                                                                             cttttggtcatgaaagtgaaTTTTGCCTTCAACATCAGCAGTAATGTGGCCCCGGCCATGCCCGTCCGTGCGGCTGGCA‑TTCGCG‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCC‑GCCGC‑GAAAGCCTGTTAGCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGGTGTGGTTCCTCCTACGCTACTTTCTTCTTTTCAGTTGCCGCAAAGTAGGCTGCTTACGCTCa                                                                             <  1:916032/221‑1 (MQ=255)
                                                                                                                                                                    ccccaaaGGTCTTTTTTTCTTTAAACATCAGCAGTAATCTGGCCCGTGCAATGCCCGTCCGTTCGGCGGCCC‑TTTCAC‑GGGTGAGGAGG‑GGTTAGGTTACATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCCACGCTGCTTTCTGCTTTTCAGTTGCCGCAAAGTAGCTTGCTTACGCTCACGCAata                                                                      <  1:504444/235‑1 (MQ=255)
                                                                                                                                                                              cAAATTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGGCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCt                                                             <  2:26583/238‑1 (MQ=255)
                                                                                                                                                                                  ttATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTa                                                        <  2:591378/239‑1 (MQ=255)
                                                                                                                                                                                       cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCa                                                                         <  2:505970/217‑1 (MQ=255)
                                                                                                                                                                                       cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                   <  2:540941/239‑1 (MQ=255)
                                                                                                                                                                                           aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑Ttctc                                                                                                                                                                                                                                             >  2:893793/1‑54 (MQ=255)
                                                                                                                                                                                           aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑Ttctc                                                                                                                                                                                                                                             <  1:893793/54‑1 (MQ=255)
                                                                                                                                                                                           aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGGGT‑GTTTTTGTTCAATTTTGCAGTCA‑GCCGC‑GGAAGCCTGTTTTCA‑GGATTTCAGTGGTTTTTTTTTATTTTTTGCCACACGGTTTTGGTTCCTCCTACGCTACTTTTTTCTTTTCCGTTTCCCCCAAAGAACCCGCTTACCCTCCCGCCAAATGCCTTTTTTTCTccccccccgg                                               >  1:253071/1‑232 (MQ=255)
                                                                                                                                                                                              aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCt                                                             <  2:290117/222‑1 (MQ=255)
                                                                                                                                                                                                cagcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTTTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGGTTTTCTTTTTTTTTTTCAACACGGTTGTGGGTTCTCCTACCCTACTTTCTGCTTTTTAGTTGTCGCAAAGATGCCTGCTTTCTCTCACGCAATATGCTTTTTTTTCTTCCGCACTGGCCt                                            >  2:67718/1‑237 (MQ=255)
                                                                                                                                                                                                      tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGgt                                                                                                                                                                                                                                  <  1:365818/53‑1 (MQ=255)
                                                                                                                                                                                                      tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGgt                                                                                                                                                                                                                                  >  2:365818/1‑53 (MQ=255)
                                                                                                                                                                                                       aaTCTGGCCCCTGACATTCCCTTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTc                                   <  2:184135/239‑1 (MQ=255)
                                                                                                                                                                                                        aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTTAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGc                                                                                                                                                                                           >  2:660662/1‑87 (MQ=255)
                                                                                                                                                                                                        aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGc                                                                                                                                                                                           <  1:660662/87‑1 (MQ=255)
                                                                                                                                                                                                        aTCTGGCCCCTGAAATGCCCGTCCGTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCa                                                                             <  2:830714/196‑1 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTTTCC‑GGTTAAGGGGT‑GGTTTTGGTCAATTTTGCAGTTA‑GCCGCGGAAAGCC‑GTTATCA‑GGGTTTTAGTGGTTTTTTTTTTTTTTTTTGTACACCGGGTTTGTTTCTCCCACCCGCCTTTATCCTTTTTCGTTGGCCCAAAGTAGGCCGCTTTCGCTCACGCCAAATCTGTTTTTTTTCACCCCACCGGCCTGCTTCCCCCAGGGGt                             >  1:772325/1‑239 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCCATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                          >  2:285115/1‑130 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAac                                                                                                                                                <  2:582755/124‑1 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                          <  1:285115/130‑1 (MQ=255)
                                                                                                                                                                                                                  tGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTTTCA‑GGGTTTCAGTGGTTTTTCTTTTTTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTCCTTTTTGCTTTTTAGTTTCCGCCAAGTAGCCTGCTTACGCTCACGCCATATTCTTTTTTTTTTACCCCCCTGGCCTCCTTCGCTCAGGGGTGTTgg                        >  1:1018101/1‑238 (MQ=255)
                                                                                                                                                                                                                   gCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCa                                                                                                                                                                         <  2:806870/93‑1 (MQ=255)
                                                                                                                                                                                                                       tGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGGTTAGGTGT‑GTTTATGTTCAATTTTGCCGTCA‑GCCGC‑GAAAGCCTGGTATTA‑GGGTTTCCGTGGTTTTTTTTTATTTTTTTCAACACGGTTGTgg                                                                                                                                      >  1:582755/1‑124 (MQ=255)
                                                                                                                                                                                                                             ttcTTGCGGCGGTCA‑TTTTCC‑GGTTAAGGTTT‑GTTTATGTTCAATTTTGCCGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg             <  2:901329/237‑1 (MQ=255)
                                                                                                                                                                                                                             tCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg             <  2:1018101/239‑1 (MQ=255)
                                                                                                                                                                                                                              ttttGGGGGGGGCT‑TTC‑CC‑GGTTTAGTTGT‑GTTTTTGTTCAATTTTGC‑GTGAGGCCGC‑GAAAGCCTGTTTTCA‑GGATTTCAGTTGTTTTTCTTTATTTTTTGCAACCCGGTTGGGGTTCCTCTTACGCTACTTTCTGCTTTTCAGTGGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg           <  1:999791/237‑1 (MQ=255)
                                                                                                                                                                                                                              ccTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGGGT‑TTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACTGTTGTGGTTCCTCCTACTCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg             <  2:138527/238‑1 (MQ=255)
                                                                                                                                                                                                                               cTTGCGGCGGGCA‑TTTGCC‑TGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg           <  2:272999/239‑1 (MQ=255)
                                                                                                                                                                                                                               cTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg           <  2:377637/239‑1 (MQ=255)
                                                                                                                                                                                                                                        ggCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTTCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCCCTGGCCTGCTTCGCTCAGGCGTTTTGCGGTTATCGTCCCGCtattattg  >  2:510437/1‑237 (MQ=255)
                                                                                                                                                                                                                                        ggCA‑TTCTCC‑GGTTAAGGGGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTTGCCTGCTTCGCTCAGGCGTTTTTCGGGCATCTTCGCGCcatttttg  >  2:695521/1‑235 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
TTCGTCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCA‑TTCTCC‑GGTTAAGGTGT‑GTTTATGTTCAATTTTGCAGTCA‑GCCGC‑GAAAGCCTGTTATCA‑GGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGG  >  NC_000913/540308‑540778

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