Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,188,087 T→C T110A (ACG→GCG)  yehA ← putative fimbrial adhesin YehA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,188,0870TC98.2% 190.6 / NA 56T110A (ACG→GCG) yehAputative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (28/27);  minor base A (1/0);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTG  >  NC_000913/2187845‑2188330
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAAt                                                                                                                                                                                                                                              >  2:188743/1‑251 (MQ=255)
 ttGTATGACAATGTCAGATTAAAGGAATAGACACCTGCACTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCCACTGAAt                                                                                                                                                                                                                                              >  1:638453/1‑250 (MQ=255)
               cAGATTTAAGGAATAGACAACTGCTTTTTGGGCCATAGAATTTCTGAGGCTATGTACACATGCATCCCATGACTTGCCCCCTAATAAAAACAGAAGAAATTGTACTGCGGCCCTCCCGCGTGTGGATGATGAAGCACTGGATAATGTCTCAACCCCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGaaa                                                                                                                                                                                                                                <  2:858471/250‑1 (MQ=255)
                       aGGAAAAGACCCCTGCACTTTGTGCCAATGCATTTCTGTAGCTATATACACATGCACCCCATGACTTGCCACCTAAAAAACAAAGTAGAAATTGTACTGCGGCTCTCCCGCGGGTCTATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGa                                                                                                                                                                                                                       <  2:638453/251‑1 (MQ=255)
                                        ttttGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGAt                                                                                                                                                                                                     >  1:104590/1‑251 (MQ=255)
                                                    gCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTAtt                                                                                                                                                                                         <  1:195112/251‑1 (MQ=255)
                                                    gCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTAtt                                                                                                                                                                                         <  1:671969/251‑1 (MQ=255)
                                                        ttCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACAGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTc                                                                                                                                                                                     >  1:418810/1‑251 (MQ=255)
                                                             tAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTc                                                                                                                                                                                >  1:586914/1‑251 (MQ=255)
                                                               gCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCaa                                                                                                                                                                              <  1:369058/251‑1 (MQ=255)
                                                                  aTTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAAt                                                                                                                                                                           <  1:247334/251‑1 (MQ=255)
                                                                   tttACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGTTCCATTCAAAATg                                                                                                                                                                          >  2:15977/1‑251 (MQ=255)
                                                                      acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTa                                                                                                                                                                        <  1:166796/250‑1 (MQ=255)
                                                                         caTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaatt                                                                                                                                                                    <  2:331823/251‑1 (MQ=255)
                                                                           tGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACa‑a                                                                                                                                                                                                           >  1:361783/1‑210 (MQ=255)
                                                                           tGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACa‑a                                                                                                                                                                                                           <  2:361783/210‑1 (MQ=255)
                                                                            gCATCCCAGGACTTGCCACCTAATAAACACCGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaatttaa                                                                                                                                                                 <  1:16779/251‑1 (MQ=255)
                                                                                     gACTTGCCAACTAAAAAACAAAGTAGAAATTGTAATGCGGCACCCCCGCGTGTCGATGATGAAGCACTGGATAAGGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATGTTATgg                                                                                                                                                         <  2:586914/250‑1 (MQ=255)
                                                                                          gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGCAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAACCTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGt                                                                                                                                                   >  1:931825/1‑251 (MQ=255)
                                                                                              ccTAATAAACAAAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATACCCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTtat                                                                                                                                                >  1:919866/1‑250 (MQ=255)
                                                                                                  ataaACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCACCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACTTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATAGAATCGTTATAAAGt                                                                                                                                           >  2:11740/1‑251 (MQ=255)
                                                                                                       cacaGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCcaac                                                                                                                                      <  1:188743/251‑1 (MQ=255)
                                                                                                              gAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGACTGATTATTTGCCA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACaataat                                                                                                                                >  1:882492/1‑250 (MQ=255)
                                                                                                               aaaTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTAAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAatat                                                                                                                              <  2:936767/251‑1 (MQ=255)
                                                                                                               aaaTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGAAAATGGCGCAACACCATTAACAAGCAATGAACCGCCGCTCCCCGCAGAAACATTCACTTTTGCATTACTGGCTGCAGTTTTCACCGCAACTGAATAACCACTTGTGAAACTTGTCTGATTATGTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGCATCGTTATAAAGTCCAACAATAatat                                                                                                                              <  2:332772/251‑1 (MQ=255)
                                                                                                                   tgttcTGCGGCACTCCCGCGTGTCGATGATGATGCCCTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCACACCGGCGAAACATGGACTTTTGCTTTCCTGGCAGCAGGTGTCACCCCAACTTAATAACAACTGGGGAAACGGGTCTGATTATTTG‑CACACGTGCGATTGTTTTTATTGGTCCATTCAACATATCAATTTAATATTTGGGTATCTTTTTAAAGACCAACaataatcaccgc                                                                                                                          >  2:668016/5‑245 (MQ=255)
                                                                                                                   tGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAATTCCAACAATAATATCCg                                                                                                                           <  1:827164/250‑1 (MQ=255)
                                                                                                                   tGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCg                                                                                                                           >  1:548774/1‑250 (MQ=255)
                                                                                                                    gatgttcGGCACTCCCGCGTTTCGATGATAAAACCCGGGATAATGTCTCATCACCATTAATTATCACTGAACCGCCGCCCCCCGCAGAAACATTCACGTTTTCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGAc                                                                                                                         <  2:192450/245‑1 (MQ=255)
                                                                                                                     aaCTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGGGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGCTTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATAGTTATAAAGTCCAACAATAATATCCGACg                                                                                                                        >  2:742850/2‑251 (MQ=255)
                                                                                                                                       tgtCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGAt                                                                                                        >  1:455464/1‑249 (MQ=255)
                                                                                                                                               atgaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCACCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGATATCGTTATAAAGTCCAACAATAATATCCGACAGGCTCATAGGGGAGATAGTATTTGt                                                                                               >  1:393767/1‑250 (MQ=255)
                                                                                                                                                tgaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTAc                                                                                             >  1:551437/1‑251 (MQ=255)
                                                                                                                                                 gaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTAc                                                                                             <  1:4448/250‑1 (MQ=255)
                                                                                                                                                        tGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGATCATAGGGGGGATAGTATTTGTACTTACGCAt                                                                                     >  1:949392/1‑251 (MQ=255)
                                                                                                                                                               tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAAt                                                                              >  2:969926/1‑251 (MQ=255)
                                                                                                                                                                ggccaaaaacCATTAATAATCACTGAACCGCCGCCCCCCGAAGAAACATTCCCTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGTTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAAGTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGTCGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATc                                                                             <  2:586877/243‑1 (MQ=255)
                                                                                                                                                                  cGCAACCCCATTAAACATCACGGAACCGCCGCCCCCCGCAGAAACATACACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATGTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCAt                                                                           <  2:104590/251‑1 (MQ=255)
                                                                                                                                                                              aaTCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGc                                                                <  1:969926/250‑1 (MQ=255)
                                                                                                                                                                               atcatcACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTc                                                              <  1:752/251‑1 (MQ=255)
                                                                                                                                                                                 catcaCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGCTTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTCAAATTATGCTCat                                                            >  2:978462/1‑251 (MQ=255)
                                                                                                                                                                                    caCTGAACCGCCGCCCCCCGCAGAAACATTCACTGTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATTACGAAGCATTAAAATTATGCTCATATg                                                         <  2:903413/251‑1 (MQ=255)
                                                                                                                                                                                         aaCCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGcc                                                    >  2:102451/1‑251 (MQ=255)
                                                                                                                                                                                           cccccgccCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATGTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAt                                                  <  2:455464/248‑1 (MQ=255)
                                                                                                                                                                                            cgccgcCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAtt                                                 >  1:281716/1‑251 (MQ=255)
                                                                                                                                                                                                       gCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGt                                      >  2:190229/1‑251 (MQ=255)
                                                                                                                                                                                                        cAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTa                                     <  1:455793/251‑1 (MQ=255)
                                                                                                                                                                                                         aGAAACATTCCCTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTaa                                    <  1:227993/251‑1 (MQ=255)
                                                                                                                                                                                                                        tgccttACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCt                     >  1:810468/5‑251 (MQ=255)
                                                                                                                                                                                                                        tGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCt                     >  1:810480/1‑251 (MQ=255)
                                                                                                                                                                                                                          cATTACTGTCCGCAGGTGTCACCGACACTGAATAACCACTGGTGAAACTGGTCTGATTCTTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGCCGGGCTCATAGGGGTGATAGGATTTGTACTTACGCATGACGACTCATTAAAATTATGCTCTTATGTCGCCATTTGTGTTTGCGTAAGTGCAAAATCCCCCTg                    >  2:440361/1‑250 (MQ=255)
                                                                                                                                                                                                                             aaCTGGCCGCAGGTGACACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGaaa                 >  2:39596/2‑249 (MQ=255)
                                                                                                                                                                                                                                    cGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCa                                                                                                                  >  1:349418/1‑146 (MQ=255)
                                                                                                                                                                                                                                    cGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCa                                                                                                                  <  2:349418/146‑1 (MQ=255)
                                                                                                                                                                                                                                           gtCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTg  >  1:452803/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ggcACCGCAACTGAATAACCACTGGTGAAACTGGTCTGACTATTTGACA‑ACTTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGGAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTg  <  2:548774/249‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACA‑ACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTG  >  NC_000913/2187845‑2188330

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