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,2300AG96.7% 212.2 / NA 60L511P (CTG→CCG) hsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (26/32);  minor base T (0/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCC  >  NC_000913/4584950‑4585514
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                            
aCGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGACGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCgagaga                                                                                                                                                                                                                                                                           >  2:544739/1‑300 (MQ=255)
aCGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATACGTCAGCCCTTTAATGTCGAAAATGCCGTTGAAGGTGTCGCAGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAagcga                                                                                                                                                                                                                                                                           >  2:120052/1‑297 (MQ=255)
        cAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGt                                                                                                                                                                                                                                                                   >  2:34706/1‑300 (MQ=255)
                  aTCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGaa                                                                                                                                                                                                                                                         >  1:9206/1‑300 (MQ=255)
                   tCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAag                                                                                                                                                                                                                                                        <  1:586561/300‑1 (MQ=255)
                                      gCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCt                                                                                                                                                                                                                                     >  2:521121/1‑300 (MQ=255)
                                       cAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCt                                                                                                                                                                                                                                     >  1:16374/1‑299 (MQ=255)
                                           tCGCCTTCGGTCTGCTCTTTATCTAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGAc                                                                                                                                                                                                                                <  2:528485/300‑1 (MQ=255)
                                             gCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACt                                                                                                                                                                                                                               <  2:609356/299‑1 (MQ=255)
                                                     tCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCTTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGc                                                                                                                                                                                                                      <  1:15359/300‑1 (MQ=255)
                                                     tCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGc                                                                                                                                                                                                                      <  1:16609/300‑1 (MQ=255)
                                                               aTCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAAcatcat                                                                                                                                                                                                            <  1:614761/300‑1 (MQ=255)
                                                                  gagaATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGc                                                                                                                                                                                                         <  1:340926/300‑1 (MQ=255)
                                                                   agaATATAACCGCGCTGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCg                                                                                                                                                                                                        <  2:115920/300‑1 (MQ=255)
                                                                                gATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCANCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACggg                                                                                                                                                                                           <  1:464662/300‑1 (MQ=255)
                                                                                       ttCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTtacc                                                                                                                                                                                    <  1:494356/300‑1 (MQ=255)
                                                                                       ttCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTtacc                                                                                                                                                                                    <  2:635200/300‑1 (MQ=255)
                                                                                                 acgaTACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtc                                                                                                                                                                          >  1:219057/1‑300 (MQ=255)
                                                                                                         gTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatc                                                                                                                                                                    >  2:398758/1‑298 (MQ=255)
                                                                                                         gTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTACGGGAATTTATCCGTCAGCACTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatcga                                                                                                                                                                  >  2:440075/1‑299 (MQ=255)
                                                                                                               aCGGGCCACCGGCATCGGTGCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCcagcag                                                                                                                                                            <  2:521005/300‑1 (MQ=255)
                                                                                                                     cACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGAtctc                                                                                                                                                      >  1:305829/1‑300 (MQ=255)
                                                                                                                      aCCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGAtctc                                                                                                                                                      <  2:391391/299‑1 (MQ=255)
                                                                                                                            aTCGGTTCATCGGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTTATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttgttg                                                                                                                                                <  2:635279/299‑1 (MQ=255)
                                                                                                                                     tCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGAc                                                                                                                                       >  2:67434/1‑299 (MQ=255)
                                                                                                                                     tCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGAc                                                                                                                                       >  1:320894/1‑299 (MQ=255)
                                                                                                                                     tCTGATTGCAGGGTGCGTTTCAACAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCCGGATCTCTTGTTGCCCCTTGAc                                                                                                                                       >  2:92324/1‑299 (MQ=255)
                                                                                                                                               gggTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCtttt                                                                                                                            >  1:10539/1‑300 (MQ=255)
                                                                                                                                                   gCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAAc                                                                                                                        <  1:24661/300‑1 (MQ=255)
                                                                                                                                                   gCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAAc                                                                                                                        >  2:436295/1‑300 (MQ=255)
                                                                                                                                                   gCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTGATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATTGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAAc                                                                                                                        <  2:16374/300‑1 (MQ=255)
                                                                                                                                                               cGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGc                                                                                                            >  1:266142/1‑300 (MQ=255)
                                                                                                                                                               cGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTACGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCCACCGCGCGGACGGc                                                                                                            >  2:108732/1‑300 (MQ=255)
                                                                                                                                                                  cTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATc                                                                                                         >  1:604583/1‑300 (MQ=255)
                                                                                                                                                                         gTGGCAACGTGAATTTTGGTACTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAAAATCATGGCGATTGCAGTACGGGTTTTACCGGTACCGGTCGCCATCTCCAGCAGGATCTCTTGTTGCCCCTTGCCGATTGCCTCTTCAAACGCGCGGACGGCATCTTCCt                                                                                                    >  2:348529/1‑298 (MQ=255)
                                                                                                                                                                                        ttGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAGGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGccc                                                                                   <  1:37038/300‑1 (MQ=255)
                                                                                                                                                                                           gTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCNGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTACTGATAATAACGCAGGCCCAgc                                                                                >  2:190190/1‑300 (MQ=255)
                                                                                                                                                                                            tGCTTTCTTCCGGGAATTTATCCGTCAGCCCTTTAATTTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgct                                                                               <  2:121506/300‑1 (MQ=255)
                                                                                                                                                                                            tGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgct                                                                               <  1:583770/300‑1 (MQ=255)
                                                                                                                                                                                                       gggAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCTTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGcc                                                                     <  1:445210/299‑1 (MQ=255)
                                                                                                                                                                                                       gggAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTAGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGcc                                                                     <  1:66892/299‑1 (MQ=255)
                                                                                                                                                                                                          aaTTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAggg                                                                 <  1:577131/300‑1 (MQ=255)
                                                                                                                                                                                                           aTTTATTCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGt                                                                <  2:305829/300‑1 (MQ=255)
                                                                                                                                                                                                           aTTTATCCGTCAGCCCTTTACTGTCGAAGATGCTGTTGGAGTTGTCGCCGTTAATACGCATATCTTCAAACACGCCCGGCGCCTGTTCGCCAAGAGAATGGCGGTCGACAAGGAAGAGAGTGCGTTTAAAACGCTTGGACTGGATCAGGCGGAACATCATTGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAGCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGt                                                                <  2:288489/300‑1 (MQ=255)
                                                                                                                                                                                                           aTTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGTGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGt                                                                <  2:288501/300‑1 (MQ=255)
                                                                                                                                                                                                           aTTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCCTCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGt                                                                <  2:566125/300‑1 (MQ=255)
                                                                                                                                                                                                           aTTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGt                                                                <  1:49721/300‑1 (MQ=255)
                                                                                                                                                                                                            tttATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTTAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGNTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGt                                                                <  2:474278/299‑1 (MQ=255)
                                                                                                                                                                                                                   gTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGc                                                         <  1:583109/299‑1 (MQ=255)
                                                                                                                                                                                                                       gCCCTTTAATGTCGAAAATGCTGTTTAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCTGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAAcc                                                    <  2:219057/300‑1 (MQ=255)
                                                                                                                                                                                                                       gCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAAc                                                     >  1:218842/1‑299 (MQ=255)
                                                                                                                                                                                                                       gCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCATTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAAcc                                                    >  1:219186/1‑300 (MQ=255)
                                                                                                                                                                                                                              aaTGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGAtt                                             <  2:309827/300‑1 (MQ=255)
                                                                                                                                                                                                                                 gTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTANCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCAATGATTCt                                           >  2:209206/1‑299 (MQ=255)
                                                                                                                                                                                                                                             gTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTc                              >  1:445974/1‑300 (MQ=255)
                                                                                                                                                                                                                                                          ccGTTAATACGCGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCAtt                 >  2:385307/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                aTACGTGTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCcagc           <  2:318021/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTATCTTCAAACGCGCCCTGCGCTTGTTCGCCAAGAGAACGGTGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGctc      <  1:495627/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTATCTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGctc      <  1:495414/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cTTCAAACGCGCCCGGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCACAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTcc  >  2:230365/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                            
ACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGATACAGTCGTAACGGGCCACCGGCATCGGTTCATCTGATTGCAGGGTGCGTTTCACCAGCGACTGTACGGTGGCAACGTGAATTTTGGTGCTGTCTTCCGGGAATTTATCCGTCAGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCC  >  NC_000913/4584950‑4585514

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