Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,127,981 G→A P14L (CCC→CTC)  glcD ← glycolate dehydrogenase, putative FAD‑linked subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,127,9810GA96.0% 240.9 / NA 74P14L (CCC→CTC) glcDglycolate dehydrogenase, putative FAD‑linked subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/48);  minor base C (3/0);  total (26/48)
Fisher's exact test for biased strand distribution p-value = 4.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACCAGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGGGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGCGTGTTTTGCGAG  >  NC_000913/3127741‑3128227
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
ccaccaGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGcc                                                                                                                                                                                                                                              <  1:176981/251‑1 (MQ=255)
  accaGAAAGGCCGGTGCCGGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGGCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAAGGGACGCGTGCGATACGCGCTCAACCCGTCCCACTCGTAAGGAATGATCTCCTCATCGGTATTCAGGATCTCAAGGCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGAAATCGAGTAAAGCGCCAt                                                                                                                                                                                                                                            <  2:258938/251‑1 (MQ=255)
   ccaGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATc                                                                                                                                                                                                                                           <  2:70584/251‑1 (MQ=255)
   ccaGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATAAGTACCGATGTGCGGTCGACATCGAGTCAAGCGCCATc                                                                                                                                                                                                                                           >  1:572416/1‑251 (MQ=255)
    caGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCa                                                                                                                                                                                                                                          <  2:313626/251‑1 (MQ=255)
    caGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCa                                                                                                                                                                                                                                          <  2:439082/251‑1 (MQ=255)
        aaGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTGAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGAc                                                                                                                                                                                                                                      <  1:570934/251‑1 (MQ=255)
         aGCCCGGTGCCTGCCCCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGGGCGGGCGACATCGAGTAAAGCGCCATCAAGACg                                                                                                                                                                                                                                     <  2:299672/251‑1 (MQ=255)
           cccGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTAACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGct                                                                                                                                                                                                                                   >  2:80114/1‑251 (MQ=255)
            ccGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGCCAGACAGAATCGCTGTCACCTGTTCCATTTGCTTCGGCAGCACAACCCGTAATGGACGCGTGCGCTACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCCAGTCCAGGGCCATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGctc                                                                                                                                                                                                                                  >  1:580262/1‑251 (MQ=255)
               gTGCCTGCCCCACGGGTCACCACCGGTACACGCATGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCCCCCGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTc                                                                                                                                                                                                                               <  2:137103/251‑1 (MQ=255)
               gTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTc                                                                                                                                                                                                                               <  1:78087/251‑1 (MQ=255)
                 gCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCg                                                                                                                                                                                                                              <  1:320885/250‑1 (MQ=255)
                   cTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTAc                                                                                                                                                                                                                           >  1:495185/1‑251 (MQ=255)
                            gggTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGAGCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATc                                                                                                                                                                                                                                           <  1:539480/226‑1 (MQ=255)
                            gggTCACCACAGGTACACGCAGGCGATGGAAGAAAGCCAGAATCGCTGTCCCCTGTTCCATTTGCTTAGGCAGAACAACCAGAAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCCGGATCTCCAGTCAAGGGACATGCTCACGCAGTGACATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCAATc                                                                                                                                                                                                                                           >  2:539480/1‑226 (MQ=255)
                                 accaccGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATg                                                                                                                                                                                                               >  1:9280/1‑249 (MQ=255)
                                    accGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAg                                                                                                                                                                                                          <  1:183452/251‑1 (MQ=255)
                                      cGGTACACGCAGGCGATGGCAGCCAGCCAGAATCGCTGACACCTGTTCAATTTGCTTAGGCAGATCACCCCGTAATGGACGCGTGCGACACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTAATCGGTATGCAGGATCTCACGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGCGTAAAGCGCCATCAAGCAGCTCTTCGTACAAGATGCTCATGCGTAGGc                                                                                                                                                                                                        >  2:360941/1‑251 (MQ=255)
                                       ggTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCGCATGAGTAGGCt                                                                                                                                                                                                       <  2:66018/251‑1 (MQ=255)
                                          acacGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGCTACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGTTCGACATCGCGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCg                                                                                                                                                                                                    >  2:209061/1‑251 (MQ=255)
                                          acacGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCg                                                                                                                                                                                                    >  1:416416/1‑251 (MQ=255)
                                           cacGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGc                                                                                                                                                                                                   <  2:201884/251‑1 (MQ=255)
                                             cGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTATGCTTCGCtt                                                                                                                                                                                                 <  2:108004/251‑1 (MQ=255)
                                                       ggCAGACAGCCAGAATCGCTGTCACCTGTTCCAGTTGCTGAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACTCTCGTAAGGAATGATCTCCTCATCGTTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCATTTCCGATGTGCGGTCGACATCGAGTAAAGCGCCATCTAGTCGCTCTTCGTACAAGATGCGCATGAGGAGGGTTCGCTTtgttgtgttg                                                                                                                                                                                       <  2:300643/251‑1 (MQ=255)
                                                           gACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGAGCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTGGTGTtgtgtg                                                                                                                                                                                   <  2:395154/251‑1 (MQ=255)
                                                              aGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCa                                                                                                                                                                                <  1:78520/251‑1 (MQ=255)
                                                                ccAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCAGCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGc                                                                                                                                                                              <  2:382744/251‑1 (MQ=255)
                                                                    aagcGCTTTCACCTGTTCAATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCTATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACTCAGTGCCATCAGTTCCGATGTGCGGTCGACATTGAGTAAAGCTCCATCAAGACGCTCTTCTTGCAAGATGCTCATGAGTAGGCTGCGCTTTGTTGTGTTGTGTTGCAGCTGAt                                                                                                                                                                          <  2:465455/248‑1 (MQ=255)
                                                                       cGCTGTCACCTGTTCCAGTTGATTATGCAGAACAACCAGGAAGTGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGAGCTCCTCATCGGTAGGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTCCCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTTTTGTGTGGCAGCTGAtttt                                                                                                                                                                       <  2:330938/251‑1 (MQ=255)
                                                                           gTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTgc                                                                                                                                                                    <  2:150237/250‑1 (MQ=255)
                                                                           gTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCAGCGGTATGCATGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCTAGACGCTCTTCGTACAAGATGCTCATGTGTAGGCTGCGCTTTGGTGTTTTGTGTGGCGGCTGAGTTTTgc                                                                                                                                                                    <  2:59993/250‑1 (MQ=255)
                                                                                    tCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCtt                                                                                                                                                           <  1:80114/250‑1 (MQ=255)
                                                                                      cATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGAGCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCGTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCtg                                                                                                                                                        <  1:378133/251‑1 (MQ=255)
                                                                                        tttGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCtgtg                                                                                                                                                      >  1:433199/1‑251 (MQ=255)
                                                                                                   aGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGttatt                                                                                                                                            >  1:283681/1‑250 (MQ=255)
                                                                                                     aacaacCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGttatta                                                                                                                                           >  2:259913/1‑249 (MQ=255)
                                                                                                     aacaacCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAg                                                                                                                                         <  2:325183/251‑1 (MQ=255)
                                                                                                       caacCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCg                                                                                                                                       >  2:256669/1‑251 (MQ=255)
                                                                                                       caacCAGTAATGGACGCGTGCGATACGCGCGCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACAGCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCg                                                                                                                                       <  2:439606/251‑1 (MQ=255)
                                                                                                                 tGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCttttcgt                                                                                                                             >  1:64196/1‑251 (MQ=255)
                                                                                                                  ggACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCttttcgt                                                                                                                             <  1:566812/250‑1 (MQ=255)
                                                                                                                   gACGCGCGCGATACGCGCGCAACCTGTCACACTCGTAAGGAATGAGCTCCTCATCGGTATGCAGGATCTCAAGCCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCGGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCttttcgttt                                                                                                                           <  2:184414/251‑1 (MQ=255)
                                                                                                                           cGATACGCGCTCAACCCGTCACACTCTGAAGGAATGATCTCCGCATCGGTATGCATGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACTTCGAGTAAATCGCGATCAATTTGCTCTTCGTAGAAGATGCTCATGAGGAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGGTGCGTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTAt                                                                                                                   <  1:348219/251‑1 (MQ=255)
                                                                                                                             aTACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGt                                                                                                                                                                                                            >  2:412082/1‑160 (MQ=255)
                                                                                                                             aTACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGt                                                                                                                                                                                                            <  1:412082/160‑1 (MQ=255)
                                                                                                                              tACGCGCTCCACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTAtctc                                                                                                                >  1:49101/1‑251 (MQ=255)
                                                                                                                               aCGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTAtctc                                                                                                                <  1:248890/250‑1 (MQ=255)
                                                                                                                                cgcgCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCtt                                                                                                              >  2:429039/1‑251 (MQ=255)
                                                                                                                                cgcgCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCtt                                                                                                              <  1:346821/251‑1 (MQ=255)
                                                                                                                                   gcTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAg                                                                                                           <  2:405619/251‑1 (MQ=255)
                                                                                                                                   gcTCAACCCGTCAAACTCGTAAGGAATGATCTCCTCATCGGTATGCCGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGCGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCttt                                                                                                             >  2:81861/1‑249 (MQ=255)
                                                                                                                                          ccGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCgg                                                                                                    >  1:257158/1‑251 (MQ=255)
                                                                                                                                                    cGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCAtt                                                                                          <  1:143657/251‑1 (MQ=255)
                                                                                                                                                      tAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCAtttt                                                                                        >  2:124311/1‑251 (MQ=255)
                                                                                                                                                             tGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTcc                                                                                  >  2:570368/1‑250 (MQ=255)
                                                                                                                                                              gATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCa                                                                                 <  1:120632/250‑1 (MQ=255)
                                                                                                                                                                   cctcATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGctc                                                                             <  1:383111/249‑1 (MQ=255)
                                                                                                                                                                      cATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTCCCGGTCAGACCATTTTTTTTCCCGCTCtgtg                                                                         >  2:413233/1‑250 (MQ=255)
                                                                                                                                                                          ggTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCtt                                                                    <  1:259913/251‑1 (MQ=255)
                                                                                                                                                                           gTATGCAGGATCTCAAGTCCAGGGACAGGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGTCGCTCTGCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTg                                                                   <  1:418857/251‑1 (MQ=255)
                                                                                                                                                                             aTGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTc                                                                 >  1:77149/1‑251 (MQ=255)
                                                                                                                                                                                     tctcAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTa                                                          >  1:383805/1‑250 (MQ=255)
                                                                                                                                                                                        cAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGTCATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTTTATCTCTTTAGCTACCGGTCAGCCCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACta                                                      <  2:433199/251‑2 (MQ=255)
                                                                                                                                                                                             ccAGGGAAATGCTAACGCAGTGACCTCAGTACCGATGTGCGGTCGACATCGAGTAACGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCACTGt                                                 >  2:29715/1‑251 (MQ=255)
                                                                                                                                                                                               aGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTGCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTa                                               <  2:495185/251‑1 (MQ=255)
                                                                                                                                                                                                              cAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATa                                <  1:570368/251‑1 (MQ=255)
                                                                                                                                                                                                                gTGCCATCAATACCGATGGGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTTTGAACAGTTATTAAGCGGGCTTTTCGTGTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGTTTAACTCAATGTTAAATTGATGTAACATAAt                              <  2:390516/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCCATCAGTACCGATGTGCGGGCGACATAGAGGAAAGCGCCATCAAGACGCTCTTCGTTCAAGATGCTCATGAGTAGGCTTCGCTTTGTTGGGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTGCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATc                             <  2:269/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGTCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGTTTTTTCGTTTTCGTCTATCTCTTTAGCGACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGt                     <  1:360941/251‑1 (MQ=255)
                                                                                                                                                                                                                              gATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGAtgt                <  1:566543/251‑1 (MQ=255)
                                                                                                                                                                                                                               aTGTGCGGTCGACATCGAGAAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGAtgt                <  1:312737/250‑1 (MQ=255)
                                                                                                                                                                                                                                          aCATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGCGTGTTTTGCg    >  1:179752/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aTCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGCGTGTTTTGCGAg  <  1:439327/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
CCACCAGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGGGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGCGTGTTTTGCGAG  >  NC_000913/3127741‑3128227

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