Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,587,589 A→G 14.0% E214G (GAG→GGG)  mrr → methylated adenine and cytosine restriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,587,5890AG14.0% 185.9 / 3.1 71E214G (GAG→GGG) mrrmethylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base A (23/37);  new base G (9/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  NC_000913/4587353‑4587826
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
tGCAAAGAGCAGCCCGGGCGATCGGTTAAATCAGGCGGTAAAAGAGCTTCGTGATGCAGTGGCTGATGAGGTTCTGTAAAACTTATTGCATGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGgaga                                                                                                                                                                                                                                                 <  1:462349/239‑1 (MQ=255)
tGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTTGCGGGACTTGAGATGGGGGGATCGCTGGTGTGGTATCCCTTTGTAAACTTGGCCTGggga                                                                                                                                                                                                                                                 >  1:158948/1‑239 (MQ=255)
            cccGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGGGATGGTTTTCAGCGTGTTGGGGGGACTGGGGATGGGGGGATCCCTGGGGTGATAACGCTGGTAAAAAAAGGCCTGGGGAAAAATTATGttt                                                                                                                                                                                                                                     >  2:319648/1‑238 (MQ=255)
             ctggACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGtt                                                                                                                                                                                                                                      <  1:560156/235‑1 (MQ=255)
                     tCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTTCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGGTAAAATTGGCCTGGAGAAAGGTTATGTTTAGGCAAAAc                                                                                                                                                                                                                            >  1:128861/1‑239 (MQ=255)
                      cGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGGTGGGGTGGTATCGCTTGGTAAACTTGGCCTGGGGAAAGTTTTTGTTGTGGGAAAACg                                                                                                                                                                                                                           >  1:393618/1‑239 (MQ=255)
                       ggTTTGAGCTGGCGCTAAAAGAGCTTCGTGACGCGGTGGCTGAGGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGGTTTGAAGGCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGt                                                                                                                                                                                                                          <  2:421022/239‑1 (MQ=255)
                       ggTTAGATCAGGCGCTAAAAGGGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAATCTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGt                                                                                                                                                                                                                          <  1:913755/239‑1 (MQ=255)
                        gttaagtCAGGCGCTAAATGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGGCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCGACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTTAGGCAAAACGt                                                                                                                                                                                                                          <  1:227767/232‑1 (MQ=255)
                            gATCAGGCGCTAAAAGTGCTTCGTGATGCGGTGGTTTAGTGGGTTCTGGAAATCTTATTGCAGGGTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCa                                                                                                                                                                                                                     <  1:435386/239‑1 (MQ=255)
                              gtaGGGGCTAAAAGAGGTTCGGGTTTCGGTGGCTGGTGAGGTTCTGGAAAACTTATGGCAGGTTTCTGCTTGGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGGGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                                                                                   <  2:667103/237‑1 (MQ=255)
                               cAGGCGCTAAAAGGTGTTCGTTATGCGGTGGCTGATGAGGTTCTGGAAAACTTAGTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                                                                                   <  1:797719/238‑1 (MQ=255)
                               cAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCACCGCCTGGGGAATGGCGGCCACCGTGTGGATTTGCAGGCTCTTGGCGGTACTGGGGACGGGGGCATGAAGGGGGGGATCTGGCGTGATAAACTTGCACTGGAGGAAGTGTTGGTTGGGGCAAAAAATTGGCGGaa                                                                                                                                                                                                                  >  2:762077/1‑239 (MQ=255)
                                       aaaaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTTGAGATGGTGGCATCGATGGGGTGATATCGCTTGTTAAACTTGGCCTGGGGAAAGTTTTTGTTCAGGGAAAACTGTGGCAGAATACTGTag                                                                                                                                                                                                          >  1:225937/1‑239 (MQ=255)
                                        aaaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCAGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGGATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTagg                                                                                                                                                                                                         <  1:103904/239‑1 (MQ=255)
                                                tCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGGGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTTGCCTGGAGAAAGTTTATGGTCAGGCAAAACGTTGGCGGAATACTGTAGGCAGGCCAg                                                                                                                                                                                                 >  1:182732/1‑239 (MQ=255)
                                                      tggggtgggtgGTGGGGTTTTGGAAAATTTTTTGCAGGTTTCTTCTTCGGGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTAc                                                                                                                                                                                          <  2:771922/236‑1 (MQ=255)
                                                           tGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGGGGCATCGATGGGGTGATATCGCTTGGTAAACTTGGCCTGGGGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCa                                                                                                                                                                                     >  1:82110/1‑239 (MQ=255)
                                                             gCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTTTTTGCGGTACTGGAGATGGTGGACTCGCAGGTGGGGTTTCCCTTGTTTAAATTGTGCTTGGGGAAGATTAGTTTTCGGCAAAAAGTGGGACGAAAAAGTTAGGGGGGCCAGAAATAAAAGCAtt                                                                                                                                                                                   >  2:685232/1‑239 (MQ=255)
                                                                  tggggTTCTGGAAAACTTATTGGAGGTTTCCCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACg                                                                                                                                                                              <  1:319648/236‑1 (MQ=255)
                                                                  tgaGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTAc                                                                                                                                                                               <  1:583342/238‑1 (MQ=255)
                                                                                 cTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTTCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAATTGGCCTGGGGAAAGTTTATGTTCAGGCAAAAAGTTGGCCGAATACTGTAGGCAGGCCCGAATTACAGGCATTTTACGGGGCCCCGGCTGGGc                                                                                                                                                              >  1:152139/1‑239 (MQ=255)
                                                                                           gTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGGTATCGCCTGATAAACTTTGCCTGGGGAAAGGTTATGTTCCGGCACAACGTTTGCAGAATACTGTAGGCAGGCCCAAATTACAGGCAATTTAAGGCGCCCTGGCTTGGGAAAAAACGaa                                                                                                                                                  >  1:892692/1‑239 (MQ=255)
                                                                                            ttttcccTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGaa                                                                                                                                                  <  1:519468/233‑1 (MQ=255)
                                                                                                cccTTCGGGTTTTGTAGTCTTTTTTTTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTGCAGCGTGGTGGCGGTTCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTg                                                                                                                                             <  2:393618/238‑1 (MQ=255)
                                                                                                     cgcgCTTTGAAGTCCTTGTTCTGGATGTTTTGCATCGGCTGGGGTAAGGAGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGgtg                                                                                                                                        <  1:717542/239‑1 (MQ=255)
                                                                                                             gAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTAc                                                                                                                                <  1:59638/239‑1 (MQ=255)
                                                                                                                           ggATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACCGCGCACCTGCTGGGGAAAAAGGCGAAAGTGGGGTGGTTAATACCAATTTCGGGTTTa                                                                                                                  >  1:714162/1‑239 (MQ=255)
                                                                                                                            gATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGGGGGATATCGCTTGGAAAACTTGGCCTGGGGAAAGTTTTTGTTTATGCAAAACGTTGGCAGAAAACTGTATGAAGGCCAGAATTTAAGGGAATTTTATGCGCACCTGTCGGGCAAAAAAGGGAAACAGGGGGGGTTATTTACCCCTCTTGGTTTTa                                                                                                                  >  2:777436/1‑239 (MQ=255)
                                                                                                                                 tttGCATCGCCTGGGGTATGGCGGCCACCGTGGTGATTTGCAGCGTGTTGGCGGTACTGGGGGTGGTGGCATCGCTGGGGGGGTATCTCTTGTGAAAATTTGCCTGGGGAAAGGTTTTGTTCCGGGAAAACTGTGGCCGCATACTGTAGGCCGGACAGCATTACCAGGCTTTTACGGCGGCCGACACGGGCAAAAAGCAGAAAAGGGGGGGGTTTATAACCACTTTTgttttattttta                                                                                                             >  2:470632/1‑228 (MQ=255)
                                                                                                                                             ggggTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGTCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACt                                                                                                <  1:794784/239‑1 (MQ=255)
                                                                                                                                                    ggcggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGGGTGGTATCGCTTGATAAAATTTGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAAAACTGTAGGCCGGCCAGACTTACAGGCATTTTTCGGCGCACTTGCCGGGGAAAAAGCGGAACATGGGGGGTTTTTTACCCCTTTTGGATTTACTTTTCAGGCGCGGGACCTTGCCCa                                                                                         >  2:783477/1‑239 (MQ=255)
                                                                                                                                                      cggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAGGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCaa                                                                                        <  1:648243/238‑1 (MQ=255)
                                                                                                                                                          cACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:470632/239‑1 (MQ=255)
                                                                                                                                                           aCCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATAGCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:184541/238‑1 (MQ=255)
                                                                                                                                                             cGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                                                                                <  1:186069/239‑1 (MQ=255)
                                                                                                                                                                   gttttGCCGCGTGGTGGGGGTTCTGGAGAGGGTGGCATCGGTTGGGTGATAACGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGTACTGGCTGGTCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTa                                                                           <  1:507668/236‑1 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGgtg                                                                      <  1:783477/239‑1 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGGGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATgggg                                                                      >  1:354448/1‑237 (MQ=255)
                                                                                                                                                                         cAGCGTGTTGGCGGTCCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGCATACTGTAGGCAGGCCAGATTTACAGGCATTTTACGGCGCCCTGGCTGGGCAAAAAACGAAACGTGGGGGGGTTATTACCACTCTTGGATTTAATTCTCAGGCGCGGGACTTTTCCCACACCGCCGCGGGTATGgttgt                                                                    >  2:836302/1‑236 (MQ=255)
                                                                                                                                                                             gttttGGGGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTgttggt                                                                 <  2:619747/235‑1 (MQ=255)
                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGGGTTTATTACCACTTCTGGGTTTACTTTTCAGGCGCGTGCCTTTGCCCAAACCGCCGAGGGTGTGGGgtgggttg                                                               >  2:552278/1‑239 (MQ=255)
                                                                                                                                                                               gtTGGCGGTACTGGAGAGTGTGGCATCGATGGTGTGATAGCGCTTGAGTAACTTGGCCGGGAGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGa                                                              <  1:739068/239‑1 (MQ=255)
                                                                                                                                                                                gttgCGGTACGGGTGTTGGGGGCATCGATGGTGTGATTTCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAGAACGTTGGCAGAATTCTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTTGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGAt                                                             <  2:182732/236‑1 (MQ=255)
                                                                                                                                                                                     ggTATTGGGGATGGGGGGATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  2:152139/238‑1 (MQ=255)
                                                                                                                                                                                     ggTACGGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTTGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  1:762077/238‑1 (MQ=255)
                                                                                                                                                                                      gttcTGGAGATGGGGGAATCGATTGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGa                                                        <  2:128861/235‑1 (MQ=255)
                                                                                                                                                                                        aagggAGATGGTCGCATCGATGGTTCGATATCGTTTGATAAACTTGGCCTGGGGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCGCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                      <  1:735418/235‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTggg                                                                                                                                           <  1:301937/152‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTggg                                                                                                                                           >  2:301937/1‑152 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATTGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 <  2:56524/206‑1 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 >  1:56524/1‑206 (MQ=255)
                                                                                                                                                                                                tggtggGATCGATGGTGTGATATCGGTTGAAAAACTTGGCCTGGAGAAAGTTTAAGTTCAGGCAAAACGGTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  2:212083/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:552278/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCCCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:304931/238‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  >  2:744421/1‑201 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  <  1:744421/201‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGGGACTTTGCCCAATCCGTCGCGGGAATGGTGTTGGTTGATGGGGAACGCCTGGGGCa                                            >  1:198335/1‑239 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGGGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGc                                             >  1:417280/1‑238 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGAGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGGGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAGGCCGGGTg                                              >  1:450720/1‑237 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGCGAGGGGTTTATTACCCCTTCTGGATTTCCTTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATGGGGTTGGTTGATGGGGGACGCCGGGGGc                                             >  2:703988/1‑238 (MQ=255)
                                                                                                                                                                                                               ttggATATCGCTTGATAAACTTGGCCGGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTaa                                       <  1:843142/227‑1 (MQ=255)
                                                                                                                                                                                                               tgtgATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGCATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGaaa                              >  1:862724/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGGAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCAGTTTACGGCCCCCTGGCTGGGCAAAAAGCGAACCGTGGGGTGTTTACTACCACTTCTGGATTTACTTCTCAGGAGCGTCACTTTGCCCCATCCGTCGGGGGTATGGCGTTTGTGACAGGGGAAACCCTGGTGCCCCTACAGAACCAAAACGAaggagg                    >  2:37168/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCCTGCCAGAATTACAGGCATTTTACGGCTCACTGGCTGGGCAAAAATCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTTCTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGGCGGTCTGTGTTTGGTTGGTGGGGCACGCCTGGTGCACTTAATGATAGAAAACGAAGTAgg                    >  2:776788/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAgg                    >  1:301741/1‑239 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTgggg                                                                                                                                          <  2:606338/117‑1 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTgggg                                                                                                                                          >  1:606338/1‑117 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGtt                >  1:687926/1‑239 (MQ=255)
                                                                                                                                                                                                                                    tGGCCTGGGGAAATTTTTTTTTCAGGTAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCAc         <  2:354448/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGGGAGAAAGGTTAAGTTCGGGCAGAACGTTGGCAGAATCCTGTAGGCAGGCCAGGATTACAGGCATTTTGCGGCGCACTGGCTGGGCAAAAAGCGAAACGGGGGGTGTTTATTACCTGTTCTGGATTTACTTCTCAGGTGCGTGACTTTGCGCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGt       <  1:37168/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGtt      >  1:140634/1‑239 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGTTTGGTTGATGGGGAACGCCTGGTGCCCTTAATGATCGAAAACGAAAAAGGGGTGTTTTCAAGtt      >  2:710044/1‑239 (MQ=255)
                                                                                                                                                                                                                                           gagaAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATTGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGCAAACGAACTAGGGGTTTCTTCCCGTttgtt  >  1:226147/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  NC_000913/4587353‑4587826

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