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,5640GA94.5% 199.6 / 0.9 73Q456* (CAA→TAA) yicJputative xyloside transporter YicJ
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (40/29);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTGAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCC  >  NC_000913/3834319‑3834808
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
cTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAgg                                                                                                                                                                                                                                                  >  1:192414/1‑251 (MQ=255)
  gAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGt                                                                                                                                                                                                                                                 <  1:207783/250‑1 (MQ=255)
   aaaaTGCTCAATCCGCACACCGACAATACCGTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTTAATCAAATTGAGGCCAGGGTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGa                                                                                                                                                                                                                                               <  2:61115/251‑1 (MQ=255)
             aTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTg                                                                                                                                                                                                                                    >  2:232838/1‑251 (MQ=255)
                cGCACACCGACAAGACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCTCATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCa                                                                                                                                                                                                                                 <  1:453472/251‑1 (MQ=255)
                 gCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCCTTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGACTCAACCAGTTTCCATCGCTCATTTTCATTGCCGTTCCTACGTTCTGCACTT‑CTTGAGAGGTGAATTGTTG‑CTGGcc                                                                                                                                                                                                                                 >  1:285343/1‑249 (MQ=255)
                    cacCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACgg                                                                                                                                                                                                                             >  2:279564/1‑251 (MQ=255)
                        gacaatcccgtccttgggggAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCATGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAActgagtt                                                                                                                                                                                                                         <  2:627163/235‑7 (MQ=255)
                        gACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACATCGAACACCTGAAGCGGGTGACCCAAATTGAGGCCAGGTTGCATCCACCAGTTTACATCGCCAATTTTCATTGCCGTACCTTAGTTCTGAACTT‑CTTAGGAGGTGCATTTTTG‑CTGGCAACGGCGtt                                                                                                                                                                                                                         >  2:274631/1‑251 (MQ=255)
                          aaaTACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCTTCGCTAATTTGCATTGCCGTTCCTTATTTCTGCACGT‑CTTAAGAGGTGAATTGTTG‑CTGTCAACGGCGTTTa                                                                                                                                                                                                                       <  1:223864/250‑1 (MQ=255)
                               ccTTCCTGGGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTAGCCTGCTGTTCAGCCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTg                                                                                                                                                                                                                  <  1:458156/251‑1 (MQ=255)
                                     tgtgGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTGAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCttatcc                                                                                                                                                                                                              >  2:207465/1‑244 (MQ=255)
                                      gtgGGGAGAAAAAGCGCAACGGAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCTCGCACATCACGTGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCCTTTTCCATCGCTAATTTTCATTGCCGTTCCTTATTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGTCATCGGCGTGTACCCTGAGCCa                                                                                                                                                                                                             <  2:47806/249‑1 (MQ=255)
                                         gtgagaAAAAGCGCAACGTAAATAAAGGCGTATCAATCTGCCAGGAACGTTCACGCACATCACGGGGGGAAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCAAACACCTGAAGCGGGGGAATCAAATTGCGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCTGTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCAAGCTGt                                                                                                                                                                                                        >  2:439420/3‑251 (MQ=255)
                                               aaaaGCGCAACGTAAATAAAGGCGTCTCAAGATGCCAGGTACGTTCACGCACATCACGGGAAGCAGCATAGACCACCATTTCATTATCCTGCTCACCTACCTAGAACACCTGATGCGGGTGAATCAAAATAAGGCCAGGTTGAATCAACCAGTTTCCATCGCTACTTTTCATTGCCATTCCTTAGATCAGCACTT‑CTAAAGAGGTGAATTGTTG‑CTGGCAACGGAGTTTAccccgaacacactgttccat                                                                                                                                                                                                   >  1:316903/1‑233 (MQ=255)
                                                  aGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGACCACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCCGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACg                                                                                                                                                                                               >  2:157800/1‑251 (MQ=255)
                                                            aaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGTTTTTc                                                                                                                                                                                         >  2:19179/1‑247 (MQ=255)
                                                                           aaGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCgtgaggg                                                                                                                                                                      >  1:550743/1‑249 (MQ=255)
                                                                            aGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCgtgagtg                                                                                                                                                                      >  1:439077/1‑250 (MQ=255)
                                                                             gcagtcagttACGTTCACTCACATCACTGGGGGACGCATAGACCGCCATGTCATTATCCTGCTGTGTAATCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAAGCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGTTGAATTGTTG‑CTTGCTACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGGGTCgtgagtgaa                                                                                                                                                                    <  2:237428/242‑2 (MQ=255)
                                                                              cgtccagtgaCGTTCACGCACATCACGGGTGGCAGCATAGACCACCACTTCATTATCCTGCTGTTGTACCTCGAACACCTGAAGCGGGTGAATCTAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTGAAGAGGTGAATTGTTG‑CTGGCATCGGTGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATCGTGCGTCGTgagtgagt                                                                                                                                                                   <  2:1058/242‑1 (MQ=255)
                                                                              cTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCAGAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTaagtgagt                                                                                                                                                                   >  2:130056/1‑251 (MQ=255)
                                                                                   aggttcgtccacgctcgtcccGGGGGTCAGCATAGACCACCATTTCATTAGCCTGCTGTACTACCTCGAACACCGTAAGAGGGTGAATCAAATTTAGGATAGGTTGAATCAACCTGTTGCCATCTCTAATTTTCATTGCCGTGCCTTTGTTCTGTACTT‑CTTAATAGGTGAATTGTTG‑CTGGCAACGGCGTTTACACTGTGCCGGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGATtagtag                                                                                                                                                              <  2:534026/231‑1 (MQ=255)
                                                                                   aGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGagtagt                                                                                                                                                                >  1:362305/1‑249 (MQ=255)
                                                                                   aGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGtagtag                                                                                                                                                              <  1:333661/251‑1 (MQ=255)
                                                                                   aGGTAAGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTCCCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGagtagtgg                                                                                                                                                              >  2:148098/1‑249 (MQ=255)
                                                                                          ttCACGCCCACCACTGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTGCAACCTCGAACACCTGAAGCTGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTATCAGATTGTGCGTCGTTAGTGAGTAGTAGCGTTTa                                                                                                                                                        <  2:368283/250‑1 (MQ=255)
                                                                                           tCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAg                                                                                                                                                       <  1:151565/250‑1 (MQ=255)
                                                                                            cACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCg                                                                                                                                                     <  2:362305/251‑1 (MQ=255)
                                                                                                   aTCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTGAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATcg                                                                                                                                              <  1:207465/251‑1 (MQ=255)
                                                                                                     cACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATcgc                                                                                                                                             <  1:279564/250‑1 (MQ=255)
                                                                                                        ggggggCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTc                                                                                                                                         >  1:181414/1‑251 (MQ=255)
                                                                                                          ggggCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAg                                                                                                                                       <  1:98614/251‑1 (MQ=255)
                                                                                                            ggCAGCATAGACCACCATGTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCATTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCTCGCTCAGCa                                                                                                                                     <  1:471065/251‑1 (MQ=255)
                                                                                                                              ttCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGa                                                                                                                    <  2:181414/250‑1 (MQ=255)
                                                                                                                                ccttATCCTACTGTTCAACCTCGAACACCTGAAGCTGGTGAATCAAAATGAGGCCATGTTGAATCCACCAGTTTCCATCGCTAATTTTCATTTCCTTTCCTTAGTTCTGCACTT‑CGTAAGATGTGAATTGTTT‑CTGGCAACTGCGTTTACCCTGAGCCAGCGGTTCCATAACGGTTTTCATATTGTGCGTCGTGAGTTAGTAGTAGCTTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGcttgaaca                                                                                                                 <  2:445184/249‑8 (MQ=255)
                                                                                                                                cagtATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAAc                                                                                                                  <  1:373294/247‑1 (MQ=255)
                                                                                                                                cATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCCGCAAATAACAGATCGCCGGAAc                                                                                                                  >  1:577742/1‑250 (MQ=255)
                                                                                                                                 aTTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGCGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACa                                                                                                                 >  1:169106/1‑250 (MQ=255)
                                                                                                                                      ccTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCggg                                                                                                           >  1:119725/1‑249 (MQ=255)
                                                                                                                                         gctgTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGTTAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGGGAAt                                                                                                        >  2:15422/1‑251 (MQ=255)
                                                                                                                                         gctgTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAAt                                                                                                        <  1:235953/251‑1 (MQ=255)
                                                                                                                                              tCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑ATTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATATCGGGTTTCAGATTGTGAGTCGTGAGTGTGTAGTAGCGGTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAgcgc                                                                                                   >  2:500658/1‑251 (MQ=255)
                                                                                                                                                     cGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCCGGTTCCAGTGCTAATTTTCATTGCCGGTCCTTAGTTCTGCACTT‑CGTGAGAGGTGAATTGTTG‑CTGGCAACGGCGGTTTCCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGGCGGGAGTGAGTAGTAGCGTATAGCGATAATCGCTGTCAGCAAAGAACTGACCGCCGGAACAATCGTGAATAGCGCaataat                                                                                             <  2:285343/250‑1 (MQ=255)
                                                                                                                                                     cGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCaataat                                                                                             <  1:183322/250‑1 (MQ=255)
                                                                                                                                                      gAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatg                                                                                            >  1:636106/1‑250 (MQ=255)
                                                                                                                                                      gAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatg                                                                                            >  1:192995/1‑250 (MQ=255)
                                                                                                                                                       aaCACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatgat                                                                                          <  1:10646/251‑1 (MQ=255)
                                                                                                                                                        acacCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatgatg                                                                                         <  1:157800/251‑1 (MQ=255)
                                                                                                                                                            cTGAAGCGGGGGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CAGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGAGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCTCCGGAACAATCGTGAATAGAGCAATAATGATGCTa                                                                                      >  2:255582/1‑250 (MQ=255)
                                                                                                                                                                tgCGGGTGAATCAAATTGCGGCCAGGTTGAATCACCCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGGGAATTGCTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTAAATAGCGCCATAATGATGCTAATCGt                                                                                 >  2:155303/2‑251 (MQ=255)
                                                                                                                                                                aGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATc                                                                                                              <  2:258891/222‑1 (MQ=255)
                                                                                                                                                                aGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATc                                                                                                              >  1:258891/1‑222 (MQ=255)
                                                                                                                                                                   gggTGAATCCAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGgc                                                                              >  2:31384/1‑251 (MQ=255)
                                                                                                                                                                    ggTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGACGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTCGCGTTTAGCGATCATAGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGCATAGCGCAATAATGATGCTAATCGTGgcg                                                                             >  1:642543/1‑251 (MQ=255)
                                                                                                                                                                      tGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCt                                                                           >  1:348937/1‑251 (MQ=255)
                                                                                                                                                                          tCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCa                                                                                                                                                                                                    <  2:529734/126‑1 (MQ=255)
                                                                                                                                                                          tCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCa                                                                                                                                                                                                    >  1:529734/1‑126 (MQ=255)
                                                                                                                                                                                 gAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCtt                                                                >  1:340712/1‑251 (MQ=255)
                                                                                                                                                                                        ggTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTTGCGCTGTTCTGCGCTTTTTCTGc                                                         >  2:115422/1‑251 (MQ=255)
                                                                                                                                                                                         gTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGcc                                                        <  1:19179/251‑1 (MQ=255)
                                                                                                                                                                                             aaTCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGcat                                                    >  1:224855/1‑251 (MQ=255)
                                                                                                                                                                                                     aGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCaat                                                                                                <  1:260552/199‑1 (MQ=255)
                                                                                                                                                                                                     aGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCaat                                                                                                >  2:260552/1‑199 (MQ=255)
                                                                                                                                                                                                      gTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGc                                           >  1:663695/1‑251 (MQ=255)
                                                                                                                                                                                                      gTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTCCCGCATCATATCCGc                                           >  2:540684/1‑251 (MQ=255)
                                                                                                                                                                                                      gTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGACGGAACAATCGTGAATAGCGCACTAATGATGCAAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGc                                           >  2:535272/1‑251 (MQ=255)
                                                                                                                                                                                                                 tAATTTTCATTGCCGTAACTTAGTACTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAGCGGCGTTTACCCTGAGCCAGATGTTCCATAAGGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGAACGCAGGAAAACTCGTGAATAGCGCAATAATGATGCTAAACGTGGCGATGTTCTGCGCTTGTTCTGCCGCATCATATCCGCCATCAGACAGc                                >  2:268718/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCCGTTCCTTAGTCCT‑AACTTTCTTAAGAGGT‑AAATGTTGGCGGGCAACGGTGTGTCCCTTGGGCCTGCTGTCCCATAACGTTTTTCGGTTTGTGCGTCGTGAGTGAGTATTAGCTTTTAGCGATAAGCTCGTTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGAAATAATGATGCTAATCGTGGCGCTGTTCTGCCCTTTTTCTGCCGCATCATTTCCGCCATAAGCCAGCATCCAGCCGa                      <  2:636106/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aGTTCTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGcgccg           >  1:651647/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gtCTGCACTT‑CTTAAGAGGTGACTTGTTG‑CTGGCAAAGGCGTTTACACTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGAAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGACGCTGTTCTGCGCTTTTTCTGCCGAATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCgccgcc         >  1:25701/2‑251 (MQ=255)
                                                                                                                                                                                                                                          cTGCACTT‑CTTAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGccgccg        >  1:337471/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ctt‑cttAAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCAGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGcc  >  2:620114/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTT‑CTTGAGAGGTGAATTGTTG‑CTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCC  >  NC_000913/3834319‑3834808

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