Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 540,558 T→G 29.7% intergenic (+50/‑7) 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,5580TG29.7% 128.0 / 16.9 70intergenic (+50/‑7)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base T (12/35);  new base G (20/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 2.80e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCT  >  NC_000913/540333‑540794
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
gcTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTTTGTTCa                                                                                                                                                                                                                                       >  1:407094/1‑239 (MQ=255)
gcTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCAAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGTGTGTTTATGTTCa                                                                                                                                                                                                                                       >  1:416152/1‑239 (MQ=255)
   aTGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGTTTAGGGGTGTTTTTGtttcatt                                                                                                                                                                                                                                    >  2:75333/1‑234 (MQ=255)
   aTGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAAtt                                                                                                                                                                                                                                    <  1:289146/239‑1 (MQ=255)
           aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCTGGGCTTCTTCAGTTATGGGGGGTTTTTGCTCCAGTTTTCAGt                                                                                                                                                                                                                            >  2:205179/1‑238 (MQ=255)
           aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTTCGGTTAAGGGGTGTTTATGTTCAATTTTTGAGTc                                                                                                                                                                                                                           >  1:1124856/1‑239 (MQ=255)
           aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTc                                                                                                                                                                                                                           >  1:72885/1‑239 (MQ=255)
           aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAAAAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGAAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGAGGGAATTTCCCGAGTACTGTGGGTATTTTTTCAAttttttcggg                                                                                                                                                                                                                            >  2:46527/1‑233 (MQ=255)
                tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATCGGCGCGCGAATCACGAAAACCATCTTACGTGGCGATGTGATTTACGACATTGAACAGTGATTCCCTGTTGCGCCGAACGGTCAATTTATCCTTAAAAATCAGCAGTACTCTGGCCCCTGCAATGCCCGTCCTTGCTGCGGGCATTATCCCGGTTAGGGGTGCTTTCTTTCCATTTTGCAGttgacc                                                                                                                                                                                                                       >  2:682860/1‑233 (MQ=255)
                   cGATCTGGAATATCGCCGCAAAGGCAACCCGTATGTTGGCCGTACCATTGGCGCGCTTATCACGAAAACCATCTTACGTCGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCTTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAATCAGCcgcg                                                                                                                                                                                                                    <  2:1124856/238‑1 (MQ=255)
                       cTGTAATATCGCCACAAAGTCAGCCCTTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                                                                                                                                                                               <  2:645385/239‑1 (MQ=255)
                                  ccACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGCCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTCCAATGCCCGTCCTTGCTGCGGGCATTCTTCGGGTAGGGTGGGTTTATGTTCAATTTTTCCGTCAGGAGCGCACGCCCGTGATCAgg                                                                                                                                                                                                   >  2:548019/1‑239 (MQ=255)
                                     cacaGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATACCCGTCCTTGCGGCGGGCATTCTCTGGTTTAGGTGTGTTTTTGTTCACTTTTTCAGTTCGACGCGGAAGACTGTTTTTCGGAtt                                                                                                                                                                                                >  2:342427/4‑239 (MQ=255)
                                        aGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                                                                                                                                                                                            <  1:912099/239‑1 (MQ=255)
                                           cAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAg                                                                                                                                                                                           <  1:208596/237‑1 (MQ=255)
                                                         gCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTTTTTTTACTTTTGCAGTCAGCCGTGTAAACCTTTTTCCCGGATTTCTGTGGGtttttttttttt                                                                                                                                                                           >  1:251888/1‑239 (MQ=255)
                                                           cGTACCAGTGGCGCGCGTATCACGAAAACCAGCTTACGGGGTGATGTGATTTACGGCATTGAACAGGGCTTCCCTGTGGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATGCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAttttt                                                                                                                                                                         <  2:310292/239‑1 (MQ=255)
                                                            gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCGGTTAAGGGGTGTTTATGTTTAATTTTTTAGGCAGCCGCGAAAGCCTGTTATTAGGATTTCGGGGGGTTTTCtttttttttt                                                                                                                                                                        >  1:387234/1‑239 (MQ=255)
                                                              aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGTCTTCTCCCGGTCAGGGGTGCTTTTGTTCTATTTTTGCGTCATCACGGCACGCCTGTTATGCAGATTTCTCGTGATTTTTTTTAtttttt                                                                                                                                                                        >  2:38303/1‑237 (MQ=255)
                                                                cATTGGCGCGCGTATCACGAAAACCATCTTACGTGTGGATGTGATTGACGACAATTAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCaa                                                                                                                                                                    <  2:48268/239‑1 (MQ=255)
                                                                cATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCaa                                                                                                                                                                    <  1:891536/239‑1 (MQ=255)
                                                                  ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                                                  <  2:791298/239‑1 (MQ=255)
                                                                  ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                                                  <  1:784178/239‑1 (MQ=255)
                                                                   tGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac                                                                                                                                                                 <  1:91988/239‑1 (MQ=255)
                                                                    ggCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTTCGGTTAAGGGGTGTTTTTGTTCCATTTTTCAGGCAGCCGCGAAAGCCTGTTATCAGGATTTTAGTGGGTTTTCTTTTTTTTTTTTAacaag                                                                                                                                                                >  1:193984/1‑237 (MQ=255)
                                                                          cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTTCGGTTAAGGTGTGTTTATTTTCAATTTTGCAGTCAGCCGCGGAAACCTGTTATCAGGATTTCAGTGGTTTTTTTTTTTTTTTTGCAACACGGTtgtg                                                                                                                                                          >  1:638226/1‑239 (MQ=255)
                                                                          cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGTGGGTTTTTGTTCAATTTTTCAGTCCGACGGGGAAGACTTGTATCAGGGTTTCCATGGGTTTTTTTTTTTTTTTTTTACACGGTtgtt                                                                                                                                                          >  2:310554/1‑238 (MQ=255)
                                                                                      aCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                                             <  1:26767/239‑1 (MQ=255)
                                                                                      aCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                                             <  1:315070/239‑1 (MQ=255)
                                                                                        cATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGtt                                                                                                                                                                                                         <  1:374497/179‑1 (MQ=255)
                                                                                          tCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTTTGTTCCATTTTGCCGTCCGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTTCAACCCGGTTTTGGTTCCTCCTACGCt                                                                                                                                           >  1:516457/1‑237 (MQ=255)
                                                                                          tCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTAAGGGGTGTTTTTGTTCAATTTTTCAGTCAGCCGCGAAAGACCGTTATCAGGGTTTTTGTGGTTTTTTTTTTTTTTTTTGAACACGGGTGTGGGtcctcccctcgt                                                                                                                                            >  2:912167/1‑231 (MQ=255)
                                                                                             tACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACttt                                                                                                                                      <  1:342427/239‑1 (MQ=255)
                                                                                              aCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                                                  <  1:327277/211‑1 (MQ=255)
                                                                                                    tgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGTGTGTTTATGGTTAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTTCGCTTCTTTTTGCttt                                                                                                                               >  1:43826/1‑239 (MQ=255)
                                                                                                             ttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGcc                                                                                                                      <  1:46527/239‑1 (MQ=255)
                                                                                                              tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCCGGTTTGGGGGGTTTATGTTTAAGTTTGTCAGCAGCCGCGAAAACCCGTTTTCAGGCTTTCATTGGTTTTTTTTTATTTTTTGTCACACAGTTGTGGGTACTCCTACGCTGCTTTCTGTTTTTCATTTGcg                                                                                                                      >  2:800285/1‑237 (MQ=255)
                                                                                                              tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTTAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTCCGTGGGTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTTAGTTGcc                                                                                                                      >  1:918886/1‑238 (MQ=255)
                                                                                                                             ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGGAAACCCGTTATCAGGTGTTCCGTGGGTTTTCTTTATTTTTTGCAACACGGTTTTGGGTCCTACTACCCTACTTTTTGCTTTTCTGTTTCCGCAAAGTAGCCTGCtt                                                                                                      >  1:590311/1‑239 (MQ=255)
                                                                                                                             ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTTCCGTCAGCCCCGAAAGCCTGTTATCAGGGTTTCCGCGGGTTTTCTTTATTTTTTTTCACACCGGTGTGGTTCATCCTACACTACTTTCTTCTTTTCATGTGCCGCAAAATTGCCTCCtt                                                                                                      >  2:1158565/1‑239 (MQ=255)
                                                                                                                             ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCAATCTCCGGGTTAGGGGGGTTTATGTTTAATTTTGCAGTCCGCCGCGGAAGCCTGTTATCAGGATTTCCGTGGTTTTTCTTTATTTTTTTCAACACGGTTGTGGTTCCTCCCTCGCTACTTTTTGCTTTTTAGTTTCCGCACAGttaccctctt                                                                                                      >  2:281537/1‑230 (MQ=255)
                                                                                                                                  cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTCCCCGTTAAGGGGTTTTTATGTTCTATTTTTCAGTTACGCGCGCAAACCTGTTAATCGGGTTTTTGTGGTTTTTTTTTTTTTTTTGCAAAACCGCTGTGGTTCCTCCTCCCGCCATTTTTTCTTTTCAGGTTGCGCGAAGTAGCCtcgtttccc                                                                                                  >  2:1055784/1‑230 (MQ=255)
                                                                                                                                  cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCt                                                                                                <  2:72885/239‑1 (MQ=255)
                                                                                                                                   ccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTc                                                                                               <  1:310554/239‑1 (MQ=255)
                                                                                                                                        tGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCCATTTTGCAGTCAGCCCCGAAAGCCTGTTATTAGGATTTCAGTGGGTTTTCTTTTTTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTTCCGCAAAGGAGCCTGCTTTAGCTCACCCa                                                                                          >  1:764962/1‑239 (MQ=255)
                                                                                                                                           gcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAata                                                                                       <  1:895751/239‑1 (MQ=255)
                                                                                                                                               aaaGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCt                                                                                   >  1:1135113/1‑239 (MQ=255)
                                                                                                                                                     aaaTTTATACTTAAACATCAGCAGTAAGCTGGCTCCTGCATAGCCCGTGCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTGTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTGGCAACACGGTTGTGGTTCCTCCTACGCTGCTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCt                                                                              <  2:638488/237‑1 (MQ=255)
                                                                                                                                                       aTTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAACGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGTTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGc                                                                           <  2:1031917/239‑1 (MQ=255)
                                                                                                                                                         ttATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTa                                                                         <  1:548019/239‑1 (MQ=255)
                                                                                                                                                              cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                    <  1:682860/239‑1 (MQ=255)
                                                                                                                                                                  aaCATCAGCAGTAATCTGGCCCCTGCAATGCACGTCCTTGCGGCGGCCATTCTCCGGTTAAGGAGGGTTTATGTTAAATTTTTCACTCACACGCCAGAGGCTGCTATCAGGAGTTCTGTGGTTTTTTTTTTTTTTTTTGAACACGGTTGTGGTTCCCCCTTCGCTACTTTCTTCTTTTCAGTTGCCCCGAAGTAGACTCCTTAACACCACGCCATATGCTTTTTTTGTTCCCGccccg                                                                 >  2:123972/1‑234 (MQ=255)
                                                                                                                                                                   aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGc                                                               <  1:188188/239‑1 (MQ=255)
                                                                                                                                                                               aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCa                                                                                          <  2:12011/200‑1 (MQ=255)
                                                                                                                                                                               aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCa                                                                                          >  1:12011/1‑200 (MQ=255)
                                                                                                                                                                                    gCCCCTGCAATGCCCGTCCTTGCGGCGGGGCTTCTCCGGGTAAGGGGTGGTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCCGTTATCAGGGTTTTAGGGGTTTTTTTTTATTTTTTGCAACACGGGTGTGGTTCCCCCCACGCTTCTTTTTGCTTTTTTGTTTCCGCAAAATAGCCTGCTTTCGCACACGCCACATGCTTTTTTTTCTACCCCACTGGCCCCCCTCGCTCCGGCGt                                             >  2:915343/1‑239 (MQ=255)
                                                                                                                                                                                      cccTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGc                                                               <  1:125461/220‑1 (MQ=255)
                                                                                                                                                                                       ccTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTTTGTTTACTTTTGCAGTCAGCCGCGGAAGCCTCGTTTTCGGGTTTCAGTGGGTTTTCTTTTTTTTTTTTAACACGGGTGTGGGTCGTCCTACCCTACTTTCTGCTTTTCAGTTTCCCCAAAAAGACCTGCTTAAGCCCCCCGCAAATTCTTTTTTTGCTTACGCACTGGCCTCGCTTGGCCCGGGCtttt                                           >  2:833242/1‑239 (MQ=255)
                                                                                                                                                                                         tGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGGTTATGTTCCATTTTGCAGTCAGCCGCGGAAGCCTGTTATCAGGGTTTCAGGGGGTTTTTTTTATTTTTTTCAACACGGCTGTGGTTACTTCTACGCTACTTTCTGTTTTTCAGTTGACGCAAAGTAGCCTTCTTTCGGTTCCGCCATATGCTTTTCTTTCTCCCGCCCTCGCCTGCCTCGCTCAGGCCTTTTGc                                        >  2:244767/1‑239 (MQ=255)
                                                                                                                                                                                            aaTGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGt                                     <  1:75333/239‑1 (MQ=255)
                                                                                                                                                                                              tGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCa                                   <  1:205179/239‑1 (MQ=255)
                                                                                                                                                                                                     ccTTGGGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                             <  1:1055784/238‑1 (MQ=255)
                                                                                                                                                                                                     ccTTGCGGCGGGCATTCTCCGGTTAAGGGTTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCATGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgc                            <  2:43826/239‑1 (MQ=255)
                                                                                                                                                                                                      cTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                           <  1:833242/239‑1 (MQ=255)
                                                                                                                                                                                                      cTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                           <  1:915343/239‑1 (MQ=255)
                                                                                                                                                                                                                 cATTCTCCGGTTAAGGTGTGTTTATGTTCAAGTTCGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCATTAGTTTTTCTTTAGTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTACTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTCCCTCACTGGCCTGCGTCGCTCAGGCGTTTTGCGGTCAGCGTCGCGCTATTATGGa                 <  2:251888/238‑1 (MQ=255)
                                                                                                                                                                                                                        cGGTTAGTGGGTGTTTATGTTGAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTTGTTTTTCTTTATTTTTTGCAACCCGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAGCGTCGCGCTATTATGGAAGGGCCAg         <  2:193984/239‑1 (MQ=255)
                                                                                                                                                                                                                        cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAg         <  1:38303/239‑1 (MQ=255)
                                                                                                                                                                                                                              aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                            <  1:229574/136‑1 (MQ=255)
                                                                                                                                                                                                                              aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                            >  2:229574/1‑136 (MQ=255)
                                                                                                                                                                                                                                gtgtgtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTgg     >  2:170656/1‑235 (MQ=255)
                                                                                                                                                                                                                                gtgtgtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTCCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTgg     >  1:31735/1‑235 (MQ=255)
                                                                                                                                                                                                                                 tgtgtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCt  <  1:170656/237‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
GCTATGTTCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCT  >  NC_000913/540333‑540794

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