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

CTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCC  >  NC_000913/3834319‑3834808
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
cTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAgg                                                                                                                                                                                                                                                  >  1:192414/1‑251 (MQ=255)
  gAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGt                                                                                                                                                                                                                                                 <  1:207783/250‑1 (MQ=255)
   aaaaTGCTCAATCCGCACACCGACAATACCGTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTTAATCAAATTGAGGCCAGGGTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGa                                                                                                                                                                                                                                               <  2:61115/251‑1 (MQ=255)
             aTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTg                                                                                                                                                                                                                                    >  2:232838/1‑251 (MQ=255)
                cGCACACCGACAAGACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCTCATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCa                                                                                                                                                                                                                                 <  1:453472/251‑1 (MQ=255)
                 gCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCCTTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGACTCAACCAGTTTCCATCGCTCATTTTCATTGCCGTTCCTACGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGcc                                                                                                                                                                                                                                 >  1:285343/1‑249 (MQ=255)
                    cacCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACgg                                                                                                                                                                                                                             >  2:279564/1‑251 (MQ=255)
                        gacaatcccgtccttgggggAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCATGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAActgagtt                                                                                                                                                                                                                         <  2:627163/235‑7 (MQ=255)
                        gACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACATCGAACACCTGAAGCGGGTGACCCAAATTGAGGCCAGGTTGCATCCACCAGTTTACATCGCCAATTTTCATTGCCGTACCTTAGTTCTGAACTTCTTAGGAGGTGCATTTTTGCTGGCAACGGCGtt                                                                                                                                                                                                                         >  2:274631/1‑251 (MQ=255)
                          aaaTACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCTTCGCTAATTTGCATTGCCGTTCCTTATTTCTGCACGTCTTAAGAGGTGAATTGTTGCTGTCAACGGCGTTTa                                                                                                                                                                                                                       <  1:223864/250‑1 (MQ=255)
                               ccTTCCTGGGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTAGCCTGCTGTTCAGCCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTg                                                                                                                                                                                                                  <  1:458156/251‑1 (MQ=255)
                                     tgtgGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCttatcc                                                                                                                                                                                                              >  2:207465/1‑244 (MQ=255)
                                      gtgGGGAGAAAAAGCGCAACGGAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCTCGCACATCACGTGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCCTTTTCCATCGCTAATTTTCATTGCCGTTCCTTATTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGTCATCGGCGTGTACCCTGAGCCa                                                                                                                                                                                                             <  2:47806/249‑1 (MQ=255)
                                         gtgagaAAAAGCGCAACGTAAATAAAGGCGTATCAATCTGCCAGGAACGTTCACGCACATCACGGGGGGAAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCAAACACCTGAAGCGGGGGAATCAAATTGCGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCTGTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCAAGCTGt                                                                                                                                                                                                        >  2:439420/3‑251 (MQ=255)
                                               aaaaGCGCAACGTAAATAAAGGCGTCTCAAGATGCCAGGTACGTTCACGCACATCACGGGAAGCAGCATAGACCACCATTTCATTATCCTGCTCACCTACCTAGAACACCTGATGCGGGTGAATCAAAATAAGGCCAGGTTGAATCAACCAGTTTCCATCGCTACTTTTCATTGCCATTCCTTAGATCAGCACTTCTAAAGAGGTGAATTGTTGCTGGCAACGGAGTTTAccccgaacacactgttccat                                                                                                                                                                                                   >  1:316903/1‑233 (MQ=255)
                                                  aGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGACCACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCCGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACg                                                                                                                                                                                               >  2:157800/1‑251 (MQ=255)
                                                            aaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGTTTTTc                                                                                                                                                                                         >  2:19179/1‑247 (MQ=255)
                                                                           aaGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCgtgaggg                                                                                                                                                                      >  1:550743/1‑249 (MQ=255)
                                                                            aGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCgtgagtg                                                                                                                                                                      >  1:439077/1‑250 (MQ=255)
                                                                             gcagtcagttACGTTCACTCACATCACTGGGGGACGCATAGACCGCCATGTCATTATCCTGCTGTGTAATCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAAGCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGTTGAATTGTTGCTTGCTACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGGGTCgtgagtgaa                                                                                                                                                                    <  2:237428/242‑2 (MQ=255)
                                                                              cgtccagtgaCGTTCACGCACATCACGGGTGGCAGCATAGACCACCACTTCATTATCCTGCTGTTGTACCTCGAACACCTGAAGCGGGTGAATCTAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTGAAGAGGTGAATTGTTGCTGGCATCGGTGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATCGTGCGTCGTgagtgagt                                                                                                                                                                   <  2:1058/242‑1 (MQ=255)
                                                                              cTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCAGAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTaagtgagt                                                                                                                                                                   >  2:130056/1‑251 (MQ=255)
                                                                                   aggttcgtccacgctcgtcccGGGGGTCAGCATAGACCACCATTTCATTAGCCTGCTGTACTACCTCGAACACCGTAAGAGGGTGAATCAAATTTAGGATAGGTTGAATCAACCTGTTGCCATCTCTAATTTTCATTGCCGTGCCTTTGTTCTGTACTTCTTAATAGGTGAATTGTTGCTGGCAACGGCGTTTACACTGTGCCGGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGATtagtag                                                                                                                                                              <  2:534026/231‑1 (MQ=255)
                                                                                   aGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGagtagt                                                                                                                                                                >  1:362305/1‑249 (MQ=255)
                                                                                   aGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGtagtag                                                                                                                                                              <  1:333661/251‑1 (MQ=255)
                                                                                   aGGTAAGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTCCCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGagtagtgg                                                                                                                                                              >  2:148098/1‑249 (MQ=255)
                                                                                          ttCACGCCCACCACTGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTGCAACCTCGAACACCTGAAGCTGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTATCAGATTGTGCGTCGTTAGTGAGTAGTAGCGTTTa                                                                                                                                                        <  2:368283/250‑1 (MQ=255)
                                                                                           tCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAg                                                                                                                                                       <  1:151565/250‑1 (MQ=255)
                                                                                            cACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCg                                                                                                                                                     <  2:362305/251‑1 (MQ=255)
                                                                                                   aTCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATcg                                                                                                                                              <  1:207465/251‑1 (MQ=255)
                                                                                                     cACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATcgc                                                                                                                                             <  1:279564/250‑1 (MQ=255)
                                                                                                        ggggggCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTc                                                                                                                                         >  1:181414/1‑251 (MQ=255)
                                                                                                          ggggCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAg                                                                                                                                       <  1:98614/251‑1 (MQ=255)
                                                                                                            ggCAGCATAGACCACCATGTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCATTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCTCGCTCAGCa                                                                                                                                     <  1:471065/251‑1 (MQ=255)
                                                                                                                              ttCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGa                                                                                                                    <  2:181414/250‑1 (MQ=255)
                                                                                                                                ccttATCCTACTGTTCAACCTCGAACACCTGAAGCTGGTGAATCAAAATGAGGCCATGTTGAATCCACCAGTTTCCATCGCTAATTTTCATTTCCTTTCCTTAGTTCTGCACTTCGTAAGATGTGAATTGTTTCTGGCAACTGCGTTTACCCTGAGCCAGCGGTTCCATAACGGTTTTCATATTGTGCGTCGTGAGTTAGTAGTAGCTTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGcttgaaca                                                                                                                 <  2:445184/249‑8 (MQ=255)
                                                                                                                                cagtATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAAc                                                                                                                  <  1:373294/247‑1 (MQ=255)
                                                                                                                                cATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCCGCAAATAACAGATCGCCGGAAc                                                                                                                  >  1:577742/1‑250 (MQ=255)
                                                                                                                                 aTTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGCGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACa                                                                                                                 >  1:169106/1‑250 (MQ=255)
                                                                                                                                      ccTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCggg                                                                                                           >  1:119725/1‑249 (MQ=255)
                                                                                                                                         gctgTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGTTAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGGGAAt                                                                                                        >  2:15422/1‑251 (MQ=255)
                                                                                                                                         gctgTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAAt                                                                                                        <  1:235953/251‑1 (MQ=255)
                                                                                                                                              tCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTATTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATATCGGGTTTCAGATTGTGAGTCGTGAGTGTGTAGTAGCGGTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAgcgc                                                                                                   >  2:500658/1‑251 (MQ=255)
                                                                                                                                                     cGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCCGGTTCCAGTGCTAATTTTCATTGCCGGTCCTTAGTTCTGCACTTCGTGAGAGGTGAATTGTTGCTGGCAACGGCGGTTTCCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGGCGGGAGTGAGTAGTAGCGTATAGCGATAATCGCTGTCAGCAAAGAACTGACCGCCGGAACAATCGTGAATAGCGCaataat                                                                                             <  2:285343/250‑1 (MQ=255)
                                                                                                                                                     cGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCaataat                                                                                             <  1:183322/250‑1 (MQ=255)
                                                                                                                                                      gAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatg                                                                                            >  1:636106/1‑250 (MQ=255)
                                                                                                                                                      gAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatg                                                                                            >  1:192995/1‑250 (MQ=255)
                                                                                                                                                       aaCACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatgat                                                                                          <  1:10646/251‑1 (MQ=255)
                                                                                                                                                        acacCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatgatg                                                                                         <  1:157800/251‑1 (MQ=255)
                                                                                                                                                            cTGAAGCGGGGGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCAGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGAGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCTCCGGAACAATCGTGAATAGAGCAATAATGATGCTa                                                                                      >  2:255582/1‑250 (MQ=255)
                                                                                                                                                                tgCGGGTGAATCAAATTGCGGCCAGGTTGAATCACCCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGGGAATTGCTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTAAATAGCGCCATAATGATGCTAATCGt                                                                                 >  2:155303/2‑251 (MQ=255)
                                                                                                                                                                aGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATc                                                                                                              <  2:258891/222‑1 (MQ=255)
                                                                                                                                                                aGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATc                                                                                                              >  1:258891/1‑222 (MQ=255)
                                                                                                                                                                   gggTGAATCCAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGgc                                                                              >  2:31384/1‑251 (MQ=255)
                                                                                                                                                                    ggTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGACGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTCGCGTTTAGCGATCATAGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGCATAGCGCAATAATGATGCTAATCGTGgcg                                                                             >  1:642543/1‑251 (MQ=255)
                                                                                                                                                                      tGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCt                                                                           >  1:348937/1‑251 (MQ=255)
                                                                                                                                                                          tCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCa                                                                                                                                                                                                    <  2:529734/126‑1 (MQ=255)
                                                                                                                                                                          tCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCa                                                                                                                                                                                                    >  1:529734/1‑126 (MQ=255)
                                                                                                                                                                                 gAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCtt                                                                >  1:340712/1‑251 (MQ=255)
                                                                                                                                                                                        ggTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTTGCGCTGTTCTGCGCTTTTTCTGc                                                         >  2:115422/1‑251 (MQ=255)
                                                                                                                                                                                         gTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGcc                                                        <  1:19179/251‑1 (MQ=255)
                                                                                                                                                                                             aaTCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGcat                                                    >  1:224855/1‑251 (MQ=255)
                                                                                                                                                                                                     aGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCaat                                                                                                <  1:260552/199‑1 (MQ=255)
                                                                                                                                                                                                     aGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCaat                                                                                                >  2:260552/1‑199 (MQ=255)
                                                                                                                                                                                                      gTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGc                                           >  1:663695/1‑251 (MQ=255)
                                                                                                                                                                                                      gTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTCCCGCATCATATCCGc                                           >  2:540684/1‑251 (MQ=255)
                                                                                                                                                                                                      gTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGACGGAACAATCGTGAATAGCGCACTAATGATGCAAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGc                                           >  2:535272/1‑251 (MQ=255)
                                                                                                                                                                                                                 tAATTTTCATTGCCGTAACTTAGTACTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAGCGGCGTTTACCCTGAGCCAGATGTTCCATAAGGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGAACGCAGGAAAACTCGTGAATAGCGCAATAATGATGCTAAACGTGGCGATGTTCTGCGCTTGTTCTGCCGCATCATATCCGCCATCAGACAGc                                >  2:268718/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCCGTTCCTTAGTCCTAACTTTCTTAAGAGGTAAATGTTGGCGGGCAACGGTGTGTCCCTTGGGCCTGCTGTCCCATAACGTTTTTCGGTTTGTGCGTCGTGAGTGAGTATTAGCTTTTAGCGATAAGCTCGTTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGAAATAATGATGCTAATCGTGGCGCTGTTCTGCCCTTTTTCTGCCGCATCATTTCCGCCATAAGCCAGCATCCAGCCGa                      <  2:636106/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGcgccg           >  1:651647/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gtCTGCACTTCTTAAGAGGTGACTTGTTGCTGGCAAAGGCGTTTACACTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGAAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGACGCTGTTCTGCGCTTTTTCTGCCGAATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCgccgcc         >  1:25701/2‑251 (MQ=255)
                                                                                                                                                                                                                                          cTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGccgccg        >  1:337471/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cttcttAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCAGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGcc  >  2:620114/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCC  >  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: