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,5640GA96.5% 173.9 / ‑1.9 58Q456* (CAA→TAA) yicJputative xyloside transporter YicJ
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (30/25);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

GAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCC  >  NC_000913/3834321‑3834808
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAgg                                                                                                                                                                                                                                                 <  1:160572/249‑1 (MQ=255)
          ataccGCACACCGACAATACCTTACTGTTGTGAGAAAAAGCGCAACGTAAATTAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACAACACGGGGGGCAGCATAGACAACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTGCCATCGCTAAATTCCATTGCGTTTCCTTAGTTGTGCTCTTCTTAAGACGTGAAtggttgtt                                                                                                                                                                                                                                     <  2:424650/248‑3 (MQ=255)
          aaTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACGTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCt                                                                                                                                                                                                                                     <  2:374665/251‑1 (MQ=255)
          aaTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGATGTGAATTGTTGCt                                                                                                                                                                                                                                     >  1:145879/1‑251 (MQ=255)
                     cGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCAAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGt                                                                                                                                                                                                                          >  1:441811/1‑251 (MQ=255)
                                  cTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGcc                                                                                                                                                                                                              >  2:201737/1‑250 (MQ=255)
                                      ggggAGAAAAAGCGCAACGTAAATAAAGGCGTATAAAGCTGCCAGGTACGTTCCCGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTCATGTTCAACCTCGAACACCTGACGCGGGCGAATCAACTTGAGACCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTAATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGTAGACTTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTg                                                                                                                                                                                                         >  1:467451/1‑251 (MQ=255)
                                      ggggAGAAAAAGCGCAACGTAAATAAAGGCGGATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTg                                                                                                                                                                                                         <  2:21380/251‑1 (MQ=255)
                                           gaAACAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGAAGCATAGACCACCATTTCATAATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCAATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCCCTTCTTCAGAGGTGACTTGGTGCTGGCAACTGCGTTTACCGTGAGCCAGCAGttac                                                                                                                                                                                                     >  2:391594/1‑248 (MQ=255)
                                                          aaataaaGGCTTAACAATCTGCCAGGTACGTTCACGCCCATCACGGGGGTCAGCATAGACCACCATTTCATGTTCCCGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCATGTTGAATCAACCAGTTTCCAGCGCTAATTTTCATTGACGTTCCTTAGTTCTGCCGTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTTTTCCATAACGGTTGTCAGAt                                                                                                                                                                                     <  2:288978/251‑1 (MQ=255)
                                                          aaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTCCCCTGAGCCAGCTGTTCCATAACGGTTTTCAGa                                                                                                                                                                                      <  1:405786/250‑1 (MQ=255)
                                                              aaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATtgtg                                                                                                                                                                                 <  1:357385/251‑1 (MQ=255)
                                                               aaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGc                                                                                                                                                                                >  2:227514/1‑251 (MQ=255)
                                                                aGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGcg                                                                                                                                                                               >  1:203770/1‑251 (MQ=255)
                                                                  gCGTATCAAGCTGCCATGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGcgtc                                                                                                                                                                             <  1:363084/251‑1 (MQ=255)
                                                                         aaGCTGCCAGGTACGTTCACGCACATAACGGGGGGCAGCATAGAACACCATTTCATTATCCTGCTGTTAAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGTGGTGAATTGTTGCTGGTAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCgtgagtg                                                                                                                                                                      >  1:179336/1‑251 (MQ=255)
                                                                                                          ggCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTTAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCa                                                                                                                                     <  1:337326/251‑1 (MQ=255)
                                                                                                          ggCAGAATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGAGAATCAAATTGAGGCCAGGTTGAATCAACCAGATTCCATCGCTAATTTTCATTGCCGTTCCTTAGATCTGCACTTCTTAAGAGGTGCATTGTTGCTGGCAACGGCGTTTACACTGAGCCAGCTGTTCCATAACGGTTATCAGATTGTGCGTCGTGAGTGAGTAGTATCGTTTAGCGATAATCGCGCTCAGc                                                                                                                                      >  1:318506/1‑250 (MQ=255)
                                                                                                               caTAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATa                                                                                                                                 >  2:373454/1‑250 (MQ=255)
                                                                                                                   gACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGa                                                                                                                            <  1:345288/251‑1 (MQ=255)
                                                                                                                     ccaccaTTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATTAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATc                                                                                                                          <  1:433059/251‑1 (MQ=255)
                                                                                                                     ccaccaTTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATc                                                                                                                          <  1:273798/251‑1 (MQ=255)
                                                                                                                        ccaTTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGcc                                                                                                                       >  2:154747/1‑251 (MQ=255)
                                                                                                                                 tATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAAt                                                                                                               <  2:167472/250‑1 (MQ=255)
                                                                                                                                    ccTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGt                                                                                                            <  1:154747/250‑1 (MQ=255)
                                                                                                                                      tgctgTTCAACCTCGAACACCTGACGAGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGACTTGTTGCTGGCAACGTCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGa                                                                                                          >  1:87686/1‑250 (MQ=255)
                                                                                                                                         tgTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAAt                                                                                                        >  1:6398/1‑249 (MQ=255)
                                                                                                                                                    gAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGGGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAACGCAATAatg                                                                                            >  2:91419/1‑250 (MQ=255)
                                                                                                                                                      acacCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAatgatg                                                                                         >  1:413229/1‑251 (MQ=255)
                                                                                                                                                            gAAGCGGGTGAATCAAATTGCGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGTTAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGTTGCTAATc                                                                                   >  2:376946/1‑251 (MQ=255)
                                                                                                                                                              aGCGGGTGAATCAAATTGAGGCCATGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTTCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGt                                                                                 <  1:454438/251‑1 (MQ=255)
                                                                                                                                                              aGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGt                                                                                 >  1:231673/1‑251 (MQ=255)
                                                                                                                                                               gctggttaaTCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTGCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGGAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTTAATAGCGCACTAAGGATGCTAATCGTg                                                                                <  2:274677/244‑1 (MQ=255)
                                                                                                                                                                    tGAATCAAATTGAGGCAAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCt                                                                           >  2:481124/1‑251 (MQ=255)
                                                                                                                                                                      aaTCAAATTGAGGCCAGGTTGAATCAACCATTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGTTTTTCAGATTGTGCGTCTTGAGTGAGTAGTATCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGt                                                                         <  1:397367/251‑1 (MQ=255)
                                                                                                                                                                         cAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCt                                                                      >  2:253687/1‑251 (MQ=255)
                                                                                                                                                                          aaaTTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGAAGCTAATCGTGGCGCTGTTCTg                                                                     >  2:376755/1‑251 (MQ=255)
                                                                                                                                                                                     aGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTg                                                          >  1:449087/1‑251 (MQ=255)
                                                                                                                                                                                       gTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGcc                                                        <  1:99542/251‑1 (MQ=255)
                                                                                                                                                                                         tGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATtgt                                                                                                                                                                                  <  1:401834/127‑1 (MQ=255)
                                                                                                                                                                                         tGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATtgt                                                                                                                                                                                  >  2:401834/1‑127 (MQ=255)
                                                                                                                                                                                         tGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCACTAATGATGCTAATCGt                                                                                 <  2:356662/224‑1 (MQ=255)
                                                                                                                                                                                         tGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGt                                                                                 >  1:356662/1‑224 (MQ=255)
                                                                                                                                                                                                aCCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCatat                                               >  1:387517/1‑251 (MQ=255)
                                                                                                                                                                                                        cctgcGCTAATTTTCATTGCCGTTCCTTATTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGCCGTGAGTGAGTAGTAGCTTTTAGCGATAATCGCGCTCAGAAAATAACAGATCGCCGGAACAATCGTAAATAGCGCAATATTGATGCAAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATa                                       <  1:271105/247‑1 (MQ=255)
                                                                                                                                                                                                        ccATCGTTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGGGGTGATTTGTTGCTGGCAACGGCGTTTGCCCTGAGCCTGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGAGATAATCGCGCTCAGCATATAACAGATCGCCGGAACAATCGTGAAGTGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGATTTTTCTGCCGcatcattgccgcc                                          <  2:467451/248‑8 (MQ=255)
                                                                                                                                                                                                         cATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATa                                       >  2:293909/1‑250 (MQ=255)
                                                                                                                                                                                                            cGCTAATTTTCATTGCCGTTCCTTAGTTCTGCCCTTCTTAAGAGGTGAATTGTTGCTGTCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGTTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGAAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTTTTCTGCTCTTTTTCTGCCGCATCATATCCGCCATAAGc                                    <  1:110982/250‑1 (MQ=255)
                                                                                                                                                                                                                 aTTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCa                               >  1:108816/1‑250 (MQ=255)
                                                                                                                                                                                                                 aTTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGGGAATTGTTGCTGGCAACGGCGTTTACCCGGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGATTGAGTAGGAGCTTTTAGCGATAATCCCGCTAAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAAAAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCa                               <  2:108816/250‑1 (MQ=255)
                                                                                                                                                                                                                   tttCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGcctcat                                                 >  2:198665/1‑230 (MQ=255)
                                                                                                                                                                                                                   tttCATTGCCGTTCCTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGcatcat                                                 <  1:198665/230‑1 (MQ=255)
                                                                                                                                                                                                                                acTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCCACGGCGTTTACCCTGAGCCAGCTGTTCCCTAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTAAGCAAATAACAGATCGCCGGAACAATCGTGACTAGCGCAATAATGATGCTAATCGTGGCCCTGTTCTGCGCTTTTTCTGCCGCATCCTATCCGCACTAAGCCAGCATCCCGCCGATAATAgc               >  2:335529/2‑251 (MQ=255)
                                                                                                                                                                                                                                 cTTAGTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAgcg              >  1:463355/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gTTCTGCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCCATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCCGCATCCAGCCGATAAGAGcgccgc          >  2:445954/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGa       >  2:49178/1‑249 (MQ=255)
                                                                                                                                                                                                                                          gCACTTCTTAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGACATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGa       >  1:221539/1‑249 (MQ=255)
                                                                                                                                                                                                                                             cttcttAAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGcc  >  1:146427/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTTCTGCACTTCTTGAGAGGTGAATTGTTGCTGGCAACGGCGTTTACCCTGAGCCAGCTGTTCCATAACGGTTTTCAGATTGTGCGTCGTGAGTGAGTAGTAGCGTTTAGCGATAATCGCGCTCAGCAAATAACAGATCGCCGGAACAATCGTGAATAGCGCAATAATGATGCTAATCGTGGCGCTGTTCTGCGCTTTTTCTGCCGCATCATATCCGCCATAAGCCAGCATCCAGCCGATAAGAGCGCCGCCGAAGGCC  >  NC_000913/3834321‑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: