Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 540,546 C→T 12.5% 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.5% 165.0 / 3.5 71intergenic (+38/‑19)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base C (17/43);  new base T (6/3);  total (24/47)
Fisher's exact test for biased strand distribution p-value = 5.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGG  >  NC_000913/540308‑540778
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
ttCGTTTTCCTTTAGCCGAATAGCAGGGTTGTTCTTACCAATGCCTATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCAGACCAGTGGCGTGCGTATCACGAAAACCATCTTGCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATtc                                                                                                                                                                                                                                               <  2:320190/239‑1 (MQ=255)
ttCGTCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATtc                                                                                                                                                                                                                                               <  1:287940/239‑1 (MQ=255)
ttCGGCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAAGTACGATCTGGAATAGCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAATCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATtc                                                                                                                                                                                                                                               <  1:483997/239‑1 (MQ=255)
        tattcaGCCGAATAGCAGTTATGTTCTTACCAATTACGATCTGGAATATCGCCACAAAGACAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                                                                                                                                                                      <  2:353385/238‑1 (MQ=255)
        cattccCCCGAATAGCAGCTTTGTTCTTACCAAGGCCGTTCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                                                                                                                                                                      <  2:121293/239‑1 (MQ=255)
        cattcaGCGGAATGGCAGCTTTTTTCTTACCAATGGCGATCGGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCCTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                                                                                                                                                                      <  1:1025519/239‑1 (MQ=255)
        cattcaGCCGAATAGCCCCTATGTCTGTCCGAATGTCTGTCTGGAATATCGCCACAAAGTCATGCCTTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                                                                                                                                                                      <  2:676556/239‑1 (MQ=255)
        cattcaGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                                                                                                                                                                      <  1:173022/239‑1 (MQ=255)
        cattcaGCCGAAATACACTTATGTTCTTACCAATGCCGATCTGGAATATCGCCACAAAGTCAGCCCGGAGGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                                                                                                                                                                      <  2:877885/239‑1 (MQ=255)
          ttcaGCCGAATAGCAGCTATGTTCTTCCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAg                                                                                                                                                                                                                                    <  1:634566/239‑1 (MQ=255)
           tcaGCCGAATAGCAGCTAGGTTCTTTCCAAGGACGATCTGGAATATCGCCAAAAAGGCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAgg                                                                                                                                                                                                                                   <  1:854905/239‑1 (MQ=255)
            caGCCGAATAGCAGCTATGTTCCTACCAATGACGATCTGGAATATCCCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                                                                                                                                                                  <  1:572210/239‑1 (MQ=255)
                                         tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGCAAGGCCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTTCCGGTTTAAGTGGGTTTTTGTTCCAtttttgcgttcgcg                                                                                                                                                                                                     >  1:834427/1‑228 (MQ=255)
                                         tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGGGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGACGTAATCTGGGCCCTGCCATGCCCGTCCTTGCGGCGGGGATTTTTCCGATAAAGGGGGTTTTTTTTTCATTTTTGCAggagccg                                                                                                                                                                                                     >  2:809357/1‑233 (MQ=255)
                                           gcctatgggAATATCGGCCATAATTCAGGCCGTATGTTTGGCGTACCATTGGGGAGCGTATCACGAAAGCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCATTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                                                                                                                                                                <  2:128007/232‑1 (MQ=255)
                                             gATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGaa                                                                                                                                                                                              <  1:448170/239‑1 (MQ=255)
                                              aTCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGCTTTTCCGGGTTAGGGGGGTTTATTTTGAATTTTTCAGCTAGCCCCGaaa                                                                                                                                                                                             >  1:31975/1‑239 (MQ=255)
                                                cTGGAATATCGCCAAAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                                                                                                                                                           <  1:967691/239‑1 (MQ=255)
                                                       agcGCCACAAATTCATCCCGTATGTTTGCCGTACCATTGGCGCGCGTAACACGAAAACCATCTTACGGGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAAGCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                                                                                                                                                    <  1:809357/237‑1 (MQ=255)
                                                         cGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCCGGCCCCTGCAATGCCCGTCCTTGCGGCGGACAATATCCCGGTCAGGGGGGTTTATGTTAAATTTTTTAAGCATCCGCGGAAGccctgtttca                                                                                                                                                                                  >  2:562041/1‑231 (MQ=255)
                                                                aaGTCAGCCCGTTTGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTc                                                                                                                                                                          <  1:281732/239‑1 (MQ=255)
                                                                      gCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGAGGGTTTGTGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTCatgggg                                                                                                                                                                    >  1:397691/1‑234 (MQ=255)
                                                                       cccGTATGTTGGCCGTACCATTGGCGCGGGTTTCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGCCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGTCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGtt                                                                                                                                                                   <  1:809444/239‑1 (MQ=255)
                                                                           tATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCTGACATTCTCCGGTTAAGGGGTGGTTTATTTCAATTTTTCAGTCAGCCGCGAAAACCCGTTATTAGGATTTCATGGGTTTTTCt                                                                                                                                                              >  1:504444/1‑239 (MQ=255)
                                                                            aTGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTTCGGCGGGCATTCTCCGGTTTAGGGGGGTTTTTGTTTCATTTTTCAGTCAGCCCCCCAAACCTGTTTTTAGGAGTTTAGTGGgtttttt                                                                                                                                                              >  2:591378/1‑237 (MQ=255)
                                                                                ggcctttccctttGGCGGGCGACCACGAAAACCCGCTTACGGGGTGATGTTATTTACGACATTTAACCGGGCTTCCATGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAtt                                                                                                                                                         <  2:808662/228‑1 (MQ=255)
                                                                                     gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAttttt                                                                                                                                                      <  1:228955/238‑1 (MQ=255)
                                                                                           ttGGCGCGCGTTGCACGAAAACCATCTTACGTGGTGATGTGATTTACGGCATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                               <  1:452637/239‑1 (MQ=255)
                                                                                           ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTTCCGGGTTAAGGGGGCTTTTGTTTCATTTTTCGGTCAGCCCGGAAAACAGGTTATTAGGGGTTTCGAGGGTTTTTTTTTTTTTTTGctaaa                                                                                                                                               >  2:228955/1‑235 (MQ=255)
                                                                                           ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCATCTTCCGGTAAGGTGTGTTTTTGTTCCATTTTTCCGTGCGCCGCGAAAGCCCGGTTTCAGGGTTTTGAGGGGTttttttttttttttgacaaa                                                                                                                                                >  2:377637/1‑233 (MQ=255)
                                                                                           ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                               <  1:562041/239‑1 (MQ=255)
                                                                                            tGGCGCGCGTTTCACTAAAACCATCTTACGTGGGGAGTTGATTTACCACAGTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGGTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac                                                                                                                                              <  2:407283/239‑1 (MQ=255)
                                                                                                  gcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCaa                                                                                                                                                 <  1:752731/230‑1 (MQ=255)
                                                                                                      ttcgcGAAAGCCATCTTACGTGGTTATTTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                                                                                                     <  1:717188/234‑1 (MQ=255)
                                                                                                            aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGTGTGTTTATGTTTCATTTTTCTGTCAGCCGCCAAAACCTGTTTTCAGGGTTTTAGTGGGTTTTTTTTTTTTTTTGCAACAAGGGTGTgggtggttct                                                                                                                              >  2:752731/1‑231 (MQ=255)
                                                                                                                  aTCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGGGGTTTTGTTTAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTCAGTGGTTTTTTTTTTTTTTTTTTAACACCGTTTTGGTTTCCCCTAcccca                                                                                                                        >  1:500929/1‑235 (MQ=255)
                                                                                                                         ggtgtgatgttatTTCACACATTGAACGGGGCTTCCCTGTTGCGCCGAAAGTTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTg                                                                                                                 <  2:31975/236‑1 (MQ=255)
                                                                                                                               atgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCGGTCAGCCGCGAAAGCCTGTTTTCAGGATTTCACGGGGTTTTCTTTATTTTTTTGAACACGGTTGTGGTTCTTCCTACCCGACTTTTTGCttttt                                                                                                           >  1:108902/1‑238 (MQ=255)
                                                                                                                                       tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGGTTTTCTTTATTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGcc                                                                                                    >  1:999791/1‑238 (MQ=255)
                                                                                                                                              ttGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                           <  1:346120/208‑1 (MQ=255)
                                                                                                                                                taccaggctttCCCGGTTGCCCCGAAAGGTCAATTTTTCCTTAAACCTCCGCAGTAATCTGGCCCATGCAATGCCCGTCCTTGCGGCGGGCATTGTCCGGTTAATGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                          <  2:1023516/230‑1 (MQ=255)
                                                                                                                                                taacctgggttCCCTTGTGGGCCGAAAGGTCAATTTTTCCTTAAACATCAGCAGTAATCTGGCCTCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                          <  2:500929/230‑1 (MQ=255)
                                                                                                                                                       gCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTTTGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACCCTACTTTTTGCTTTTCAGTTGCCGCAAAGTAGCCTGCtt                                                                                    >  1:471421/1‑238 (MQ=255)
                                                                                                                                                       gCTTCCCTGTTGCGCAGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCt                                                                                     <  1:899792/237‑1 (MQ=255)
                                                                                                                                                            ccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACCTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTc                                                                              <  1:194568/239‑1 (MQ=255)
                                                                                                                                                                    ccccaaaGGTCTTTTTTTCTTTAAACATCAGCAGTAATCTGGCCCGTGCAATGCCCGTCCGTTCGGCGGCCCTTTCACGGGTGAGGAGGGGTTAGGTTACATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCCACGCTGCTTTCTGCTTTTCAGTTGCCGCAAAGTAGCTTGCTTACGCTCACGCAata                                                                      <  2:504444/235‑1 (MQ=255)
                                                                                                                                                                              cAAATTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGGCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCt                                                             <  1:26583/238‑1 (MQ=255)
                                                                                                                                                                                  ttATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTa                                                        <  1:591378/239‑1 (MQ=255)
                                                                                                                                                                                       cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCa                                                                         <  1:505970/217‑1 (MQ=255)
                                                                                                                                                                                       cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                   <  1:540941/239‑1 (MQ=255)
                                                                                                                                                                                           aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                                                                                                                                                                             >  1:893793/1‑54 (MQ=255)
                                                                                                                                                                                           aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                                                                                                                                                                             <  2:893793/54‑1 (MQ=255)
                                                                                                                                                                                           aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTTTGTTCAATTTTGCAGTCAGCCGCGGAAGCCTGTTTTCAGGATTTCAGTGGTTTTTTTTTATTTTTTGCCACACGGTTTTGGTTCCTCCTACGCTACTTTTTTCTTTTCCGTTTCCCCCAAAGAACCCGCTTACCCTCCCGCCAAATGCCTTTTTTTCTccccccccgg                                               >  2:253071/1‑232 (MQ=255)
                                                                                                                                                                                              aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCt                                                             <  1:290117/222‑1 (MQ=255)
                                                                                                                                                                                                cagcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGGTTTTCTTTTTTTTTTTCAACACGGTTGTGGGTTCTCCTACCCTACTTTCTGCTTTTTAGTTGTCGCAAAGATGCCTGCTTTCTCTCACGCAATATGCTTTTTTTTCTTCCGCACTGGCCt                                            >  1:67718/1‑237 (MQ=255)
                                                                                                                                                                                                      tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                                                                                                                                                                  >  1:365818/1‑53 (MQ=255)
                                                                                                                                                                                                      tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                                                                                                                                                                  <  2:365818/53‑1 (MQ=255)
                                                                                                                                                                                                       aaTCTGGCCCCTGACATTCCCTTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTc                                   <  1:184135/239‑1 (MQ=255)
                                                                                                                                                                                                        aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                                                                                                                                                           >  1:660662/1‑87 (MQ=255)
                                                                                                                                                                                                        aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                                                                                                                                                           <  2:660662/87‑1 (MQ=255)
                                                                                                                                                                                                        aTCTGGCCCCTGAAATGCCCGTCCGTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCa                                                                             <  1:830714/196‑1 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGTGGTTTTGGTCAATTTTGCAGTTAGCCGCGGAAAGCCGTTATCAGGGTTTTAGTGGTTTTTTTTTTTTTTTTTGTACACCGGGTTTGTTTCTCCCACCCGCCTTTATCCTTTTTCGTTGGCCCAAAGTAGGCCGCTTTCGCTCACGCCAAATCTGTTTTTTTTCACCCCACCGGCCTGCTTCCCCCAGGGGt                             >  2:772325/1‑239 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCCATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                          >  1:285115/1‑130 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAac                                                                                                                                                <  1:582755/124‑1 (MQ=255)
                                                                                                                                                                                                             gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                          <  2:285115/130‑1 (MQ=255)
                                                                                                                                                                                                                  tGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTCAGTGGTTTTTCTTTTTTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTCCTTTTTGCTTTTTAGTTTCCGCCAAGTAGCCTGCTTACGCTCACGCCATATTCTTTTTTTTTTACCCCCCTGGCCTCCTTCGCTCAGGGGTGTTgg                        >  2:1018101/1‑238 (MQ=255)
                                                                                                                                                                                                                   gCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                                                                                                                                                                         <  1:806870/93‑1 (MQ=255)
                                                                                                                                                                                                                       tGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAGGTGTGTTTATGTTCAATTTTGCCGTCAGCCGCGAAAGCCTGGTATTAGGGTTTCCGTGGTTTTTTTTTATTTTTTTCAACACGGTTGTgg                                                                                                                                      >  2:582755/1‑124 (MQ=255)
                                                                                                                                                                                                                             ttcTTGCGGCGGTCATTTTCCGGTTAAGGTTTGTTTATGTTCAATTTTGCCGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg             <  1:901329/237‑1 (MQ=255)
                                                                                                                                                                                                                             tCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg             <  1:1018101/239‑1 (MQ=255)
                                                                                                                                                                                                                              ttttGGGGGGGGCTTTCCCGGTTTAGTTGTGTTTTTGTTCAATTTTGCGTGAGGCCGCGAAAGCCTGTTTTCAGGATTTCAGTTGTTTTTCTTTATTTTTTGCAACCCGGTTGGGGTTCCTCTTACGCTACTTTCTGCTTTTCAGTGGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg           <  2:999791/237‑1 (MQ=255)
                                                                                                                                                                                                                              ccTTGCGGCGGGCATTCTCCGGTTAAGGGGTTTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACTGTTGTGGTTCCTCCTACTCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg             <  1:138527/238‑1 (MQ=255)
                                                                                                                                                                                                                               cTTGCGGCGGGCATTTGCCTGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg           <  1:272999/239‑1 (MQ=255)
                                                                                                                                                                                                                               cTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg           <  1:377637/239‑1 (MQ=255)
                                                                                                                                                                                                                                        ggCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTTCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCCCTGGCCTGCTTCGCTCAGGCGTTTTGCGGTTATCGTCCCGCtattattg  >  1:510437/1‑237 (MQ=255)
                                                                                                                                                                                                                                        ggCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTTGCCTGCTTCGCTCAGGCGTTTTTCGGGCATCTTCGCGCcatttttg  >  1:695521/1‑235 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
TTCGTCTTCATTCAGCCGAATAGCAGCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGG  >  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: