Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 540,558 T→G 24.6% 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,5580TG24.6% 81.0 / 5.6 52intergenic (+50/‑7)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base T (7/30);  new base G (12/3);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 6.46e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGA  >  NC_000913/540340‑540779
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                
tCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGAATTCCCCGTTAAGGGGTGgtttttt                                                                                                                                                                                                                       >  2:230111/1‑225 (MQ=255)
    aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGCTTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGGCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGTTGGTTTTGTTCa                                                                                                                                                                                                                    >  1:687399/1‑228 (MQ=255)
    aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGATCAATTTATCCTTAAACATCAGCAGTAATCTGGCACCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTATGGTGTGTTTATGTTCaa                                                                                                                                                                                                                   >  1:368607/1‑229 (MQ=255)
         tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGACAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGTCCCCTGCAATGCCCGTCCTTGCGGCGGGAATTCTCCCCGTAAGGGGTGGTTCTATTCagacttg                                                                                                                                                                                                              >  1:127040/1‑223 (MQ=255)
                cTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGTGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                                                                                                                                                                       <  1:283155/229‑1 (MQ=255)
                          cgccaaaaGGCCGCCCGTATGTTGGCCGTCGCATGGGCGCGCGTATAACGAAAACCATCTTACGGTGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                                                                                                                                                             <  2:687398/224‑1 (MQ=255)
                          cccaccaATTCAGCCCGTATGTTGCCCGTACCAGTGGCGGGCGTATCACGAAAACCATCTTACGTGGTGATGGGATTTACTACATTGAACAGGGCGTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                                                                                                                                                             <  1:230308/228‑1 (MQ=255)
                          cccACAAATGCCGCCCGTTTTGTGGCCGTACCATTGGCGACCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                                                                                                                                                             <  2:653894/228‑1 (MQ=255)
                             acaAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTg                                                                                                                                                                                          <  1:234286/229‑1 (MQ=255)
                                    cAGCCCGTATGTTGGCCATTCCATTGGCGGGCGTATCCCGAAGTCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAg                                                                                                                                                                                   <  1:213170/229‑1 (MQ=255)
                                     aGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGCGGTGGTTATTTTTAATTTTTCGTGCAGCcgcggagaccgggtttctgg                                                                                                                                                                                  >  1:56887/1‑213 (MQ=255)
                                       cccGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGTGGTTTATTTTACTTTTTAAGTCAGcccgcaagacacgtttctaggat                                                                                                                                                                                >  1:47478/1‑208 (MQ=255)
                                                             ggCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCTGCGGCGCTTCTCCTGTTTAGGTGGGATTTTGTTCAATTTTGCAGTCCGCCCCCCAGGCCTGTTTCCAGGATTTTAGTAGGTtttttttttttt                                                                                                                                                         >  1:435125/1‑229 (MQ=255)
                                                                   cgTATTACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCACCTGCAATGCCCGTCCTTGCGGCGGGAATTCTCCGGTTAAGGGGTGGTTATGTTTATTTTTGCTGTCATCCGCGAAAGCCGCTTTTTAGGATTTCTGTGGTTTTTTTTTTTTTTTTTCaa                                                                                                                                                   >  2:726795/1‑229 (MQ=255)
                                                                         ccGAAATCCATCGTCCGTGTGGATGTGATTTTCGTCATGGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACAGCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCCTGCGGCGGGCATTCTCCGGTTAAGGTTTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGt                                                                                                                                             <  2:127040/228‑1 (MQ=255)
                                                                          cGAAACCCATTTTGCGTGGTGATGTTATTTACGACATTGAACAGGGCTTCCCGGTTGCGCCGAAAGGGCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                            <  2:223781/229‑1 (MQ=255)
                                                                                ccATCTTGCGTGGTGATGTTATTTACGACATTTAACAGGGGTTCCCTGTTGCGGCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATTTGGCCCCTGCAATGCCCGTTCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                                                                                                       <  1:80312/228‑1 (MQ=255)
                                                                                ccATCTTGCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                                                                                                       <  1:191847/228‑1 (MQ=255)
                                                                                   tCTTACGTGGTGATGTGATTTACGACATTGCACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCAGTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCCGTCAGCCCCGCAAGCCTGTTATACAGAGTTTACTGGTTTTTCTTTATTTTTTTGACAATGGTTGTGGTtc                                                                                                                                     >  1:579177/1‑227 (MQ=255)
                                                                                          tggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGTGGTTTTTTTTCATTTTTCCGTCGACGCCCAAAACCCGGTTTCAGGTTTTCTCTTGTGTTTTTTTTTTTTTTTGCACAACGTTTTGGTTTTTCCTatcc                                                                                                                            >  2:78637/1‑226 (MQ=255)
                                                                                                attgtttttcgcaattggacgcggCTTCCCTGGAGCGCCCAAAGATCAACTTATCCTTAACAATCCGCCGTAATCTGGCCCCGGAAATGCCCGTCCGTGCGGGGGGCATTCTCCGGTTAAGGGGTGTTTATGTGCAATTTTTCAGTCAGCCGCGAAAGCGTGTTATCAGGATTTCAGTGGGTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACttt                                                                                                                      <  2:368607/207‑1 (MQ=255)
                                                                                                       tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGTCATTCTTCGGTTAGGGGGTGCTTAGTTTCCATTTTGCAGTGCTCCCCGAAAGCCTGTTATCTGGATTTTTGGGGGTTTTTTTTTTTTTTTGCAAACCGGTTTTGGGtcctccccgcccactttttttttt                                                                                                               >  2:691910/1‑211 (MQ=255)
                                                                                                            cATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGt                                                                                                         <  2:648017/229‑1 (MQ=255)
                                                                                                              ttGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCCCGAAAGCCTGTTATCAGGGTTTCAGTGGTTTTTCTTTTTTTTTTTCAACACGGGTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGttt                                                                                                       >  1:523993/1‑228 (MQ=255)
                                                                                                              ttGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTAAGGGGTGGTTTTTTTCAATTTTTCAGGCCGCCCCGCAAGGCTGTTATCAGGGTTTTATGGGTTTTTTTTTTTTTTTTTGACAACCGGTTTGGTTCCTCCCACCCCACttttttttttttcgttt                                                                                                        >  2:746926/1‑215 (MQ=255)
                                                                                                                      ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAAGTGTGTTTTTGTTCCATTTTGCCGTCAGCCGCGAAAGCCTGTTATCAGGATTTCTGGGGTTTTTTTTTTTTTTTTGCAAAACGGTTGTGGGTCCTCCTACGCTCCTTTCTTTTTTTCAGTTTCCGCcaa                                                                                                >  1:596229/1‑228 (MQ=255)
                                                                                                                          tCCCTGTTGCGCCGAAAGGTCAATTTACCCTTAATCATCAGCAGTAATCTGGCCCCTGCAATGCCCGGCCTTGCGGCTGGCATTCTCCGGTTAAGTTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                           <  2:575972/229‑1 (MQ=255)
                                                                                                                          tCCCTGTTGCCCCGAAAGGTCAATTTATCCTTAAACATCAACAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAGTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                           <  1:124084/229‑1 (MQ=255)
                                                                                                                          tCCCTGTCAGGCCGAAAGTTTACTCTACCCGTAAACATCATCAGTTACCTGGCCCCTGCAATGCCCGTCCTGGCGGGGGGCATTCACCGGTTAAGGTGTGGTTATGTTGAATTTTTCAGTAAGCCGCGAACGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGATTTTCAGTTGCCGCAAAGTCGc                                                                                           <  2:398761/229‑1 (MQ=255)
                                                                                                                           cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATAAGCAGAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTCAGGGGTGTTTATGTTCAATTTTTCAGTCAGCCGCCCGAGCCCTTTTTACGGATTTCACTGGTTTTTTTTTTTTTTTTTCTAACCGGTTGTGGTTCCTCCTACTCTTATTTTTTCTTTTTACTTTCCCCaaaaaaaccg                                                                                         >  2:684771/1‑222 (MQ=255)
                                                                                                                                     ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTGAAGGTGTGTTTATGTTTTCTTTTTCAGTCACGCCCGTAAGCCGTTTATCTCGATTTAAGGGGTTTTTCTTTATTTTTTTCAACAAGGTTGTGGTGCTTCCTACGCTCCTTTTTTCTTTTCATTTTCCCGCACAGAGCCTGCTTTCGcc                                                                                >  2:395513/1‑228 (MQ=255)
                                                                                                                                      cGAAAGGTCAATTTATCCTTAAACATCCGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTc                                                                               <  1:691911/229‑1 (MQ=255)
                                                                                                                                            gTCAATTTATCCTTAAACCTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCaa                                                                         <  1:307815/229‑1 (MQ=255)
                                                                                                                                              cAATTTTTCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCACCGGTTAAGGTGTTTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAata                                                                       <  2:406247/229‑1 (MQ=255)
                                                                                                                                               aaTTTATCGTTAAACAGCAGCAGTAAGCTGGCCCCTGCACTGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAata                                                                       <  1:721403/228‑1 (MQ=255)
                                                                                                                                                  ttATCCTAAAACATTAGCAGTAATTGTGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCt                                                                   <  1:726796/229‑1 (MQ=255)
                                                                                                                                                  tgtgtaaaagagacgggcagTAATCTTGCCCCTGCAATGCCCGTCCTTGGGGCGGGCATTCTCCGTTGAAGGGGTGTTTATGTTCAAGTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGGTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCGTACGCTCACGCAATATGCt                                                                   <  1:986/213‑1 (MQ=255)
                                                                                                                                                          aaaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTTAGGGGTGTTTTTGTTCAATTTTTCCGTCAGCCGCGCAAGCCCTCTGTCAGGCTTTCAGTGGGTTTTCTTTATTTTTTTCAACCCGGTCTTTGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTTCCGCAAAAAAGACTTCTTTCGCGCCCCCCATATGCTtttttttc                                                           >  2:421834/1‑227 (MQ=255)
                                                                                                                                                            catctagcagTAATTTGGCCCCTGCAATGCCCTTCCTTGGGGGGGGCATTCTCCGGTTAAGGTTTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTa                                                         <  2:219482/224‑1 (MQ=255)
                                                                                                                                                                 cgctgTAATCTGGCCCCTGATATGCCCTCCCCTGCGGCTGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                    <  2:102910/225‑1 (MQ=255)
                                                                                                                                                                      ttatcggtccccgtgattGCCCTTTCTTGGGGCTGGGATTTGCCGGTTAAGGGGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGc                                               <  2:579176/212‑1 (MQ=255)
                                                                                                                                                                       aacttGGCCCCTGTAATGCCCGTCCTTGCGGCGGGCATTCCCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGc                                               <  1:699250/224‑1 (MQ=255)
                                                                                                                                                                                                       gggCATTCTCCGGTTAAGGTGTGTTTATGTCAAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg              <  2:596228/229‑1 (MQ=255)
                                                                                                                                                                                                        ggttttCTGCGTTTAAGTTGTTTTTTTTTTCCATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGTTTTTCAGTTGTTTTTCTTTATTGTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgc             <  2:47478/225‑1 (MQ=255)
                                                                                                                                                                                                        ggCATTCTCCGGTTAAGGTGTTTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCGTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg              <  1:185999/228‑1 (MQ=255)
                                                                                                                                                                                                         gaaTTCCCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGTTTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTTCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg            <  1:746927/227‑1 (MQ=255)
                                                                                                                                                                                                         gCATTCCCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTGTCTTTATTTTTTGAAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg            <  1:729639/229‑1 (MQ=255)
                                                                                                                                                                                                         acattcgtcgtgtAAGGTGTTTTTATGCTCAAGTTTGCATGCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCGGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg            <  1:432498/217‑1 (MQ=255)
                                                                                                                                                                                                               tcCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCtatt       <  1:78637/228‑1 (MQ=255)
                                                                                                                                                                                                                   gTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGa  >  1:520247/1‑229 (MQ=255)
                                                                                                                                                                                                                       aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                    >  1:622854/1‑175 (MQ=255)
                                                                                                                                                                                                                       aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCCAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                    <  2:622853/175‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                
TCTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGA  >  NC_000913/540340‑540779

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