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,5330CT98.2% 167.1 / NA 58Y213Y (TAC→TATserBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (31/25);  minor base G (0/1);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAA  >  NC_000913/4625289‑4625774
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
aaCTGGCGGGAACGGGCGAGTTGGGGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCGGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTgcgc                                                                                                                                                                                                                                             <  1:446502/251‑1 (MQ=255)
             gggCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTCATTTCTTTGCTGGATATCTGCGCGCCAAGTTGCGcc                                                                                                                                                                                                                                >  2:386321/1‑251 (MQ=255)
                        gggcgaaataaacggacgggggggggcgcGGCAAACTCGATTTTACCGCCAGCCTGCGCCGCCGCGTGGCGACGCTGAAAGGCGCTGACGCAAATATCCTGCAACGGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGCCCGCCGTg                                                                                                                                                                                                                     <  2:390732/226‑1 (MQ=255)
                           cGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGt                                                                                                                                                                                                                   >  1:325151/1‑250 (MQ=255)
                            ggAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGttg                                                                                                                                                                                                                 >  2:400614/1‑249 (MQ=255)
                                          gggCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTgg                                                                                                                                                                                                     >  1:366724/1‑249 (MQ=255)
                                            gCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATAGCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGAt                                                                                                                                                                                                 <  1:391503/251‑1 (MQ=255)
                                                    cggcgAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACg                                                                                                                                                                                         >  1:82008/1‑251 (MQ=255)
                                                       ggaactggtttttCCCGCCAGCCTGCGCAGCGGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTAAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACGTTCGTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCAGGGACGGTa                                                                                                                                                                                      <  2:247665/242‑1 (MQ=255)
                                                          atcctaTTTTCCCGCCAGCCGGCGCAGCCGCGGGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAt                                                                                                                                                                                   <  2:133707/246‑1 (MQ=255)
                                                              gATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAc                                                                                                                                                                               >  1:93001/1‑251 (MQ=255)
                                                                  ttcccGCCAGCCTGCGCAGCCGCGTGGCGCAGCGGAAAGGCGCTGACGCCAATATCCTGCAACAGTTGCTTGAAAACCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGTTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATTTCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCgg                                                                                                                                                                            <  2:468331/247‑1 (MQ=255)
                                                                  ttACCGCCAGCCTGAGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTACCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGACAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGATCTGGAGATCATGGACGGTAAATTTACCGGc                                                                                                                                                                           >  1:19903/1‑251 (MQ=255)
                                                                   tACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGc                                                                                                                                                                           >  2:251398/1‑250 (MQ=255)
                                                                     ccGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAt                                                                                                                                                                        >  1:386981/1‑251 (MQ=255)
                                                                           gCCTGCGAGGCCGCGTGGCGATGCTAAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCTGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGATGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATc                                                                                                                                                                  <  1:189948/251‑1 (MQ=255)
                                                                            ccTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCg                                                                                                                                                                 <  1:251398/251‑1 (MQ=255)
                                                                             atGCGCAGCAGCGTGGCGACGCTGAAAGACGCAGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAtg                                                                                                                                                                       >  1:367517/2‑244 (MQ=255)
                                                                                             ggACGCGGAAAGGGCCTGACGCCAATCTTCGGCAACAGTGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGTTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTTCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAcgcg                                                                                                                                                <  2:346208/250‑1 (MQ=255)
                                                                                               aCGCTGAAAGGCGCTGCCGCCAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTAAATATCTGCGCGACAAGTTGAGCCTGACCGCCGTTGGAGCCAATTAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGCTCGGCGGCATCGTAGAcgccc                                                                                                                                               >  2:240315/1‑248 (MQ=255)
                                                                                                                   cAATATCCGGCAAAAGGGGCGTGAAAAGCTGCCGATGATGCCAGGCTTAACGCAACTGGGTCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACGTTCTTTGCTGAATATCTGCGCGCCAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGAACGGCGACAGCGTAGACGCGCAGTACAAAGCGAAAActccga                                                                                                                          <  2:19903/251‑4 (MQ=255)
                                                                                                                     atatTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAct                                                                                                                        <  1:59292/251‑1 (MQ=255)
                                                                                                                      tatTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGActc                                                                                                                       >  1:413668/1‑251 (MQ=255)
                                                                                                                      tatTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGActc                                                                                                                       >  1:355792/1‑251 (MQ=255)
                                                                                                                             caacagtggCGTGAAAATTCGCCGCTGTTGCCAGGCTTACCCCAACTGGTGTCCATGCTGGAAACGCTGGTCTGGAAAGTGGAGATTGCCTCCGGCGGCTTTACTTTTTTTGCTGATTATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACGCTGActcgcctc                                                                                                                  <  2:354239/241‑1 (MQ=255)
                                                                                                                             cAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAcc                                                                                                                                                                              <  1:366753/189‑1 (MQ=255)
                                                                                                                             cAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAcc                                                                                                                                                                              >  2:366753/1‑189 (MQ=255)
                                                                                                                             cAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGAAGACGCGCAGTCCAAAGCGAAAACTCTGActcgcctcgc                                                                                                                >  2:122049/1‑251 (MQ=255)
                                                                                                                               aCAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCgcgc                                                                                                              >  1:220517/1‑251 (MQ=255)
                                                                                                                                 aGGTGCGTGAAATTCGGCCGCGGATGCCAGGCTTAACGCAACTGGGGCTCAAGCTGGAAACGCTGGGCTGGAAAGGGGCGATTGCCTCCGGCGGCTTTACGTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAg                                                                                                            <  2:366724/251‑1 (MQ=255)
                                                                                                                                 aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTa                                                                                                                                                      <  2:240424/209‑1 (MQ=255)
                                                                                                                                 aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTa                                                                                                                                                      >  1:240424/1‑209 (MQ=255)
                                                                                                                                 aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCgcgc                                                                                                              >  2:127366/1‑249 (MQ=255)
                                                                                                                                      cGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGt                                                                                                        >  1:156500/1‑250 (MQ=255)
                                                                                                                                       gTGAAAATCTGCTGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGCCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTAt                                                                                                      <  1:319566/251‑1 (MQ=255)
                                                                                                                                       gTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTAt                                                                                                      <  1:439629/251‑1 (MQ=255)
                                                                                                                                         aaaaaTCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAATTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACTCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATg                                                                                                     <  1:147676/249‑1 (MQ=255)
                                                                                                                                               cTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTCCAAAGCGAAAACTCTGACTCGCCTCGCGCCGGCGTATGAAATccc                                                                                              >  2:205953/1‑251 (MQ=255)
                                                                                                                                               cTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                              >  1:418433/1‑251 (MQ=255)
                                                                                                                                                 gCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                              <  2:141371/249‑1 (MQ=255)
                                                                                                                                                 gCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATc                                                                                                <  1:127366/247‑1 (MQ=255)
                                                                                                                                                  ccGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTACGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTACATTTACCGGAAATGTGATCGGCGACCTCGAAGACGCGCGGTACAAAGCGAAAAATCTGACTCGCCTCGCGCCGGAGTATGAAATCCCGCt                                                                                           >  2:465212/1‑251 (MQ=255)
                                                                                                                                                            cAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGGGTCTGAAATCCCGCTGGCGCAGAcc                                                                                 >  2:155909/1‑251 (MQ=255)
                                                                                                                                                                cTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtgg                                                                             >  1:463132/1‑251 (MQ=255)
                                                                                                                                                                        aaCTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattggc                                                                     <  1:205226/251‑1 (MQ=255)
                                                                                                                                                                             gTGCTCAAGCTAGAAACGCTGGGCTGGAAAGTGTCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCACTGAACTGGAGATCATGGACGGAAAATTTACCGGCACTGAGATCGGCGACATCGTCAACGCGCCGTACAAAGCGGAGACTCTGCCTCGCCTCACGCAGGAGTATGACATCCCGCGGGCGCAGCCCAgggcgagtggcgaacg                                                                >  2:87470/1‑248 (MQ=255)
                                                                                                                                                                                     gCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATg                                                        <  1:400614/251‑1 (MQ=255)
                                                                                                                                                                                      cTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGTCGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTAAAAAGCGAAAACTCTGACTCGCCCCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGAGATGGCGCCAATg                                                        >  2:84266/1‑250 (MQ=255)
                                                                                                                                                                                                              gcgatgtcctccgtcggtgttattttctgtGCTGAATATCGGCGCGACAAGTTGCGCTTGGCCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGTTAAATTTACCGGCAATGTGATCGGCGACATCTTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAgg                               <  1:468337/228‑1 (MQ=255)
                                                                                                                                                                                                               cGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAGTATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAggg                              >  1:131912/1‑251 (MQ=255)
                                                                                                                                                                                                                    gCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAACGCGGCAGGGCTggg                         >  1:476399/1‑251 (MQ=255)
                                                                                                                                                                                                                       tCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGAt                      >  1:29197/1‑251 (MQ=255)
                                                                                                                                                                                                                         cggcggCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGGTTg                    >  1:286884/1‑251 (MQ=255)
                                                                                                                                                                                                                            cggcTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCt                 <  1:122049/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tatttgttttGCTAAATATTTGCGCGAAAAGTTTCGCTTGACCGCCGGGGTAGCCAAGGAACTGGAGATCATGGCCGGTAAATTTACCGGCAATGTGACCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACCCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCAt            <  2:355792/243‑1 (MQ=255)
                                                                                                                                                                                                                                     tttctttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgcc         <  1:155909/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCCTgccaacccaa  >  2:73950/1‑250 (MQ=255)
                                                                                                                                                                                                                                            gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagccaa  >  1:71427/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
AACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAA  >  NC_000913/4625289‑4625774

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