Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,834,564 G→A Q456* (CAA→TAA)  yicJ ← putative xyloside transporter YicJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,834,5640GA89.8% 120.8 / 2.5 52Q456* (CAA→TAA) yicJputative xyloside transporter YicJ
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (28/16);  total (35/17)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGG  >  NC_000913/3834322‑3834806
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aaaaTGCTCAATTCGCACCCCGACAATTGCTTCCTGTGGAGAGAAAAAGCGGAAGGTAAGTAAAGGCGTATCAAGCTGCCAGGTACGTTCATGCACAGCACGGGGTGCAGTATAGACCTCCCTTTCATTATCCTTCTGTTCAGCCTCTAACACCTTAAGGGGTTGAATCATATTGAGGCCAGTTTGAATCAACCTGTTTCCAGCGCTAATCTTCATTGGCGTTCCTTAGTTCTGATCTGCTTAAGAGGTGa                                                                                                                                                                                                                                             <  2:23007/251‑1 (MQ=255)
aaaaTGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAACTAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACCTCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGACCACCTGACGCGGGAGAATCAAATTGAGGCCAGGTTGAATAAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTATGCACTTCTTAAGAGGTgc                                                                                                                                                                                                                                             >  2:247750/1‑250 (MQ=255)
   attcttaaTCCGCACAACGACAATACCTGCCATGGGTGAGAAAAAGAGCAAAGTAAAAAAAGGCGTATCAAGCTGCCATGTACGTTCCCGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTCCAACCTCAAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTGAAGAGGTGAAtt                                                                                                                                                                                                                                          <  2:100071/245‑1 (MQ=255)
         aaTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCCTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCt                                                                                                                                                                                                                                    <  1:8904/251‑1 (MQ=255)
                                        ggAGAAAAAGCGCACCGTAAATAAAGGCGTATCAAGCTGCCAGGAACGTTCAAGCACATCACGGGGGGAAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGCCGCATCATATTGAGCCCAGGTTGATTCAACCAGTTTCAATAGCTAATTTTCATTGCCGTTCCTTCGTTATGCCCTTATTTAGCGGTGAATTTATGCTGGAAACGGTGTTTAACCTCCGCCAGCGGtt                                                                                                                                                                                                      >  2:336008/1‑251 (MQ=255)
                                                  cgcAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTATGCACTTCTTCAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAAAGGt                                                                                                                                                                                            >  2:314259/1‑251 (MQ=255)
                                                        gTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACCCCTGAAGCGGGTGAATCAACTTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTc                                                                                                                                                                                        >  1:217881/1‑249 (MQ=255)
                                                              aaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATtgtg                                                                                                                                                                                <  1:61658/251‑1 (MQ=255)
                                                               aaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGc                                                                                                                                                                               >  1:277951/1‑251 (MQ=255)
                                                                aGGCGTATCAATCGGTCAGGTACGTTCACGCACAACACTGGGGCCATCATATACCACCATTTAATTATCCTGCTGTTCCAACGCGAACACCTGAAGCGGTTGAATCAAATTGAGGACAGGTTGAAGCAACCAGTTTCCATCTCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGGCAGCTGTTCCATAACGGTTGTCAGATTGTGcg                                                                                                                                                                              <  2:144905/251‑1 (MQ=255)
                                                                 ggCGTATCAAGCTGCCAGGTACGTTCACGCCCATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGCATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGTATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCAATAACGGTTTTCAGATTATGctt                                                                                                                                                                             >  2:332984/1‑249 (MQ=255)
                                                                      aTCAAGCTGCCATGTACGGTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTAGAACACCTGAAGCGGGTGAATCAAATTGAGGCCATGTTGAATCAAACAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGAACTTCTTAAGAGATGAATTGTTGCTGGCCACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGAGCGTCgtg                                                                                                                                                                         >  2:298161/1‑250 (MQ=255)
                                                                        cAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCgtgagg                                                                                                                                                                      >  1:25389/1‑250 (MQ=255)
                                                                                 aGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGagtagt                                                                                                                                                              >  1:241356/1‑249 (MQ=255)
                                                                                  ggTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGc                                                                                                                                                           >  1:374997/1‑251 (MQ=255)
                                                                                                 aTCACGGGGGGCAGCATAGACCACCATTTCATTCTCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAACTTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATcg                                                                                                                                            >  1:137347/1‑251 (MQ=255)
                                                                                                          ggCAGCATAGACCACCATTTCATGAGCCTGCTGTTCAACCTCGAACCCCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAAGCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCa                                                                                                                                   <  2:69557/251‑1 (MQ=255)
                                                                                                                aTAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAAc                                                                                                                             >  1:98118/1‑251 (MQ=255)
                                                                                                                         caTTTCATTATACTGATGTTCAACCTCAAACAACTACAGCGGGTGAATCAAATTATGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCATCAGTTCTGCACTTCTTACAAAGATAATTGTTGATGGCACCGGAGTTTACCCTGAGCCAGCAGTTCCAAAAAGGTTTTCCGATTGTGAGTAGAGAGTGGTGAGTAGCATTTAACGATAATAGCGCTCAGACAATAACCGATCGc                                                                                                                      >  2:34874/1‑249 (MQ=255)
                                                                                                                             tCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGGTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGaa                                                                                                                 <  1:336008/250‑1 (MQ=255)
                                                                                                                               aTTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGGGTGATGAATGGCGTTTAGCGTTAATCGGGCCCCGCAAATAACAAATCGCCGGaaaa                                                                                                               >  2:1144/1‑248 (MQ=255)
                                                                                                                                        ctgTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTCAATCAACCAGTTTCCATCGCTAATTTTCATTGCAGTTCCTTAGTTCTGCGCTTCTTAAGAGGTGAATTGTTGCTGGCAACGCCGTTTACCCTGAGCCAGCTGTTACATAACGGTTTTCAGCTTGTGCGTCGTGAGTGAGTAGTAACGTTTAGCGATAATCGCGCTAAGCAAATAACAGATCGCCGGAACAATCGTGAAt                                                                                                      >  2:35995/1‑250 (MQ=255)
                                                                                                                                            tCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTCGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATACCAGATCGCCGGAACAATCGTGAATAgcgc                                                                                                 >  1:17116/1‑251 (MQ=255)
                                                                                                                                                ccTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCaat                                                                                              >  1:203824/1‑250 (MQ=255)
                                                                                                                                                   cGAACACCTGAAGCGGGTTAAACAAATTTAGGCCAGGTTGAATCAACCTGTTTCCATCGCTAATTTACATTGCCGTTCCTTAGTTCTGCACTTCTTAATAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCTCaataa                                                                                            <  2:145665/249‑1 (MQ=255)
                                                                                                                                                         ccTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTa                                                                                    <  1:142286/251‑1 (MQ=255)
                                                                                                                                                            taaGCGGGAGAATCAAATTGAGGCCAGGTTGACTCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTAGTGAGTGAGTAGTAGCGTTTAGCGATACTCGCGCTCAGCAAATAACAGATCGACGGAACAATCGTGAATAGCGCAATAATGATGCTAAt                                                                                  >  2:22466/2‑250 (MQ=255)
                                                                                                                                                            gAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTTCGTCGTGATTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTaa                                                                                   <  1:273343/249‑1 (MQ=255)
                                                                                                                                                                 gggTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGgc                                                                            >  1:76415/1‑251 (MQ=255)
                                                                                                                                                                   gTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAg                                                                                                                           >  1:117284/1‑202 (MQ=255)
                                                                                                                                                                     gAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTg                                                                        >  2:230136/1‑251 (MQ=255)
                                                                                                                                                                         cAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCt                                                                    >  1:302322/1‑251 (MQ=255)
                                                                                                                                                                          aaaTTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTg                                                                   <  1:294247/251‑1 (MQ=255)
                                                                                                                                                                               gAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGAGGCGCTGTTCTGCGCtt                                                              >  2:249657/1‑251 (MQ=255)
                                                                                                                                                                                aGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCttt                                                             >  1:38655/1‑251 (MQ=255)
                                                                                                                                                                                     tggTTGAATCAACCAGTTTCCATCGCTATTTTTCATTGCCGTTCCTTAATTCTACACTTCTGAGGAGGTGAATTGTTGCTGGTAACGGCGTTTACCCTGAGCCAGCTGTTCCATAAAGAATTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGATTAGCGATAATCGCTCTCAGCAAATAACAAATCGACGGAACAATCGGGAATCGCGCAATAATGATACTAATCGAGGCCCTGTTCCGCGCTTTTTCTg                                                        >  2:298716/2‑251 (MQ=255)
                                                                                                                                                                                               aaCCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCata                                              <  1:314259/251‑1 (MQ=255)
                                                                                                                                                                                                   aGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGCTCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCATGGCGCTGTTCAGCGCTTTTTCTACCGCATCATATCCg                                          >  2:244316/1‑251 (MQ=255)
                                                                                                                                                                                                             gCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCTGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGGGAGCAGTAGCGTTTCGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGAGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATccacactaagcc                                 >  2:361557/1‑240 (MQ=255)
                                                                                                                                                                                                              cTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGCGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTAAGATTGTGCGTCGTGAGTGAGAAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGGGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTTCCCCCACATCACCGCCATAAGCCa                                >  2:203172/1‑250 (MQ=255)
                                                                                                                                                                                                                  ttttCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGTGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGAGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATc                           >  1:203598/1‑251 (MQ=255)
                                                                                                                                                                                                                         ttccGTTCCTTAGTTTTGCACTTCGTAAGAGGTGAATTTTTGCTGTCACCGGCGTTTACCTTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTGGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGAAAATAACTGATCCCGGGAACAATCGTGAATAGCGTAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCTCATCTTATCCGCCATAAGCCAGCATCCAGCCg                     <  1:189927/248‑1 (MQ=255)
                                                                                                                                                                                                                         tGCCGTTCTTTAGTTCGGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGATTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCATCATCCAGCCg                     <  1:58807/250‑1 (MQ=255)
                                                                                                                                                                                                                             gTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGc                                                                                                                      <  1:203134/149‑1 (MQ=255)
                                                                                                                                                                                                                             gTTCATTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGc                                                                                                                      >  2:203134/1‑149 (MQ=255)
                                                                                                                                                                                                                                ccTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGAGCTGTTCTGCGCTTTTTCTGCCTCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAgc             >  2:110573/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tAGTTCTGCTCTTCTTAAGAGGTGAATTGTTGCTGTCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGGTCGTCGTGAGTGAGTAGTAGCGTGTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGGGGCGCTGTTCTGCGCttttt                                                           <  2:204358/202‑1 (MQ=255)
                                                                                                                                                                                                                                   tAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCttttt                                                           >  1:204358/1‑202 (MQ=255)
                                                                                                                                                                                                                                   tAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCa                                       <  1:228776/222‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGcgccgc        >  1:214005/1‑248 (MQ=255)
                                                                                                                                                                                                                                         tGCACTTCTTAAGAGGTGAATTGTTGCTGGCAAAGGCGTTTACACTGAGCCAGCTGTTACATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGCTGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGTATCATATCCGCCATACGCCAGCATCCAGCCGATAAAAGCGccgccgc     >  1:271931/1‑249 (MQ=255)
                                                                                                                                                                                                                                           cACTTCTTGAGAGGTGAATTGTTGCTGGCCACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGCGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCATGGCGCTGTTCCGCGCTTTTTCTGCCGCATCATATCAGCCATACGCCAGCATACCGCCGATAAGCGCGCCGCCGAAgg  >  2:199362/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
AAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGG  >  NC_000913/3834322‑3834806

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: