Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,618,030 C→T G236G (GGC→GGT deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,618,0300CT98.2% 174.6 / ‑6.8 58G236G (GGC→GGTdeoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/31);  total (26/32)
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.61e-01

CTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCAT  >  NC_000913/4617786‑4618273
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
cTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTcc                                                                                                                                                                                                                                               >  1:248188/1‑251 (MQ=255)
cTGAAATCTCCATCAAAGCGGGAGCGGACTTCATCAAAACCTCTACCGTTAAAGTGGCTGTGAACGCGACGCCGGAAAGCACGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGGAATATCTCACCATCGCAGGTGAACTGTTCGGTGCTGACTGGGCAGAAGCGCGTCACTACCGATTTGGTGCTTcc                                                                                                                                                                                                                                               >  1:412618/1‑251 (MQ=255)
                         ggACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGCTGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGaa                                                                                                                                                                                                                        >  1:265940/1‑249 (MQ=255)
                          gACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACGACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAgc                                                                                                                                                                                                                     <  2:406432/251‑1 (MQ=255)
                                   aaaaCCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCCCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTgggcca                                                                                                                                                                                                            <  1:325566/251‑4 (MQ=255)
                                   aaaaCCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCACGCGCATCCTGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGCTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGAAGCGCAGAAATATCTCGCCCTCGCCGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCCCTAACGCTTTGGTGCTTCCAGCCTGCTGGCGAGCCTGATAAAAGCGCTgggccc                                                                                                                                                                                                            >  2:269478/1‑248 (MQ=255)
                                                   aaGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCAGCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAgc                                                                                                                                                                                            <  1:8151/251‑1 (MQ=255)
                                                       ggCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCca                                                                                                                                                                                        <  2:248188/251‑1 (MQ=255)
                                                          ggtgAACGCGACGACGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACCACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCcagca                                                                                                                                                                                     <  1:275757/250‑1 (MQ=255)
                                                              aaCGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGcta                                                                                                                                                                                 <  2:175420/251‑1 (MQ=255)
                                                                   gACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCCGctacta                                                                                                                                                                              >  2:303401/1‑249 (MQ=255)
                                                                      gCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACGGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACtaagtaa                                                                                                                                                                         <  2:265940/251‑1 (MQ=255)
                                                                         ggAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGGTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGAt                                                                                                                                                                      <  1:132117/251‑1 (MQ=255)
                                                                            aaGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCt                                                                                                                                                                   <  2:169116/251‑1 (MQ=255)
                                                                             aGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCtt                                                                                                                                                                  >  1:269110/1‑251 (MQ=255)
                                                                                   cgcATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGcc                                                                                                                                                            <  1:171057/251‑1 (MQ=255)
                                                                                   cgcATCATGATGGAAGGGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGGGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGcc                                                                                                                                                            <  2:58728/251‑1 (MQ=255)
                                                                                        catGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTGCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCTCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAt                                                                                                                                                        <  1:384185/250‑1 (MQ=255)
                                                                                         atgatgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgc                                                                                                                                                      >  1:248399/1‑251 (MQ=255)
                                                                                            attgAAGTGATCCGTGATATGGGCGTAGAAAAAACCGATGGTTTCAAACCAGCGGGCGGCGGGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCGGATGAGCTGTGCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTTGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgct                                                                                                                                                   <  2:12560/248‑1 (MQ=255)
                                                                                                     aTCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCCGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATc                                                                                                                                          >  2:327037/1‑251 (MQ=255)
                                                                                                        cGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAgg                                                                                                                                       >  1:323521/1‑251 (MQ=255)
                                                                                                                gggCGTAGAAAAAACCGTTGGGTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACgag                                                                                                                               <  2:229104/251‑1 (MQ=255)
                                                                                                                    gTAGAAAAAACCGTTGGGTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGTTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGt                                                                                                                           <  1:233540/251‑1 (MQ=255)
                                                                                                                          aaaaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGGAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAc                                                                                                                     >  2:266717/1‑251 (MQ=255)
                                                                                                                                   ggTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGcc                                                                                                            >  1:144794/1‑251 (MQ=255)
                                                                                                                                                cgggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGt                                                                                                 >  1:484448/1‑249 (MQ=255)
                                                                                                                                                cgggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCGCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAg                                                                                               <  1:416567/251‑1 (MQ=255)
                                                                                                                                                  ggcggcGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGAc                                                                                             >  1:285330/1‑251 (MQ=255)
                                                                                                                                                        gtgcgtACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATcc                                                                                       <  1:327037/251‑1 (MQ=255)
                                                                                                                                                         tgcgtACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCTGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCg                                                                                      >  1:59473/1‑251 (MQ=255)
                                                                                                                                                             tACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGGCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCACGCCTGTTGAAGGCGCGGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGCTGCGCTGCGATTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCa                                                                                   >  1:400856/1‑250 (MQ=255)
                                                                                                                                                                     aagaagCGCAGAAATATCCCGCCATAGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgccgcc                                                                          <  2:185888/251‑1 (MQ=255)
                                                                                                                                                                       gATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAgcgc                                                                                                                                                                                                                   <  1:195130/112‑1 (MQ=255)
                                                                                                                                                                       gATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAgcgc                                                                                                                                                                                                                   >  2:195130/1‑112 (MQ=255)
                                                                                                                                                                           cgcAGAAATATCTCGCCACCGCAGATTAACTGTTCGTGGCTTACTGGGCAGATGCGCGTCACTACCGTTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGTGCCTCGGTGACTTTAGGAGCGCCAGCAGCTACTAAGTAATATGCTTTACGCCTGAGGCGCTGCTCTTATCAGTCCTACGAGACGTATCTACACGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCAttc                                                                    <  1:311866/251‑3 (MQ=255)
                                                                                                                                                                             cAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATActc                                                                  <  1:88785/251‑1 (MQ=255)
                                                                                                                                                                               gAAATCTCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAACGCGCTGGGACACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCACAAGCCGCCTAATACTCtt                                                                >  2:324351/1‑251 (MQ=255)
                                                                                                                                                                                 aaTATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCtttt                                                              <  1:43075/251‑1 (MQ=255)
                                                                                                                                                                                       tcGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTCCCCGTAGCCCGGATAAGGCGTAGTCGCATCCGGCAACAGCCGCCTCATACTCTTTTCCTCgg                                                        >  1:167582/1‑251 (MQ=255)
                                                                                                                                                                                             ccgCAGATGAACTGTTCGGTGCTGACTGGGCATATGCGCGTCACTACCGCTTTGGTGCGTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACTCCTGATGCGCTTCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGAGAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGtt                                                  <  1:450691/249‑1 (MQ=255)
                                                                                                                                                                                               gCAGATGAACGGTGCGGTGTTGACTGGGCAGATGCGCGTCATTACCGCTTTGGTGCTTCGAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTAGCAGGCCTACGAGACGTATCGACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTAc                                                <  2:44433/251‑1 (MQ=255)
                                                                                                                                                                                                 aGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCt                                              <  1:397011/251‑1 (MQ=255)
                                                                                                                                                                                                      aaCTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGttt                                         >  1:122879/1‑251 (MQ=255)
                                                                                                                                                                                                       cctgttcgtgcatggcTGGGCAGATGCTCGTCTCTTCCGCTTTGGGTCTTCCAGCCGGCTGGTAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGGCCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCGGCGCTTATCAGGCCTACGAGACGTATCTTCCCGGAGGCCGGATAAGGCGTAGACGCATCTGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTtc                                        <  2:456359/236‑1 (MQ=255)
                                                                                                                                                                                                            tCGGTGCTGACTGGGCAGATGCGCGTCACTACCGATTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGCCGCATCCGGCAACAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGcc                                   >  2:371758/1‑250 (MQ=255)
                                                                                                                                                                                                             cGGTGCTGACTGGGCAGATGCGCGTCACTGCCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCGGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAGCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGca                                   <  2:138646/250‑1 (MQ=255)
                                                                                                                                                                                                                 gCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGcaca                                 >  1:274853/1‑248 (MQ=255)
                                                                                                                                                                                                                   tGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGaa                             >  2:147102/1‑250 (MQ=255)
                                                                                                                                                                                                                     aCTGGGCAGATGCGCGTCACTAACGCTTTAATGATTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCCGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTCCCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGAGAGGTTACCTTGTTTCTCGCACACGAAatt                          >  2:486761/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggggaagatgcgctgaccTACCGCGTTAGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAATCGCTGGGCTTCGGTGACGGGAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAattat                        <  2:340986/234‑1 (MQ=255)
                                                                                                                                                                                                                        gggCAGATGCGCGTCACTTCCGCTTTTGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTTCCTTGTTTCTCGCACAAGAAattatt                       <  1:218496/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggCAGAGGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTc                      <  2:248399/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCg                     <  2:479716/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTaa                  >  1:45908/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cgTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCAACAGAAATTATTCGTaaaaaa              >  2:8671/1‑249 (MQ=255)
                                                                                                                                                                                                                                        cTACCGCATTGGCGACTCCAGCCTGCTGGCAATCCTCTTGAACGAGCTGGACCACGGTGACGGTAAGCGCGCCAGCAGCTACTCAGAAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGCAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTATTTTCATCGGGAGGTTACCTTGATTCTCGCACAAGAAATTACTCGTAAAAAACgggctt       >  2:189000/1‑246 (MQ=255)
                                                                                                                                                                                                                                             gCTTTGGTGCTTCCATCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCAt  <  1:196459/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCAT  >  NC_000913/4617786‑4618273

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: