Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,840,287 G→T 13.7% V256L (GTG→TTG)  ascF → cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,840,2870GT13.7% 174.0 / 2.3 71V256L (GTG→TTG) ascFcellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components
Reads supporting (aligned to +/- strand):  ref base G (20/40);  new base T (9/2);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 5.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGCCGTCAGAG  >  NC_000913/2840056‑2840523
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
cTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGGCCAAAAACTTCCTCAAACCCATGCTGATTGTGTTGAt                                                                                                                                                                                                                                              >  2:948032/1‑239 (MQ=255)
 tGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGGTGCCAAAAACTTCCTCCAACCCATGCTGATTTGTTTGAtt                                                                                                                                                                                                                                             >  2:1012783/1‑239 (MQ=255)
                     tgcgtggagcctCCGGGTTTTGTTGTACGGGTGGCCAAACGGGGCCGGGGTGAACAGTTGTAATTGGCCGCGATTCCTGGCACCGGGGTGAAATGCCCCTACCAGGTGATGCCGGCGCTGGGGATGAGCTGGGGCCTGTCATATATCGAACGCTGGGTGGACAGCATGCCACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCct                                                                                                                                                                                                                         <  1:233339/228‑1 (MQ=255)
                      tgCTGGTCCACCCGAGGTTTATTGAACGGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGAAAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCctg                                                                                                                                                                                                                        <  1:810304/239‑1 (MQ=255)
                                tCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACCCCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGcc                                                                                                                                                                                                              <  2:480668/239‑1 (MQ=255)
                                tCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCCAAAACTTCCTCAAACCCATACTGATTTTGTTGATTTCCGCCCCCCTGGCAATACTTCTGATTGGc                                                                                                                                                                                                               >  2:141751/1‑238 (MQ=255)
                                cccGGGTTTTTTTGAATTTGTGGCGAAGGCGGCCCAGGGTTAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGcc                                                                                                                                                                                                              <  1:564099/238‑1 (MQ=255)
                                 ccGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCCGATTTTGTTGATTGCCGCCCCGCCGGCAATccctcctgcccctta                                                                                                                                                                                                             >  2:707494/1‑226 (MQ=255)
                                     gTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCCAACCAATTCTGTTTTTGTTGATTGCCGCACCCCGGGCAAACCCGCTACTTGGCCCGCtt                                                                                                                                                                                                         >  2:503855/1‑239 (MQ=255)
                                       tttATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTg                                                                                                                                                                                                        <  2:899823/238‑1 (MQ=255)
                                         tagagAATTGATGGCCAAAGCGGCCTAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATCCACCTACCCGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCCCCGCTGGCAATCCTGCTGATTGGCCCGATTGGTa                                                                                                                                                                                                    <  2:381788/235‑1 (MQ=255)
                                                    gggcgaaggggcgccAGGGTGATCATGTCCAATTTGCCCTTCTTCCGGTCACCGCGGGGAAAAGCCCCTTCGCGGGGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTTGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTTGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGt                                                                                                                                                                                         <  1:699414/226‑1 (MQ=255)
                                                             cGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGTTTTGCCCGATTGGTATCTGGGTCGGTAGCGCCCtt                                                                                                                                                                               >  2:205673/1‑239 (MQ=255)
                                                                   gggggGAACATGTGGAATTTGCCCGTGTTCCGGTCCCCGCGGTGAAATACACCTTCACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGgcg                                                                                                                                                                         <  1:948032/238‑1 (MQ=255)
                                                                   ggggTGAACATGGCAAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGgcg                                                                                                                                                                         <  2:633116/238‑1 (MQ=255)
                                                                    ggggGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATCCACCTACGCGGTGAGCCCGGCGCTGGTGATGGTCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGgcgc                                                                                                                                                                        <  2:589114/239‑1 (MQ=255)
                                                                      ttGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGgcg                                                                                                                                                                         <  2:601066/235‑1 (MQ=255)
                                                                          aCATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGATGATTTTGTTGATTGCCGCACCGCTGGCCATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGGGCCCTTTTGGCGCCGGgtt                                                                                                                                                                 >  2:291683/1‑239 (MQ=255)
                                                                          aCATGGCGAATTTCCCCTGATTCCGGACCCCGGGGTGAAATCCCCCTACACGGTGAACCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACTCCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGtt                                                                                                                                                                  <  2:988739/238‑1 (MQ=255)
                                                                            ctGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACCCCTACACGGTGATCCCGGCTCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTac                                                                                                                                                               <  2:735655/238‑1 (MQ=255)
                                                                             ttttgAATTTGCCCTGATTGCGGTCACCGCGGTTAAATACCCCTACTCGGTGATCCCGGCGCTGGTGATGACCTGTTGCCTGTCATATATCGAACGCTGGGTGGAAAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTaca                                                                                                                                                              <  2:197838/235‑1 (MQ=255)
                                                                             tGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACGCGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTaca                                                                                                                                                              <  2:997766/239‑1 (MQ=255)
                                                                             gttccaaTTGGCCCTGATTTCCTTCCCCGCTGTGAAATTCACCTACGCGGGGAGCCCGCCGCTGGTGATGGCCTGGGGGATGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGCCCAAAAACTTCCTCAACCCAAGGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGGTTGGCCCGATTGGTATCTGGATCTGTAGCGCCATTTCGGCGCTGGTTTaca                                                                                                                                                              <  1:563600/234‑1 (MQ=255)
                                                                              gTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCCCACCGCTTGCAATCCTGCTGGTTGGCCCCATTGGGTTCTGGATCGGTAGCGCCCTTTCCgccccgggtttcc                                                                                                                                                               >  2:752704/1‑227 (MQ=255)
                                                                               tCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCCCACCGCTGGCCATCCTGCTGGTTGGCCCGAATGGATATCGGATCGGTGGGGCCATTTCgcggtcgggttaaacc                                                                                                                                                             >  2:168007/1‑224 (MQ=255)
                                                                                      tGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCCAAAACCTCCCCAAACCAATGCTGATTGGGTTTCTTGCCGCCCCCCCGGCAATCCTGCCGATTTGCCCGATTGGGCTTCGGAGTGGGAGCGCCATTTccggcccggggtaccacagtcaag                                                                                                                                                      >  1:1035081/1‑216 (MQ=255)
                                                                                                        cgcgGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCCGTAGCGCCATTTCTGCGCTGGTTTACACCCTTCATGGTTTTTTGGGCTGGCttt                                                                                                                                >  2:788948/1‑239 (MQ=255)
                                                                                                           ggTTAAATCCACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAg                                                                                                                             <  2:792175/239‑1 (MQ=255)
                                                                                                              gAAATACCCCTTCACGGGGATCCCGGCGCTGGTTGTGACCTTGTGCCTGTAAAATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCACAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTc                                                                                                                           <  1:1017495/238‑1 (MQ=255)
                                                                                                              gAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTc                                                                                                                           <  2:1241041/238‑1 (MQ=255)
                                                                                                                     acCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGtt                                                                                                                                                                                                                                                 <  2:976916/119‑1 (MQ=255)
                                                                                                                      ccTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCCCCATTTCGGCGCCGGTTTTCACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATg                                                                                                                 >  2:615340/1‑239 (MQ=255)
                                                                                                                       cTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGAGACCCAAAACTTCCTCAAACCAAAGATGGTTGTGTTGGATTCCGCACCGCTGGGAATACTGTCGATTGGCCCGGTTGGGATCTGGGTCGGTATGGCCATTTCGGGGCTGGTTTACACCACTCAAGGGTATCGTGGGTGGCTTTCAAGCGCCATTATgg                                                                                                                >  1:129856/1‑239 (MQ=255)
                                                                                                                         acacGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGGATCTGGATCGGTAGCGCCCTTTCGGCGCTGGTTTTAACCATTCATGGTTATCTGGGCCGGCTTTCAATCGCCATTTTTGgc                                                                                                              >  2:343264/1‑239 (MQ=255)
                                                                                                                          cccGGTGATCCCGTCGCTGGTTATGCCCGTGTGCCTGTCATATATCGAACGCTGGTTGGACAGCATTCCACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGgcg                                                                                                             <  1:631307/237‑1 (MQ=255)
                                                                                                                          cacGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGCGGACAGCATTACACCGGCGGCGCTCCAAAACCTCCTCAAAACACTTGTAGTTGTGTTTATTTCCGCCACGCTTGCACTCCTGCTTGTTGGCCCGCGTTGTGTTTGGCTCCGTAACCGCCTTTCGTCGCGTGGTTTACCCCATCCGTGTGATCTGGGCGGGGGTTCTgccgcccaatgtgcgg                                                                                                              >  1:1071033/1‑229 (MQ=255)
                                                                                                                               tGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAAATTCCTCAAACCAATGCTGATTTGGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTTTCTGGATCGGTAGCGCCCTTTCGGCGCTGGGTTACACCATTTATGGTTTATCTGGCTGGGTTTACGTCGCCCTTTTTGgcgggccg                                                                                                        >  2:96728/1‑237 (MQ=255)
                                                                                                                                    gcGGCGCTGGGTATGCCCTGGTGCGTGTCATATATCGACCGCTGGTTGGTCAGGATTATGCCGGAGGGGACCAAAACCGTCCTCAACCCCATGCTGATGGTGTTGATTGCCGGCCCGGTGGCAAGCCTGCGGATGGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGcc                                                                                                   <  1:844362/238‑1 (MQ=255)
                                                                                                                                    cggtgggtgtgtttggccTGGGGCCTGTAGAAGAGCGAACCCTGGGGGGAAAGGATTACCCCGGCGGGGACCAAAAACTTCCTCAAACCAATGCTAATTGTGGTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGGCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGcc                                                                                                   <  1:141751/223‑1 (MQ=255)
                                                                                                                                          cTGGTGATGACCTGGTGCCTGGCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTctgc                                                                                              <  2:613935/238‑1 (MQ=255)
                                                                                                                                             gTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGTCCAAAAAATTCCTCAAACCAATGGTGGTTGTGTTGATTGCCGCCCCGCCTGCAATCCTGCTGATTGGCCCGGTTTGTATCTGGGGCGGGCAGGCCCCTTCGGCGCTGGGTTAAACCCATCCTGGTTTTTCGGGGTGGCTTTTAGGTGCCCTTATTGGGCGCCCTGGGGctgttggggg                                                                                           >  1:529713/1‑230 (MQ=255)
                                                                                                                                             gTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGATGGGTGGGTTGATTTGCGCCCCGCCGGGAAATATTGTGGATGTACCGCTTGGTATTTGGGTCGGTGGGCCACTTTTGGAGGTCGGTTAAACCATTTCTGGTATTGTGGGGCTGGCTTTGTTGGCCATGATGGGGGCCGTGTGGCctctctgggg                                                                                          >  1:1085423/1‑233 (MQ=255)
                                                                                                                                                 tgaCCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTTTGTTGATTGCCGCAACGCTGGCCAACCTGCTGATTGGCCCGGTTGGGTTCTGGATCGGTAGAGCCATTTTGGCGCTGGTTTACAACAtttctgg                                                                                                                                                    >  2:412640/1‑176 (MQ=255)
                                                                                                                                                  gatctgtgtcctgatatataaCTAACGCGGTGGGGAGAGCATTTCCCCGGCGGTGGACAAAATACTCCTCAAACCAATGCAGGTTGTGTTGAGTGCCGACCCGCTGGCAATCCTGCTGATGGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGa                                                                                     <  1:205673/224‑1 (MQ=255)
                                                                                                                                                  aaCCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGGCCCACCATTTCCTACAACCTAAGCCTGATGAGGTGATTGCCGCCCCCCATGCAAATCTTCTTACTGTTACCATTGGTTTTTGCAGCGGTACCGCAGTTTCGACGACTGTTTTCACCATTTATTGTTATCAGGTGTGTGCTCGCGCACGCATTATGGGCGGCCTGTGCTCTCTGCTggcacagc                                                                                     >  1:360268/2‑233 (MQ=255)
                                                                                                                                                       ggTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGGGACCAAAAACTTCCTCCAAACCATGCTGATTTTGTTGATTGCCGCCCCGCTGGCAATCCTGCTGATCTGCCCGATTGGTGTCTGGGTCGGTAGCGCCGTTTCGCTGCCGGGTTATACCCTTCCTGGTTATTTGGGCGGGCTTTCAGCCGCCCTTCAGGGCCGGCTGGGGGGCCTGCGGGGAATAACgggg                                                                                >  1:577361/1‑239 (MQ=255)
                                                                                                                                                         ggccttgaaaatatAGAACGCTGGGGGGAAAGGTTTCCCCCGGCGGGAACCAAAAACTTCCTCAACCCAAGGCTGATTGTGTTGATTGGCGGACCGCTGGCAATCCGTCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTGCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTGTGCTGGTAATGACCGGGa                                                                               <  1:752704/228‑1 (MQ=255)
                                                                                                                                                               tCATATATCGAACGCTGGGTGGACAGCATTACCCCGGCGGTGACCAAAAGCTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTctg                                                                                               <  2:282293/216‑1 (MQ=255)
                                                                                                                                                               tCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAAt                                                                                       <  2:162893/224‑1 (MQ=255)
                                                                                                                                                                tttatttctacccttGGGTGGCAAGCATTCCGCCGGCGGTTACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCGGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACcg                                                                        <  1:343264/224‑1 (MQ=255)
                                                                                                                                                                aaaataaCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGGTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACcg                                                                        <  2:1071033/235‑1 (MQ=255)
                                                                                                                                                                     atCGAACGCTGGGTGGACAGCATTACCCCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCtt                                                                  <  2:129856/239‑1 (MQ=255)
                                                                                                                                                                         aaCGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTTTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGGTTGGTATCTGGATCGGTAGCGCCCTTTCTGCGCTGGTTTACACCATTCATGGTTATCTTGGCTGGCTTTCAGTCGCCATTTTGGGCGCGCTGTGGCCCCTGCTGGTAATAACCGGGGTGTCCCGCGTCttttg                                                               >  2:119283/1‑236 (MQ=255)
                                                                                                                                                                         aaCGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTc                                                                                                                           <  2:1279146/179‑1 (MQ=255)
                                                                                                                                                                           cgtggggTGGCCAGCATTACCCCGTCGGTGACCAAACACGTCCTCAAACCAATGCGGATTGTGTGGATTGCCGTACCTGTGGCAGTCCGGGTGATGGGCCCGGTTGGTAACTGGCTCGGTAGCGCCATTTTGGTGCTGGTTTACACCATTCATGGTTATCTGGGCGGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCGTTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGcc                                                            <  1:1035090/235‑1 (MQ=255)
                                                                                                                                                                           cgcgGGGTGGACAGCATTACCCCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGcc                                                            <  2:529713/239‑1 (MQ=255)
                                                                                                                                                                               gggTGGACAGCATTCCGCCGGGGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAAcc                                                        <  2:635455/239‑1 (MQ=255)
                                                                                                                                                                                 ggtgACAGTATTACCGCGGCGGGGACCACAATCTCCCTCAAACCAATGGTGATTTTGTGGATTGCCGCTCCGCGGGCAATCCTGCTGATTGGCCCGATTGGTATGTGGGTCGGTAGCGCCCTTTCGGCGCTGGTTTACACCATGCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACca                                                       <  1:96728/235‑1 (MQ=255)
                                                                                                                                                                                 gTGGACAGCATTTCACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACca                                                       <  2:8760/238‑1 (MQ=255)
                                                                                                                                                                                 gTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACca                                                       <  2:1182250/238‑1 (MQ=255)
                                                                                                                                                                                         cATTACCCCGGCGGTGACCAAAACCGTCCTCAAACCAAGGCTGATTGGGTTGATTGCCGCCCCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACGCCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACcg                                                                        <  2:1102593/213‑1 (MQ=255)
                                                                                                                                                                                         cATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATc                                                                                                                                              <  2:840315/146‑1 (MQ=255)
                                                                                                                                                                                                   gcggTGACCAAAAACTTCCTCAAAACAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGGCCGGTTGGGATCTTGATCCGTAGCGCCATTTTGGCGCCGGTTTTCACCATTCATGGGTATCTGggctggct                                                                                                                                  >  1:840315/1‑146 (MQ=255)
                                                                                                                                                                                                              aaaCTTCCTCAAAACCATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGTCCCGATTGGGATCTGGATCGGTAGCGCCATTTCGCAGCTGGTTTACACCATTCATGGTTATTCGGGATGCGTTTCAGTCGGCATTATTGGCTCTCTGTGTCCTCTGCTGGTAATGCCCGGGTGGCACCTCGCCTTTTCGCCACCCACCCTTCAAACCATTTGCGCAAACGGCCaag                         >  1:647879/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCTTCATCCAGACCATTGCCGAAACCGGCAAAGAAGGGa                   >  2:667271/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGa                   >  2:1053537/1‑239 (MQ=255)
                                                                                                                                                                                                                     cTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGTTAGCGCCATTTCGgcg                                                                                                                                                                         >  2:123810/1‑93 (MQ=255)
                                                                                                                                                                                                                     cTCAAACCAATGATGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGgcg                                                                                                                                                                         <  1:123810/93‑1 (MQ=255)
                                                                                                                                                                                                                              aaTGGTGATTTTGGTTATTGCCGCACCGGTGGCAATCCTGCTGATTTGCCCGATTGGTATCTGGATCGGGAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGc          <  2:168761/239‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGGAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGCCGTCaga   >  2:1232290/1‑239 (MQ=255)
                                                                                                                                                                                                                                      ttGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCCAAGAAGGGATGGTCAGGCCGCCagag  >  2:805947/1‑239 (MQ=255)
                                                                                                                                                                                                                                       tgtgTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCAATCCGGCGCTGCTTTACACCATTCATGGTTATCTGGGCTGGCTGTCAGTCGCCATTCTGGCCGCGCTTTGGCCTCTGCTGGTAATGACCGGGATGTCACGCGTCTTTAAGCCACCCATCATTCAGCCCATTGCACAAACCGGCAaagaaagagaggaaaggacgacaagg  >  1:575932/1‑217 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
CTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCTCAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGGTATCTGGATCGGTAGCGCCATTTCGGCGCTGGTTTACACCATTCATGGTTATCTGGGCTGGCTTTCAGTCGCCATTATGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGCCGTCAGAG  >  NC_000913/2840056‑2840523

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