Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,533 C→T Y213Y (TAC→TAT)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,5330CT96.8% 166.7 / ‑5.6 64Y213Y (TAC→TAT) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (30/31);  minor base G (2/0);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

AACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCC  >  NC_000913/4625289‑4625767
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
aaCTGGCCGGAACGGGCGAGATGGTGGCGGAAGAAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGCAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGATTTCCTTTCTTTGCTGAATATCTgagc                                                                                                                                                                                                                                        >  2:457544/1‑248 (MQ=255)
 aCTGGCCGGAACGGGGGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCAAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGcgcg                                                                                                                                                                                                                                       <  1:109859/251‑1 (MQ=255)
    ggCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGAc                                                                                                                                                                                                                                     <  1:93958/250‑1 (MQ=255)
     gCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGAGTGACGACGCTGAAAGGCGCTGACGCCAATATTCAGCAACAGGCGCGTGAACATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCACGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGcgcgtca                                                                                                                                                                                                                                    >  1:409534/1‑247 (MQ=255)
       cggcaCGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTTCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACa                                                                                                                                                                                                                                    >  2:5985/5‑248 (MQ=255)
       cGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAg                                                                                                                                                                                                                                  <  1:5985/250‑1 (MQ=255)
        ggAACAGGCGAGATGGTAGCGGAAGTAACAGAAAGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCACAGACGCGTGGCGACGCTGAAAGGCGCTGGCGCCAATATTCTGAAACAGGTGCGTGAAAATATGCCGCTGATGACAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAATGCTGGGCTGGAAAGAGGCGCTGACCACCAGCGGATTTATTTTCTTTGCTGAATATCAGAGCGACAAgtt                                                                                                                                                                                                                                >  2:433954/1‑249 (MQ=255)
            cGGGCGAGATGGTGGCGGAATTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGTCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTgcgc                                                                                                                                                                                                                            <  1:65147/251‑1 (MQ=255)
              ggCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCt                                                                                                                                                                                                                          >  1:552770/1‑251 (MQ=255)
                cGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAACGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGAGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATTTCTGCGGGCCAAGTTGCGcccgc                                                                                                                                                                                                                        >  2:225306/1‑248 (MQ=255)
                  agaTGGAGGCGGAAGGAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCCGCCTGCGCAGCCGCGTGGCGACGCTGACAGGCGCTGACGCCAATATTCAGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCGGGAAAGTGGCGCATGCCTCCGGCAGCTTTCCTTTCTTTGATGAATATCTGcgcggcgattttggccttacc                                                                                                                                                                                                                      >  2:544623/1‑234 (MQ=255)
                          gCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGACGCGAGGCGACGCTGAAAGGCGCTGACGCCACTATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGTATGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCGGCGCGCCAAGTTGCGCCTGACCGACGTgg                                                                                                                                                                                                               >  2:526336/1‑250 (MQ=255)
                              aaGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGAGCGACAAGTTGCGCCTGACCGCCGTGGTAGcc                                                                                                                                                                                                          >  1:70454/1‑251 (MQ=255)
                                   aCCGAACGGGCGATGCGCGGCGAACTCGAGTTTACCGCCAGCAAGCGCCACCGCGAGGCGCCGCTGACAGGCGCGGACGACAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGACACGCTGGGCGGGAAAGTGGCGATTGCCTACGGCGGCTTTACTTTCTTTGCTGAATAGCTGCGCGACAAGTTGCGCCTGACAGCCGTGGTGGCCAATgg                                                                                                                                                                                                     >  2:68265/1‑250 (MQ=255)
                                       aacggtcgaggcgggtcgAACTCGGTTTCACCGCCAGCCTGCGCAGCAGCGTGGCGACGCTGAAAGGCTATGACGCCAATATTCTGCAACAGGGTCGTTAAAATCTGCCGCTGATGCCAGGCTTAACCCAACTGGTGCTCAAGCT‑GGAAACGCTGGTCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCGGACCGCCGTGGTAGCCAATGAACTg                                                                                                                                                                                                 <  2:519704/235‑1 (MQ=255)
                                          gggCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTgg                                                                                                                                                                                                >  1:588152/1‑249 (MQ=255)
                                                     ggcgAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACgg                                                                                                                                                                                   >  1:357561/1‑251 (MQ=255)
                                                        aaaCTCGTTTTTCCCGCCAGCCTGCGCAGCCGCGTGGCGCCGCTGAAAGGCGCTGACGCAAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTa                                                                                                                                                                                 <  2:357561/249‑1 (MQ=255)
                                                                    cccGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCaa                                                                                                                                                                    <  1:544003/250‑1 (MQ=255)
                                                                       cccacccttcgcAGCCGCGGGGCGACGCTGAAAGGCGTGGACGCCAATATTCTGCAAAAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAtgt                                                                                                                                                                 <  2:202756/242‑1 (MQ=255)
                                                                          aGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGAt                                                                                                                                                              >  1:317025/1‑251 (MQ=255)
                                                                           gCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATc                                                                                                                                                             <  1:656269/251‑1 (MQ=255)
                                                                              tGCGCAGCCTCGTGGCGCAGCTGAAAGGCGCTGACGCCCATATTCTGCAAAAGGTGCGTGAAATTCTGCTGCTGATGACAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGCGATCATGGACGGTCACTTTACAGGCAATGTGAtaggg                                                                                                                                                          >  2:463642/1‑247 (MQ=255)
                                                                               gcggaGCCGCGTGGCGACGCTGAAAGGCGCTGACTCCAATATGCGTCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTTCGCCTGACCGCTGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGc                                                                                                                                                          <  1:382543/246‑1 (MQ=255)
                                                                               gcgcAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCgg                                                                                                                                                           <  1:567451/249‑1 (MQ=255)
                                                                                cgcAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCg                                                                                                                                                         >  1:335640/1‑250 (MQ=255)
                                                                                      cgcgtgGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGt                                                                                                                                                  >  1:386313/4‑251 (MQ=255)
                                                                                      cgcgtgGCGACGCTGAAAGGCGCTGACGCAAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAACCGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGt                                                                                                                                                  <  1:56611/248‑1 (MQ=255)
                                                                                         gtgGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGACAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCAAGGACGGTAAATTTACCGGCAATGAGGTCGGCGACATCGt                                                                                                                                                  >  2:397273/1‑248 (MQ=255)
                                                                                          tgGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAc                                                                                                                                              >  1:497787/1‑251 (MQ=255)
                                                                                          gggCGACGCTGAAAGGCACTTACGCCAATATTCGGCAACAGGTGCGTAAAAATCTGCCGCGGATGCCCGGCTTAACGCAACTGGTGCGCAAGCT‑GAAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTGTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGTCCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAc                                                                                                                                              <  2:386313/250‑1 (MQ=255)
                                                                                                    gAAAGGCGCTGACTCCAATATTCGGCAACAGATGCGTGAAAATCTGCCGCGGATGCCAGGATTAACGCCACTTGTGCTCAAGCTTGGATA‑GCTGGTCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTTCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCCTGGACGGTCAATTTACCAGCAATGTGATCGGCGACATCGTACACTCGCCGTAc                                                                                                                                     >  2:175108/1‑250 (MQ=255)
                                                                                                    gAAAGGCCCTGGCGCCAATATTCTGCAGCAGGGGCGGTAAATTCGGCCGCTGAGGCCAGGCTTAACGCAACTGGTGGCCAAGCT‑GGAAATGCTGGGTTGGAAAGTGGGGATTGCCTCGGGCGGCTTTACTTTCTTGTCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACa                                                                                                                                    <  2:438685/251‑1 (MQ=255)
                                                                                                                 gCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGCCGCGCAGTACAAAGCGCAAActct                                                                                                                       >  2:128427/1‑251 (MQ=255)
                                                                                                                   cAATATTCTGCAACAGGCGCGTGACAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCATCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGAAAAGTTGCGCCTGGCCGCAGTGGTAGCCAATGCAATGGAGATCATGGACGGTAAATTTACCGGCAATGGGATCGGACACATCGTAAACGCGCAGTACAAAGCGGAAActctt                                                                                                                      >  2:531781/1‑249 (MQ=255)
                                                                                                                     atatTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAc‑t                                                                                                                  <  1:141575/251‑1 (MQ=255)
                                                                                                                       atTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAcagg                                                                                                                 >  2:503747/1‑248 (MQ=255)
                                                                                                                       atTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAc                                                                                                                    <  1:503747/248‑1 (MQ=255)
                                                                                                                            gCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGTGCAGTACAAAGCGAAAACTCTGAc‑tcgcct                                                                                                             <  1:245472/249‑1 (MQ=255)
                                                                                                                                      cGTGAAAAGCCGCCGCTGATGCCAGGCTTAACGCAACGGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTa                                                                                                <  2:306397/251‑1 (MQ=255)
                                                                                                                                          aaaaTCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGaa                                                                                            <  1:457544/251‑1 (MQ=255)
                                                                                                                                            aaTCGGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAAt                                                                                          <  1:305300/251‑1 (MQ=255)
                                                                                                                                                            gaGGCTTATCGAAACGGGGGCTCAGGCT‑GGAGACGCTGGGTTGGAAAGGGGCGATGGCGGCCGGCGGCTTTACGTTCTTTGCTGAAGATCTGCGCGACAAGTGGGGCCTGCCGGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTTAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGTAGAcc                                                                          <  2:640237/250‑1 (MQ=255)
                                                                                                                                                            cAGGCTTAACGCAACTGGAGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGAGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGAAGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGAcc                                                                          >  1:328748/1‑251 (MQ=255)
                                                                                                                                                             aGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCg                                                                         <  1:225306/251‑1 (MQ=255)
                                                                                                                                                             aGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCg                                                                         >  1:258777/1‑251 (MQ=255)
                                                                                                                                                               gctataCGCATCTGGTGCTCAATCT‑GGAACCGCGGGTTTCGCAAGTGGCGTTTGCCTCTGGCGGCTTTACTTTCTTGGGTGAATATCGGCGCTAAAAGTTGCGCCGGACCTCCGTGGTAGCCAATGAACTGGAGATTATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCCCGTGC‑AGTAGTATGAAATCCCGCTGGCGCAGACCGtg                                                                       <  2:414571/246‑1 (MQ=255)
                                                                                                                                                               gCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCCGACCGtg                                                                       >  2:343247/1‑251 (MQ=255)
                                                                                                                                                                  taccGCAACTGGGTCTCAAGCT‑GGAACCGCTGGGCTGGAAAGTGGCGCTGGCCTCCGGCGTCTTTACTTTCTTGGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑GCGCCTCGCGC‑AGGAGTATGAAAGCCCGCTGGCGCTGACCGtga                                                                      <  2:65851/246‑2 (MQ=255)
                                                                                                                                                                     cGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCGGCGCGACAAGTTGCGCCGGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGGAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgatt                                                                 <  1:93129/251‑1 (MQ=255)
                                                                                                                                                                          cTGGAGCTCACGCA‑GGAAACGCAGGACTGGTAAGGGGCGATTGCATACGGCGGCTTTACTTTCTTTGCTGAATAGCTGCGCGACAATTTGCGCCTGACCGCCGTGGTAACCAATGCACTGGAGATCATGGACGGTAAATTTACCGGCAATGAGACCGGCGACATGGTCGACGCGAAGAACAAAGCGAAAACTCTGAC‑TCGCCCAGCGA‑AGGTGTATGAATTCCCGCGGGCGCAAGCCGtggagcttggcgc                                                            >  2:441687/1‑250 (MQ=255)
                                                                                                                                                                             ggtgttcaggtgtgaactgttgggctggAAGGTGCGGATTGCCTCCGGCTGCGTGCCTTTCGTTGCTGAAGATCGGCGGTACAAGTTGCGCCGGACCGCCGTGGTAGCCAATGAACGGGAGCTCTGGGACGGTAAATTTACCGGCAATGTGATCTGCGACATCGTAGACGCGCAGTACAAAGCGAAATCTCTTAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattggccagg                                                          <  2:148680/231‑5 (MQ=255)
                                                                                                                                                                               gCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGa                                                        <  1:397273/250‑1 (MQ=255)
                                                                                                                                                                                         ggAAACGCTGGGCTGGAAAGTGGCGATTTCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGAc                                               <  1:167454/250‑1 (MQ=255)
                                                                                                                                                                                            aaCGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGc                                           >  1:81740/1‑251 (MQ=255)
                                                                                                                                                                                                  ggggtggAAAGTGGCGATGGCCACGGGCGGCTTTACTTTCTGTGCGGAATATCTGCGCGACAAGTGGCGCTTGACCGCCGTGGTAGCCAATGAACGGGAGATCAGGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCAGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCgatga                                     <  1:331688/247‑1 (MQ=255)
                                                                                                                                                                                                      tggAAAGTGTCGATTACCTCCGTCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGCGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAACACTCTGAC‑TCGCCTCGCGC‑CGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTAGCCATTGAGCCACTGACCTGCCACTGATCa                                  >  2:117420/1‑250 (MQ=255)
                                                                                                                                                                                                       ggAAAGTGGCGATTTCCTCCGTCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCaa                                 <  1:450844/250‑1 (MQ=255)
                                                                                                                                                                                                       ggAAAGTGGCGATTGCCTCCGGCGGCTTTCTTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCaaa                                <  1:544623/251‑1 (MQ=255)
                                                                                                                                                                                                        gAAAGTAGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGATGCGCATAACCGGCGTGGTAGCCATTGAACTGGAGATCATGGACGGAAACTTTACCGGCAAAGTGATCGGCGACATCGTAGAAGCGCAGTACAACGCGAAACCTC‑GACCTCGCCTCGCGCTAGCA‑CATGAAGTCCCGCTGGCGAAGACCGTGGCGATTGGCGATGGGGCCAATGACGTGcgcgtgcgcaa                                 >  2:496348/1‑239 (MQ=255)
                                                                                                                                                                                                                 aaTTGCCCCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAATTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGc                      <  1:369990/250‑1 (MQ=255)
                                                                                                                                                                                                                    tGCCTCCGACGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGGCGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGCCGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGCGCCCATGACAGGCCGAGGCTCAAAGCGGCCGGGcggg                   >  2:365879/1‑248 (MQ=255)
                                                                                                                                                                                                                               gcTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTAc        >  1:170007/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tttctttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgcc  >  1:197511/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
AACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCT‑GGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAC‑TCGCCTCGCGC‑AGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCC  >  NC_000913/4625289‑4625767

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