Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,587,365 C→T S139S (AGC→AGT mrr → Type IV methyl‑directed restriction enzyme EcoKMrr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,587,3650CT96.2% 148.8 / NA 52S139S (AGC→AGTmrrType IV methyl‑directed restriction enzyme EcoKMrr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/28);  minor base G (0/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

AGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGT  >  NC_000913/4587125‑4587601
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
aGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACGTGCAAAGAGCAGTCCGGACGATc                                                                                                                                                                                                                                    <  2:567680/251‑1 (MQ=255)
 ggCTGGGCGCATGACCGTTTAAAACGTGCAGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCGCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCg                                                                                                                                                                                                                                   <  1:206795/251‑1 (MQ=255)
     gggCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGt                                                                                                                                                                                                                                 <  1:14653/249‑1 (MQ=255)
                cGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCt                                                                                                                                                                                                                    >  1:70804/1‑251 (MQ=255)
                cGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCt                                                                                                                                                                                                                    >  1:534298/1‑251 (MQ=255)
                            gCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGGCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCg                                                                                                                                                                                                        <  2:390781/251‑1 (MQ=255)
                                    gTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTAGTGATGCgg                                                                                                                                                                                                >  2:14879/1‑251 (MQ=255)
                                        cAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAAACATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGc                                                                                                                                                                                            >  2:465918/1‑251 (MQ=255)
                                                gTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAg                                                                                                                                                                                     <  1:137711/250‑1 (MQ=255)
                                                 tCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTGTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGTCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAgg                                                                                                                                                                                    <  1:433288/250‑1 (MQ=255)
                                                    cggggCAAATGGTGCTTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCt                                                                                                                                                                                <  2:441891/248‑1 (MQ=255)
                                                         cAAATGGTGCCTGACTCCTGCGGGTTTTGCCTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCCTCTGGCCTTCGCTTTTGTGACTGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGat                                                                                                                                                                            >  2:252418/1‑249 (MQ=255)
                                                         cAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGaa                                                                                                                                                                            >  1:272064/1‑250 (MQ=255)
                                                          aaaTGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGaaa                                                                                                                                                                           >  1:335794/1‑250 (MQ=255)
                                                           aaTGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAAc                                                                                                                                                                         >  1:532446/1‑251 (MQ=255)
                                                            aTGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACt                                                                                                                                                                        >  1:101166/1‑251 (MQ=255)
                                                                 gCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGTCTGATGAGGTTCTGTAAAACTTATTg                                                                                                                                                                   <  2:91797/251‑1 (MQ=255)
                                                                        tCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGCGGTTCTGGAAAACTTATTGCAGGttt                                                                                                                                                            >  1:168456/1‑251 (MQ=255)
                                                                                    tGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCt                                                                                                                                                >  2:236590/1‑251 (MQ=255)
                                                                                     gACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCt                                                                                                                                                <  2:17714/250‑1 (MQ=255)
                                                                                      aCTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGATTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCgcgcgtt                                                                                                                                              <  2:316685/251‑4 (MQ=255)
                                                                                                gTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTg                                                                                                                                    <  2:434404/251‑1 (MQ=255)
                                                                                                gTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTg                                                                                                                                    <  1:14879/251‑1 (MQ=255)
                                                                                                      tCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTgg                                                                                                                              >  1:260106/1‑251 (MQ=255)
                                                                                                            gCCAATGACGGAGCATGATACGAACCATCTGGACTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCACCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGAGCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGTTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGttt                                                                                                                        <  1:24351/251‑1 (MQ=255)
                                                                                                             ccAATGACGGAGCAGGAGACGAACCCTCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCTAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGAGCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGAGGCGGTGGCTGATGAGGTTCTGGAAATCTTATTGCATGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGtttt                                                                                                                       <  2:521310/251‑1 (MQ=255)
                                                                                                             ccAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGtttt                                                                                                                       <  1:465918/251‑1 (MQ=255)
                                                                                                                  gACGGAGCAAGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGGCAATCTGAAGTCACGGGCGGATGCCTTCGATTTAGATCCGAACGCCGACTCTCCCTATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTGCGTGATGCGGTTGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGGCATTGTTCTGTATGTTTTGCATc                                                                                                                  <  2:211996/251‑1 (MQ=255)
                                                                                                                     ggAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGcc                                                                                                               <  2:70804/251‑1 (MQ=255)
                                                                                                                      gAGCAGGAGACGCACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCt                                                                                                              >  1:440328/1‑251 (MQ=255)
                                                                                                                          aGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTgggg                                                                                                          <  1:123820/251‑1 (MQ=255)
                                                                                                                             agaCGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATGGCAGGTTGCTCCTTCGCGCTTTGAAGTCATTGTTCTTGATGTTTGGCATCGCCTGGGGTAt                                                                                                       <  2:263091/251‑1 (MQ=255)
                                                                                                                               aCGAACCATCTGGCCTTCGCTCTTGTGAATGTCCAAATTACGCCACGGCCAGATGCCGTCGATTTAGATCCGCAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGCAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATgg                                                                                                     >  2:111783/1‑251 (MQ=255)
                                                                                                                                     cATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCa                                                                                               <  1:56844/251‑1 (MQ=255)
                                                                                                                                                cGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGAtt                                                                                    >  1:386092/1‑251 (MQ=255)
                                                                                                                                                               cAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGGAGCGTGTTGGc                                                                      <  1:152428/250‑1 (MQ=255)
                                                                                                                                                               cAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCg                                                                     >  1:144926/1‑251 (MQ=255)
                                                                                                                                                                     tacgTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGGACTTGCAAAGAGCAGTCAGGCCGATCGGTTAGATCAGGAGCTAAAAGAGCTTCGGGGTGCGGGGGCTGAAGAGGTTCTGGAAAACTTATTGCAGGCTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGGATGGAGGCCACCGTGGTGGTTTGCAACGAGTTGGCGGtgcgg                                                               >  2:311247/4‑247 (MQ=255)
                                                                                                                                                                     tAAGTCACGGCCGTATGCCGTCGATTTAGATCCGCAAGCCGACTCTCCCGCTCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGCATGGCGGCCACCGTGATGATTTGCAGCGGGTTGGCGGTTCTg                                                               >  2:418630/1‑251 (MQ=255)
                                                                                                                                                                          cACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTggcggt                                                          >  1:7015/1‑247 (MQ=255)
                                                                                                                                                                             ggCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGGtggg                                                       >  1:219277/1‑250 (MQ=255)
                                                                                                                                                                                            aTTTAGATCCGAAAGCCGACTCGCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTTCTGGAGATGGTGGTATCGATGGtgtg                                        <  1:496872/251‑1 (MQ=255)
                                                                                                                                                                                                    ccGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCCGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACGGGAGATGGTGGCATCGATGGTGTGATATCGc                                 <  1:200323/250‑1 (MQ=255)
                                                                                                                                                                                                     cgacaGCCGCCTCTCCCGATCAAGAAGACCTGGCAAAGAGCATGCCGGACGATCGTTTAGAACAGGCGCTAAAAGAGCGTCGTGATGCGGTGGTCGATTTGGTTCTGGAAAACTTATTGCTGGTTATTCCTTCGCGATTTGAAGTCATTGTGCTGGATGTCGTGCATCGCCTGTGGTATGGCGGCCACCGTGATGATTTGCAGAGTGTTGGCGGTACTGGAGATGTTGGCATCGATGGTGTGATATCGCtt                               <  1:578863/247‑1 (MQ=255)
                                                                                                                                                                                                     cGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCtt                               <  1:118871/251‑1 (MQ=255)
                                                                                                                                                                                                                 ctcCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCGTGGAAGTCATTGTTCTGGATGTTTTGCATCGCCGGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTg                     <  2:287885/249‑1 (MQ=255)
                                                                                                                                                                                                                   cccGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTgg                    <  1:297042/248‑1 (MQ=255)
                                                                                                                                                                                                                     cGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTgg               <  2:410297/251‑1 (MQ=255)
                                                                                                                                                                                                                           tGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGCTGATTTGCCGCGTGTTGGCGGGACTGGAGATGGGGGCATCGATGGGGTGGTATAGCTTGTTAAACTTGGCCTGGCGaa           >  2:579834/1‑249 (MQ=255)
                                                                                                                                                                                                                           tGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGga                                                             >  1:255810/1‑199 (MQ=255)
                                                                                                                                                                                                                           agaagaaCTTGCAAAGAGCAGTCCGGACGATGGGTTAGATCAGGCGCTAAAAGAGCTGCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTGCTCCTTCGCGCGTTGAAGTCATTGGTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTGGCAGCGTGTTGGCGGTACTGga                                                             <  2:255810/198‑1 (MQ=255)
                                                                                                                                                                                                                                    tGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTTGCTGATGAGGTTCTGGAAAACTTATTGCAGGTGTCTCCTTCGCGCTTTGAAGTCATTGTTCTGTATGTTTTGCAGCGCCTGGTGTAGGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTTATAAACGTGGCCTGGAGAAATTTTATGt  <  2:579641/249‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
AGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGT  >  NC_000913/4587125‑4587601

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