Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,585,230 A→G L511P (CTG→CCG)  hsdR ← type I restriction enzyme EcoKI endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,585,2300AG98.2% 192.9 / NA 56L511P (CTG→CCG) hsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/29);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGC  >  NC_000913/4584935‑4585518
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                
gTGATCGAGAATGCGAAGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTACAGTTCGCCTTCGGCCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCATCAACGACGATACAGTCGTAACGGGCCACCGACATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGAATGTACGGTGGCAAAGTGAATTTTGGTGCTGTCCTACGGAAATTTATCCGTCAGCCCTTTAATGTCTAAAACGCTGTTGAAGGCGTCGCAGTTAATACGCGCATCCTCAAACGCGACCAgcgc                                                                                                                                                                                                                                                                                              >  2:282321/1‑300 (MQ=255)
                    aGGCATAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTAGCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAAcggc                                                                                                                                                                                                                                                                          <  2:101826/300‑1 (MQ=255)
                    aGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCATTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTACGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGACGTTAATACGCGTATCTTCAAACGCGCCCGGCACCTGTTCGCCAaga                                                                                                                                                                                                                                                                                 >  2:45282/1‑293 (MQ=255)
                              gTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAg                                                                                                                                                                                                                                                                <  1:231403/300‑1 (MQ=255)
                                  tCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAag                                                                                                                                                                                                                                                            <  2:313154/300‑1 (MQ=255)
                                      gctggctgCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGaa                                                                                                                                                                                                                                                        <  1:200025/300‑1 (MQ=255)
                                        tggctgCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAGCGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCGTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAAt                                                                                                                                                                                                                                                       <  2:271713/299‑1 (MQ=255)
                                         ggctgCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGc                                                                                                                                                                                                                                                     >  1:598963/1‑300 (MQ=255)
                                           ctgCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTGTCGAGAATATAGCCGCGATGTGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCACCGGTTCATCTGATTGCAGGGTACGTTTCAACAGCGACTGTACGGTGGCAAAGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCCACCGCGCCCGGCGCCTGTTCGCACAGCGCACGGCGGTCGACAAGAAAGAGaa                                                                                                                                                                                                                                                        >  2:60870/1‑295 (MQ=255)
                                                        gTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCCCCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCAGGGAATTTATCCGTCGGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCAAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACgat                                                                                                                                                                                                                                         >  2:274180/1‑295 (MQ=255)
                                                        gTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGACTTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTCATACGCGTATCTTCAAACGCGCCCGGTGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCt                                                                                                                                                                                                                                         >  2:274004/1‑297 (MQ=255)
                                                               ttCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGAt                                                                                                                                                                                                                               <  1:542853/300‑1 (MQ=255)
                                                               ttCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGGCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGAt                                                                                                                                                                                                                               <  1:187236/300‑1 (MQ=255)
                                                                tCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATc                                                                                                                                                                                                                              <  1:9901/300‑1 (MQ=255)
                                                                tCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGGCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATc                                                                                                                                                                                                                              <  1:9708/300‑1 (MQ=255)
                                                                tCGGTCTGCTCTTGATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCAACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTGTCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGAt                                                                                                                                                                                                                               <  1:542862/299‑1 (MQ=255)
                                                                   gTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCa                                                                                                                                                                                                                             >  1:385570/1‑298 (MQ=255)
                                                                    tCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGc                                                                                                                                                                                                                          <  1:340093/300‑1 (MQ=255)
                                                                       gCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGa                                                                                                                                                                                                                       >  1:504732/1‑300 (MQ=255)
                                                                       gCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGCGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGAAAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGa                                                                                                                                                                                                                       >  2:381989/1‑300 (MQ=255)
                                                                                               gATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGGGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACggg                                                                                                                                                                                               <  2:582583/300‑1 (MQ=255)
                                                                                                        cGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTtacc                                                                                                                                                                                        >  1:227420/1‑298 (MQ=255)
                                                                                                            aactacgaTACAGTCGTAACGTGCCACCTGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTTGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTtaccggtac                                                                                                                                                                                   <  2:353221/295‑1 (MQ=255)
                                                                                                                      cAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATNCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatc                                                                                                                                                                        >  1:434271/1‑300 (MQ=255)
                                                                                                                        gTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatcgc                                                                                                                                                                      >  1:161049/1‑300 (MQ=255)
                                                                                                                            taacttgCCACCGGCATCGGTTCACCTGATTACAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCcagc                                                                                                                                                                  <  2:29324/294‑1 (MQ=255)
                                                                                                                                  gCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGAtc                                                                                                                                                            >  2:8380/1‑300 (MQ=255)
                                                                                                                                            tCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGcc                                                                                                                                                  <  1:596780/300‑1 (MQ=255)
                                                                                                                                               gTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGAGGAACATCATGGCGATTGTCGTACGGGTTTTACCCGTACCGGCCGCCATCGCCAGCAGGATCTCtt                                                                                                                                                        >  2:73693/1‑291 (MQ=255)
                                                                                                                                                    tcggATTGCGGGGTGCGTTCCACCAGCGACTGTACGGTGACAACGGGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCTTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACg                                                                                                                                          <  2:49749/297‑1 (MQ=255)
                                                                                                                                                                  gCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAAc                                                                                                                            >  2:537743/1‑300 (MQ=255)
                                                                                                                                                                  gCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAAc                                                                                                                            <  1:399190/300‑1 (MQ=255)
                                                                                                                                                                                aCTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCAt                                                                                                              >  1:454620/1‑300 (MQ=255)
                                                                                                                                                                                 cTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttg                                                                                                                                                       >  1:550536/1‑258 (MQ=255)
                                                                                                                                                                                 cTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttg                                                                                                                                                       <  2:550536/258‑1 (MQ=255)
                                                                                                                                                                                 cTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATc                                                                                                             >  2:599060/1‑300 (MQ=255)
                                                                                                                                                                                  tGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATc                                                                                                             <  1:472151/299‑1 (MQ=255)
                                                                                                                                                                                   gTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTANCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCtt                                                                                                           >  1:121640/1‑300 (MQ=255)
                                                                                                                                                                                   gTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTANCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCATTGACGATTGCCTTTTCAACCGCGCGGACGGCATCtt                                                                                                           >  1:122063/1‑300 (MQ=255)
                                                                                                                                                                                   gTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGCTTGCCTTTTCAACCGCGCGGACGGCATCtt                                                                                                           >  1:55516/1‑300 (MQ=255)
                                                                                                                                                                                   gTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCtt                                                                                                           >  1:55738/1‑300 (MQ=255)
                                                                                                                                                                                          ggCAACGTGAATTTTGGTGCTGTCTTCCTGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCGTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATTGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGata                                                                                                    <  1:423103/300‑1 (MQ=255)
                                                                                                                                                                                                 tGAATTTTTGTGCTTTCTTCCTGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTTTTGAAGTTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGTCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGTACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGTTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGc                                                                                             <  2:28672/300‑1 (MQ=255)
                                                                                                                                                                                                 tGAATTTTGGTGCTTTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGc                                                                                             <  1:221334/300‑1 (MQ=255)
                                                                                                                                                                                                           tGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTGCAAACGCGCCCGGCGCCTTTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTTTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgct                                                                                   <  2:35927/300‑1 (MQ=255)
                                                                                                                                                                                                               gTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgctcgct                                                                               >  1:199528/1‑300 (MQ=255)
                                                                                                                                                                                                                          aTTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGt                                                                    <  1:341798/300‑1 (MQ=255)
                                                                                                                                                                                                                                         cTTTAATGTCGAAAATGCTGTTGAAGGTGTCGACGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACt                                                     >  2:381801/1‑300 (MQ=255)
                                                                                                                                                                                                                                          tttAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCAc                                                      <  1:30974/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTc                   >  1:421066/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                               aTACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCATGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCcagc               >  1:338149/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         cTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGNGGAACATCATGGCGATTGCCGTACGGGTTTTANCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgct                                                                                   <  2:462133/222‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         cTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgct                                                                                   >  1:462133/1‑222 (MQ=255)
                                                                                                                                                                                                                                                                                             aaaCGCGCCCGGCTCCTGTTCGCCAAGAGAACGGTGGTCGACAAGGAAGAGAATGCGTGTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTCCGGGTTTTACCGGTACCGTTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTTACGATTGCCTTTTCAACCGCGCGGACGACATCTTCCTGATAATAACGCATGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAATCACTGATTCTGTTTTTGCGGTTCGCGGCCGAGCATTTCCAGCAGCTCTTCCGGGc  <  2:172949/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             aaaCGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc  <  2:498384/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             aaaCGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc  <  1:90517/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                
GTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGC  >  NC_000913/4584935‑4585518

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