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

CTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGCCGTCAGAG  >  NC_000913/2840056‑2840523
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
cTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGGCCAAAAACTTCCT‑CAAACCCATGCTGATTGTGTTGAt                                                                                                                                                                                                                                              >  2:948032/1‑239 (MQ=255)
 tGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGGTGCCAAAAACTTCCT‑CCAACCCATGCTGATTTGTTTGAtt                                                                                                                                                                                                                                             >  2:1012783/1‑239 (MQ=255)
                     tgcgtggagcctCCGGGTTTTGTTGTACGGGTGGCCAAACGGGGCCGGGGTGAACAGTTGTAATTGGCCGCGATTCCTGGCACCGGGGTGAAATGCCCCTACCAGGTGATGCCGGCGCTGGGGATGAGCTGGGGCCTGTCATATATCGAACGCTGGGTGGACAGCATGCCACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCct                                                                                                                                                                                                                         <  1:233339/228‑1 (MQ=255)
                      tgCTGGTCCACCCGAGGTTTATTGAACGGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGAAAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCctg                                                                                                                                                                                                                        <  1:810304/239‑1 (MQ=255)
                                tCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACCCCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGcc                                                                                                                                                                                                              <  2:480668/239‑1 (MQ=255)
                                tCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCCAAAACTTCCT‑CAAACCCATACTGATTTTGTTGATTTCCGCCCCCCTGGCAATACTTCTGATTGGc                                                                                                                                                                                                               >  2:141751/1‑238 (MQ=255)
                                cccGGGTTTTTTTGAATTTGTGGCGAAGGCGGCCCAGGGTTAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGcc                                                                                                                                                                                                              <  1:564099/238‑1 (MQ=255)
                                 ccGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCCGATTTTGTTGATTGCCGCCCCGCCGGCAATccctcctgcccctta                                                                                                                                                                                                             >  2:707494/1‑226 (MQ=255)
                                     gTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CCAACCAATTCTGTTTTTGTTGATTGCCGCACCCCGGGCAAACCCGCTACTTGGCCCGCtt                                                                                                                                                                                                         >  2:503855/1‑239 (MQ=255)
                                       tttATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTg                                                                                                                                                                                                        <  2:899823/238‑1 (MQ=255)
                                         tagagAATTGATGGCCAAAGCGGCCTAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATCCACCTACCCGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCCCCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑Ta                                                                                                                                                                                                    <  2:381788/235‑1 (MQ=255)
                                                    gggcgaaggggcgccAGGGTGATCATGTCCAATTTGCCCTTCTTCCGGTCACCGCGGGGAAAAGCCCCTTCGCGGGGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTTGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTTGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGt                                                                                                                                                                                         <  1:699414/226‑1 (MQ=255)
                                                             cGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGTTTTGCCCGATTGG‑TATCTGGGTCGGTAGCGCCC‑tt                                                                                                                                                                               >  2:205673/1‑239 (MQ=255)
                                                                   gggggGAACATGTGGAATTTGCCCGTGTTCCGGTCCCCGCGGTGAAATACACCTTCACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGgcg                                                                                                                                                                         <  1:948032/238‑1 (MQ=255)
                                                                   ggggTGAACATGGCAAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGgcg                                                                                                                                                                         <  2:633116/238‑1 (MQ=255)
                                                                    ggggGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATCCACCTACGCGGTGAGCCCGGCGCTGGTGATGGTCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGgcgc                                                                                                                                                                        <  2:589114/239‑1 (MQ=255)
                                                                      ttGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGgcg                                                                                                                                                                         <  2:601066/235‑1 (MQ=255)
                                                                          aCATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGATGATTTTGTTGATTGCCGCACCGCTGGCCATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGGGCCC‑TTTTGGCGC‑CGGgtt                                                                                                                                                                 >  2:291683/1‑239 (MQ=255)
                                                                          aCATGGCGAATTTCCCCTGATTCCGGACCCCGGGGTGAAATCCCCCTACACGGTGAACCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACTCCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGtt                                                                                                                                                                  <  2:988739/238‑1 (MQ=255)
                                                                            ctGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACCCCTACACGGTGATCCCGGCTCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTac                                                                                                                                                               <  2:735655/238‑1 (MQ=255)
                                                                             ttttgAATTTGCCCTGATTGCGGTCACCGCGGTTAAATACCCCTACTCGGTGATCCCGGCGCTGGTGATGACCTGTTGCCTGTCATATATCGAACGCTGGGTGGAAAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTaca                                                                                                                                                              <  2:197838/235‑1 (MQ=255)
                                                                             tGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACGCGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTaca                                                                                                                                                              <  2:997766/239‑1 (MQ=255)
                                                                             gttccaaTTGGCCCTGATTTCCTTCCCCGCTGTGAAATTCACCTACGCGGGGAGCCCGCCGCTGGTGATGGCCTGGGGGATGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGCCCAAAAACTTCCT‑CAACCCAAGGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGGTTGGCCCGATTGG‑TATCTGGATCTGTAGCGCCA‑TTTCGGCGC‑TGGTTTaca                                                                                                                                                              <  1:563600/234‑1 (MQ=255)
                                                                              gTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCCCACCGCTTGCAATCCTGCTGGTTGGCCCCATTGG‑GTTCTGGATCGGTAGCGCCC‑TTTCCgccccgggtttcc                                                                                                                                                               >  2:752704/1‑227 (MQ=255)
                                                                               tCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCCCACCGCTGGCCATCCTGCTGGTTGGCCCGAATGGATATC‑GGATCGGTGGGGCCA‑TTTCgcggtcgggttaaacc                                                                                                                                                             >  2:168007/1‑224 (MQ=255)
                                                                                      tGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCCAAAACCTCCC‑CAAACCAATGCTGATTGGGTTTCTTGCCGCCCCCCCGGCAATCCTGCCGATTTGCCCGATTGG‑GCTTCGGAGTGGGAGCGCCA‑TTTccggcccggggtaccacagtcaag                                                                                                                                                      >  1:1035081/1‑216 (MQ=255)
                                                                                                        cgcgGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCCGTAGCGCCA‑TTTCTGCGC‑TGGTTTACACCCTTCATGG‑TTTTTT‑GGGCTG‑GCttt                                                                                                                                >  2:788948/1‑239 (MQ=255)
                                                                                                           ggTTAAATCCACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAg                                                                                                                             <  2:792175/239‑1 (MQ=255)
                                                                                                              gAAATACCCCTTCACGGGGATCCCGGCGCTGGTTGTGACCTTGTGCCTGTAAAATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCA‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTc                                                                                                                           <  1:1017495/238‑1 (MQ=255)
                                                                                                              gAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTc                                                                                                                           <  2:1241041/238‑1 (MQ=255)
                                                                                                                     acCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGtt                                                                                                                                                                                                                                                 <  2:976916/119‑1 (MQ=255)
                                                                                                                      ccTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCCCCA‑TTTCGGCGC‑CGGTTTTCACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑Tg                                                                                                                 >  2:615340/1‑239 (MQ=255)
                                                                                                                       cTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGAGACCCAAAACTTCCT‑CAAACCAAAGATGGTTGTGTTGGATTCCGCACCGCTGGGAATACTGTCGATTGGCCCGGTTGG‑GATCTGGGTCGGTATGGCCA‑TTTCGGGGC‑TGGTTTACACCACTCAAGG‑GTATCG‑TGGGTG‑GCTTTCAAGCGCCATTA‑Tgg                                                                                                                >  1:129856/1‑239 (MQ=255)
                                                                                                                         acacGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑GATCTGGATCGGTAGCGCCC‑TTTCGGCGC‑TGGTTTTAACCATTCATGG‑TTATCT‑GGGCCG‑GCTTTCAATCGCCATTT‑TTGgc                                                                                                              >  2:343264/1‑239 (MQ=255)
                                                                                                                          cccGGTGATCCCGTCGCTGGTTATGCCCGTGTGCCTGTCATATATCGAACGCTGGTTGGACAGCATTCCACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGgcg                                                                                                             <  1:631307/237‑1 (MQ=255)
                                                                                                                          cacGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGCGGACAGCATTACACCGGCGGCGCTCCAAAACCTCCT‑CAAAACACTTGTAGTTGTGTTTATTTCCGCCACGCTTGCACTCCTGCTTGTTGGCCCGCGTTG‑TGTTTGGCTCCGTAACCGCC‑TTTCGTCGCGTGGTTTACCCCATCCGT‑G‑TGATCT‑GGGCGG‑GGGTTCTgccgcccaatgtgcgg                                                                                                              >  1:1071033/1‑229 (MQ=255)
                                                                                                                               tGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAAATTCCT‑CAAACCAATGCTGATTTGGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TTTCTGGATCGGTAGCGCCC‑TTTCGGCGC‑TGGGTTACACCATTTATGGTTTATCT‑‑GGCTG‑GGTTTACGTCGCCCTTT‑TTGgcgggccg                                                                                                        >  2:96728/1‑237 (MQ=255)
                                                                                                                                    gcGGCGCTGGGTATGCCCTGGTGCGTGTCATATATCGACCGCTGGTTGGTCAGGATTATGCCGGAGGGGACCAAAACCGTCCT‑CAACCCCATGCTGATGGTGTTGATTGCCGGCCCGGTGGCAAGCCTGCGGATGGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGcc                                                                                                   <  1:844362/238‑1 (MQ=255)
                                                                                                                                    cggtgggtgtgtttggccTGGGGCCTGTAGAAGAGCGAACCCTGGGGGGAAAGGATTACCCCGGCGGGGACCAAAAACTTCCT‑CAAACCAATGCTAATTGTGGTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGGCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGcc                                                                                                   <  1:141751/223‑1 (MQ=255)
                                                                                                                                          cTGGTGATGACCTGGTGCCTGGCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTctgc                                                                                              <  2:613935/238‑1 (MQ=255)
                                                                                                                                             gTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGTCCAAAAAATTCCT‑CAAACCAATGGTGGTTGTGTTGATTGCCGCCCCGCCTGCAATCCTGCTGATTGGCCCGGTTTG‑TATCTGGGGCGGGCAGGCCC‑CTTCGGCGC‑TGGGTTAAACCCATCCTGG‑TTTTTC‑GGGGTG‑GCTTTTAGGTGCCCTTATTGGGCGCCCTGGGGctgttggggg                                                                                           >  1:529713/1‑230 (MQ=255)
                                                                                                                                             gTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGATGGGTGGGTTGATTTGCGCCCCGCCGGGAAATATTGTGGATGTACCGCTTGG‑TATTTGGGTCGGT‑GGGCCACTTTTGGAGG‑TCGGTTAAACCATTTCTGG‑TATTGTGGGGCTG‑GCTTT‑GTTGGCCATGA‑TGGGGGCCGTGTGGCctctctgggg                                                                                          >  1:1085423/1‑233 (MQ=255)
                                                                                                                                                 tgaCCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTTTGTTGATTGCCGCAACGCTGGCCAACCTGCTGATTGGCCCGGTTGG‑GTTCTGGATCGGTAGAGCCA‑TTTTGGCGC‑TGGTTTACAACAtttctgg                                                                                                                                                    >  2:412640/1‑176 (MQ=255)
                                                                                                                                                  gatctgtgtcctgatatataaCTAACGCGGTGGGGAGAGCATTTCCCCGGCGGTGGACAAAATACTCCT‑CAAACCAATGCAGGTTGTGTTGAGTGCCGACCCGCTGGCAATCCTGCTGATGGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGa                                                                                     <  1:205673/224‑1 (MQ=255)
                                                                                                                                                  aaCCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGGCCCACCATTTCCTACAACCTAAGCCTGA‑TGAGGTGATTGCCGCCCCCCATGCAAATCTTCTTACTGTTACCATTGG‑TTTTTGCAGCGGTACCGCAG‑TTTCGACGA‑CTGTTTTCACCATTTATTG‑TTATCA‑GGTGTGTGCTCGCGCACG‑CATTA‑TGGGCGGCCTGTGCTCTCTGCTggcacagc                                                                                     >  1:360268/2‑233 (MQ=255)
                                                                                                                                                       ggTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGGGACCAAAAACTTCCT‑CCAAACCATGCTGATTTTGTTGATTGCCGCCCCGCTGGCAATCCTGCTGATCTGCCCGATTGG‑TGTCTGGGTCGGTAGCGCCG‑TTTCGCTGC‑CGGGTTATACCCTTCCTGG‑TTATTT‑GGGCGG‑GCTTTCAGCCGCCCTTC‑AGGGCCGGCTGGGGGGCCTGCGGGGAATAACgggg                                                                                >  1:577361/1‑239 (MQ=255)
                                                                                                                                                         ggccttgaaaatatAGAACGCTGGGGGGAAAGGTTTCCCCCGGCGGGAACCAAAAACTTCCT‑CAACCCAAGGCTGATTGTGTTGATTGGCGGACCGCTGGCAATCCGTCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTGCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTGTGCTGGTAATGACCGGGa                                                                               <  1:752704/228‑1 (MQ=255)
                                                                                                                                                               tCATATATCGAACGCTGGGTGGACAGCATTACCCCGGCGGTGACCAAAAGCTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTctg                                                                                               <  2:282293/216‑1 (MQ=255)
                                                                                                                                                               tCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAAt                                                                                       <  2:162893/224‑1 (MQ=255)
                                                                                                                                                                tttatttctacccttGGGTGGCAAGCATTCCGCCGGCGGTTACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCGGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACcg                                                                        <  1:343264/224‑1 (MQ=255)
                                                                                                                                                                aaaataaCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGGTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACcg                                                                        <  2:1071033/235‑1 (MQ=255)
                                                                                                                                                                     atCGAACGCTGGGTGGACAGCATTACCCCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCtt                                                                  <  2:129856/239‑1 (MQ=255)
                                                                                                                                                                         aaCGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTTTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGGTTGG‑TATCTGGATCGGTAGCGCCC‑TTTCTGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑TGGCTG‑GCTTTCAGTCGCCATTT‑TGGGCGCGCTGTGGCCCCTGCTGGTAATAACCGGGGTGTCCCGCGTCttttg                                                               >  2:119283/1‑236 (MQ=255)
                                                                                                                                                                         aaCGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTc                                                                                                                           <  2:1279146/179‑1 (MQ=255)
                                                                                                                                                                           cgtggggTGGCCAGCATTACCCCGTCGGTGACCAAACACGTCCT‑CAAACCAATGCGGATTGTGTGGATTGCCGTACCTGTGGCAGTCCGGGTGATGGGCCCGGTTGG‑TAACTGGCTCGGTAGCGCCA‑TTTTGGTGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCGG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCGTTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGcc                                                            <  1:1035090/235‑1 (MQ=255)
                                                                                                                                                                           cgcgGGGTGGACAGCATTACCCCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGcc                                                            <  2:529713/239‑1 (MQ=255)
                                                                                                                                                                               gggTGGACAGCATTCCGCCGGGGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAAcc                                                        <  2:635455/239‑1 (MQ=255)
                                                                                                                                                                                 ggtgACAGTATTACCGCGGCGGGGACCACAATCTCCCT‑CAAACCAATGGTGATTTTGTGGATTGCCGCTCCGCGGGCAATCCTGCTGATTGGCCCGATTGG‑TATGTGGGTCGGTAGCGCCC‑TTTCGGCGC‑TGGTTTACACCATGCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACca                                                       <  1:96728/235‑1 (MQ=255)
                                                                                                                                                                                 gTGGACAGCATTTCACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACca                                                       <  2:8760/238‑1 (MQ=255)
                                                                                                                                                                                 gTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACca                                                       <  2:1182250/238‑1 (MQ=255)
                                                                                                                                                                                         cATTACCCCGGCGGTGACCAAAACCGTCCT‑CAAACCAAGGCTGATTGGGTTGATTGCCGCCCCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACGCCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACcg                                                                        <  2:1102593/213‑1 (MQ=255)
                                                                                                                                                                                         cATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATc                                                                                                                                              <  2:840315/146‑1 (MQ=255)
                                                                                                                                                                                                   gcggTGACCAAAAACTTCCT‑CAAAACAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGGCCGGTTGG‑GATCTTGATCCGTAGCGCCA‑TTTTGGCGC‑CGGTTTTCACCATTCATGG‑GTATCT‑Gggctg‑gct                                                                                                                                  >  1:840315/1‑146 (MQ=255)
                                                                                                                                                                                                              aaaCTTCCT‑CAAAACCATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGTCCCGATTGG‑GATCTGGATCGGTAGCGCCA‑TTTCGCAGC‑TGGTTTACACCATTCATGG‑TTATTC‑GGGATG‑CGTTTCAGTCGGCATTA‑TTGGCTCTCTGTGTCCTCTGCTGGTAATGCCCGGGTGGCACCTCGCCTTTTCGCCACCCACCCTTCAAACCATTTGCGCAAACGGCCaag                         >  1:647879/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCTTCATCCAGACCATTGCCGAAACCGGCAAAGAAGGGa                   >  2:667271/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGa                   >  2:1053537/1‑239 (MQ=255)
                                                                                                                                                                                                                     cT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGTTAGCGCCA‑TTTCGgcg                                                                                                                                                                         >  2:123810/1‑93 (MQ=255)
                                                                                                                                                                                                                     cT‑CAAACCAATGATGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGgcg                                                                                                                                                                         <  1:123810/93‑1 (MQ=255)
                                                                                                                                                                                                                              aaTGGTGATTTTGGTTATTGCCGCACCGGTGGCAATCCTGCTGATTTGCCCGATTGG‑TATCTGGATCGGGAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGc          <  2:168761/239‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGGAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGCCGTCaga   >  2:1232290/1‑239 (MQ=255)
                                                                                                                                                                                                                                      ttGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCCAAGAAGGGATGGTCAGGCCGCCagag  >  2:805947/1‑239 (MQ=255)
                                                                                                                                                                                                                                       tgtgTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑ATCCGGCGC‑TGCTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTGTCAGTCGCCATTC‑TGGCCGCGCTTTGGCCTCTGCTGGTAATGACCGGGATGTCACGCGTCTTTAAGCCACCCATCATTCAGCCCATTGCACAAACCGGCAaagaaagagaggaaaggacgacaagg  >  1:575932/1‑217 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
CTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGCGGCCCAGGGTGAACATGTCGAATTTGCCCTGATTCCGGTCACCGCGGTGAAATACACCTACACGGTGATCCCGGCGCTGGTGATGACCTGGTGCCTGTCATATATCGAACGCTGGGTGGACAGCATTACACCGGCGGTGACCAAAAACTTCCT‑CAAACCAATGCTGATTGTGTTGATTGCCGCACCGCTGGCAATCCTGCTGATTGGCCCGATTGG‑TATCTGGATCGGTAGCGCCA‑TTTCGGCGC‑TGGTTTACACCATTCATGG‑TTATCT‑GGGCTG‑GCTTTCAGTCGCCATTA‑TGGGCGCGCTGTGGCCTCTGCTGGTAATGACCGGGATGCACCGCGTCTTTACGCCAACCATCATTCAGACCATTGCCGAAACCGGCAAAGAAGGGATGGTCATGCCGTCAGAG  >  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: ‑