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,9810GA98.5% 230.0 / NA 67P14L (CCC→CTC) glcDglycolate dehydrogenase, putative FAD‑linked subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/42);  minor base C (1/0);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCCACCAGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGGGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGCGTGTTTTGC  >  NC_000913/3127737‑3128224
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
cgcgcCACCAGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAg                                                                                                                                                                                                                                               <  1:124592/251‑1 (MQ=255)
  cgcCACCAGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGCGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAgcg                                                                                                                                                                                                                                             <  1:348632/251‑1 (MQ=255)
        caGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCa                                                                                                                                                                                                                                       <  1:160871/251‑1 (MQ=255)
             aGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACg                                                                                                                                                                                                                                  <  2:7966/251‑1 (MQ=255)
               cccGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGct                                                                                                                                                                                                                                >  1:157668/1‑251 (MQ=255)
                               cGGGTCCCCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGATATCGAGTAAAGCGCCAGTAAGACGCGCTTCGTACAAGATGCt                                                                                                                                                                                                                <  2:441566/251‑1 (MQ=255)
                                  gTCACCACCGTTACACGCAGGCGATGGCAGACAGCCAGAAGCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGTACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCAGCAAGACGCTCTTCGTACAAGATGCTCAt                                                                                                                                                                                                             <  2:387348/251‑1 (MQ=255)
                                        accGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTTCAAGATGCTCATGAGTAg                                                                                                                                                                                                       <  1:327951/251‑1 (MQ=255)
                                                acGCAGGCGATGGCAGCCAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGTCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGATTAGTCTTCGCt                                                                                                                                                                                               <  2:450220/251‑1 (MQ=255)
                                                acGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCt                                                                                                                                                                                               <  2:392346/251‑1 (MQ=255)
                                                  gCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCttt                                                                                                                                                                                             <  2:120337/251‑1 (MQ=255)
                                                   cAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTttg                                                                                                                                                                                            <  1:355051/251‑1 (MQ=255)
                                                    atgCGATAGCAGCCAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAAAAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTCAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTttgt                                                                                                                                                                                           >  2:295301/3‑251 (MQ=255)
                                                       cGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTtgttgt                                                                                                                                                                                        <  1:369973/251‑1 (MQ=255)
                                                             cAGACAGCCAGAATCGCTGTCCCCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTAGCTTTGttgtgttgtg                                                                                                                                                                                  >  1:322992/1‑251 (MQ=255)
                                                                              tGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTgc                                                                                                                                                                 >  1:77085/1‑251 (MQ=255)
                                                                                  aCCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTgcgcgc                                                                                                                                                             <  1:320649/251‑1 (MQ=255)
                                                                                   ccTGTTCCATTTGCTTAGGCAGACCAACCAGGAATGGACGCGTTCGATACGCGCTCAACCCGTCACACTCTTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCgct                                                                                                                                                            <  2:19269/251‑1 (MQ=255)
                                                                                    ctgttacctttGCTTAGGCAGAACCCCCAGTACTGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGAATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGCTGTGCGGCCGACATCGAGTAAAgc                                                                                                                                                                                                                                              >  2:267727/9‑168 (MQ=255)
                                                                                    cTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATTATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAgc                                                                                                                                                                                                                                              <  1:267727/168‑1 (MQ=255)
                                                                                     tGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCgctgc                                                                                                                                                          >  1:261883/1‑251 (MQ=255)
                                                                                      gTTCCATTTGCTTAGGCAGAACAACCAGTACTGGCCGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGGGCGGGCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGATGTGTTGTGTGGCAGCTGATTTTTTAGCgctgct                                                                                                                                                         >  1:117264/1‑251 (MQ=255)
                                                                                             ttGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGTTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGTAGCTGATTTTTGCGCGCTGCTTCtgtg                                                                                                                                                   <  1:19553/250‑1 (MQ=255)
                                                                                             ttGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCtgtg                                                                                                                                                   <  1:61140/250‑1 (MQ=255)
                                                                                               gCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAAc                                                                                                                                                <  1:4662/251‑1 (MQ=255)
                                                                                               gCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTGGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAAc                                                                                                                                                <  1:401675/251‑1 (MQ=255)
                                                                                                   aGGCAGAACAACCAGTAATGGACGTGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCGCTTCGTTCAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGTCAGCTGATTTTTGCTCGCTGCGTCTGTGAACAGtt                                                                                                                                            <  1:77947/251‑1 (MQ=255)
                                                                                                   aGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGt                                                                                                                                             <  1:129916/250‑1 (MQ=255)
                                                                                                       aGAACAACCAGTAATGTACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTTGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCGTCTGTGAACAGttatta                                                                                                                                        <  2:218895/251‑1 (MQ=255)
                                                                                                       aGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGttatt                                                                                                                                         >  1:355754/1‑250 (MQ=255)
                                                                                                        gAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGttatta                                                                                                                                        >  1:133006/1‑250 (MQ=255)
                                                                                                         aacaacCAGTAATGGACGCATGCGATACGCGCCAAACCCGTCACACTCGTACGGAATGATCTCCTCATCGGTATGCCGGCTCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGCGTAAAGCGCCCTCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGCACAGttatta                                                                                                                                        >  2:295080/1‑249 (MQ=255)
                                                                                                          acaacCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCt                                                                                                                                                                                                    <  1:296679/188‑1 (MQ=255)
                                                                                                             acCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCggg                                                                                                                                  >  1:260782/1‑251 (MQ=255)
                                                                                                                    aTGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGCAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCttttcg                                                                                                                           >  2:86215/1‑251 (MQ=255)
                                                                                                                     tGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCttttcgt                                                                                                                          >  1:301536/1‑251 (MQ=255)
                                                                                                                     tGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGAGGACGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCttttcgt                                                                                                                          >  1:202669/1‑251 (MQ=255)
                                                                                                                      ggACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCttttcgtt                                                                                                                         <  1:432106/251‑1 (MQ=255)
                                                                                                                                gATACGCGCTCAACCCGTCACACTCGTAAGGACTGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCCTGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGAGCGCTGCTTCTGTGCACAGTTATTAAGCGGGCTTTTCGTTTTCGTCtctc                                                                                                               >  2:166468/1‑248 (MQ=255)
                                                                                                                                gATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCGTTTCGTTTTCGTCTAtc                                                                                                               <  1:384390/251‑1 (MQ=255)
                                                                                                                                   aCGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATctct                                                                                                            <  2:261883/251‑1 (MQ=255)
                                                                                                                                     gcgcTCAACCCGTCACACTCGTAAGGAATTATCTCCTCAGCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCttt                                                                                                          <  2:129501/251‑1 (MQ=255)
                                                                                                                                             cccGTCACACGCGTAAGGAATGATCTCCTCATCGGTATGCAGTATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTTAGTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCg                                                                                                  <  1:232066/251‑1 (MQ=255)
                                                                                                                                                gTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTAAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTAACGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTc                                                                                               >  1:277793/1‑251 (MQ=255)
                                                                                                                                                                             cGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTTGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGAcc                                                                   <  2:133006/250‑1 (MQ=255)
                                                                                                                                                                              ggTATGCAGGATCTCAAGGCCAGGGACATGCTCACGCAGTGCCATCAGTTCCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCGCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCTTCTATCTCTTTAGCGACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCtt                                                                 <  2:336304/251‑1 (MQ=255)
                                                                                                                                                                                    cAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTg                                                           >  1:114683/1‑251 (MQ=255)
                                                                                                                                                                                     aGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTTAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTgg                                                          <  1:17056/251‑1 (MQ=255)
                                                                                                                                                                                          ctcAAGTCCAGGGCCATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAAc                                                     >  1:79769/1‑251 (MQ=255)
                                                                                                                                                                                            cAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTc                                                   >  1:6939/1‑251 (MQ=255)
                                                                                                                                                                                                 ccAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGt                                              <  2:6939/251‑1 (MQ=255)
                                                                                                                                                                                                        cATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTg                                       >  1:274929/1‑251 (MQ=255)
                                                                                                                                                                                                           gCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGCCATCGAGTAAAGCGACATCAAGACGCTCTTCGTACAAGATGCTCCTGAGTAGGCTTCGATTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCTGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTTGTTAACTCAATGTTAAATTgctt                                    >  2:40903/1‑248 (MQ=255)
                                                                                                                                                                                                           gCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATg                                    >  2:349865/1‑251 (MQ=255)
                                                                                                                                                                                                              cACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTa                                  >  2:411516/1‑250 (MQ=255)
                                                                                                                                                                                                              cACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTa                                  <  2:114683/250‑1 (MQ=255)
                                                                                                                                                                                                                cgccgTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATTCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCAGTCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTTGTTAACTCAATGTTAAATTGATGTAACa                               <  2:408607/247‑1 (MQ=255)
                                                                                                                                                                                                                 gCAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACAt                              <  1:116749/251‑1 (MQ=255)
                                                                                                                                                                                                                  cAGTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATa                             <  1:442703/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTGCCATCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAAt                           <  1:212636/251‑1 (MQ=255)
                                                                                                                                                                                                                       ccagcAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCa                         <  2:113855/246‑1 (MQ=255)
                                                                                                                                                                                                                         aTCAGTACCGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACtt                      <  2:77085/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGATGTGCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGAt               <  1:92611/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tgtgCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGc           >  1:449135/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gtgCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGc           <  1:182370/250‑1 (MQ=255)
                                                                                                                                                                                                                                      tgCGGTCGACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGCgt         >  1:52001/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gACATCGAGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGGTTTCGTCTATCTCGTTAGCTACCGGTCAGACCATTTTTTGTCCAGCTCTGTGACCTTGTCTTGGTTAACGCATTGTTAAATTGATGTAACATAATCACTTACGTGATGTGCGTGTTTTGc  <  1:295080/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CGCGCCACCAGAAAGCCCGGTGCCTGCACCACGGGTCACCACCGGTACACGCAGGCGATGGCAGACAGCCAGAATCGCTGTCACCTGTTCCATTTGCTTAGGCAGAACAACCAGTAATGGACGCGTGCGATACGCGCTCAACCCGTCACACTCGTAAGGAATGATCTCCTCATCGGTATGCAGGATCTCAAGTCCAGGGACATGCTCACGCAGTGCCATCAGTACCGATGTGCGGTCGACATCGGGTAAAGCGCCATCAAGACGCTCTTCGTACAAGATGCTCATGAGTAGGCTTCGCTTTGTTGTGTTGTGTGGCAGCTGATTTTTGCGCGCTGCTTCTGTGAACAGTTATTAAGCGGGCTTTTCGTTTTCGTCTATCTCTTTAGCTACCGGTCAGACCATTTTTTTTCCAGCTCTGTGACCTTGTCTTGGTTAACTCAATGTTAAATTGATGTAACATAATCACTTACGTGATGTGCGTGTTTTGC  >  NC_000913/3127737‑3128224

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