Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,582,257	T→G	21.8%	T265P (ACT→CCT)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,582,257	0	T	G	21.8%	139.0 / 4.7	68	T265P (ACT→CCT)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base T (15/35);  new base G (10/7);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGT  >  NC_000913/4582026‑4582493                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              gcGGCAGCCGAAGAGTGGGGGGCAGATGACACTTACCCATCAGGTCACGCCGAATTGCGGTGCCTTGGCCGCCTTCAAACGGCACGTTATCCGGCC‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGAc                                                                                                                                                                                                                                         <  1:958799/237‑1 (MQ=255)               agggTGTGCAGATGCCCCTTATCAACAAGGCCACGAACAATTGCGGTGCCTTTGCCGCCTTAAAACCGCGCGTTATCCGGCC‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCgg                                                                                                                                                                                                                          <  1:51753/237‑1 (MQ=255)                tggtggGCGGATGACACTTACCCATCAGGTCCCGACGATTGTCGGTGCCTTGGCCGCCTTCAAACAGCCCGTTACCCGGCC‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGc                                                                                                                                                                                                                         <  1:407131/239‑1 (MQ=255)                  gtgtGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGGGGCGACAATATTGGCCTTCGGCAg                                                                                                                                                                                                                       >  1:185765/1‑239 (MQ=255)                   tgtgCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGGTAGGGGCGGCAATATTGGCCTTCGGCggg                                                                                                                                                                                                                      >  1:142888/1‑239 (MQ=255)                           tGACTCTTACCCATCGGGTCACGACGATGGTCGGTCCCTTTGCCTCTTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCAc                                                                                                                                                                                                               <  1:935131/238‑1 (MQ=255)                             acacTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTCCGGCGCTGCCAAACGGCGGGGTAGGGGGGCCCAATTGGGCCCTCGGCAGGTTTTTACCg                                                                                                                                                                                                             >  1:848295/1‑238 (MQ=255)                               acTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGCCGGGTGAACAAAGGTGCGGGTAAGGTTGGTGCCTGCGGCGCTGCCAACCGGCGGGTAAGTGGGGCAAAAAAGCGCCTCGGGCAGGGTTTCACCgggg                                                                                                                                                                                                          >  2:447784/1‑236 (MQ=255)                               acTTACCCATCGGGTCCCGGCGAATGGCGGTGCCTTTGCCGCCTTCAAACAGCACGTTACCCGGCC‑CCCCCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCg                                                                                                                                                                                                          <  1:378542/239‑1 (MQ=255)                                    tcccctaGGGCAGGACGGATGTCGGTGCCTTTGCCGCATACAACTCGCCCTTTACCCGGCA‑CCACCCCCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑Acc                                                                                                                                                                                                    <  1:590645/233‑1 (MQ=255)                                    tCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAAGGTTGGTGCCTTCGGCGCGGCCCAACGGCGGGGTATGGGCGGCAATATGCCCCTTCGGCGGGGTTTTAcccggcctca                                                                                                                                                                                                      >  2:990396/1‑230 (MQ=255)                                      cactcggccaccacGATTGCGGGTCCCTTCCCCCCCGTCAAACAGCCGTTTACCCGGCC‑CCCCCCCCCCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCa                                                                                                                                                                                                  <  2:142888/231‑1 (MQ=255)                                         cAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGGAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGGTAGAGGTCGCAATATGCACACTCAGCCAGGTTTCACCGCCGCG‑ACCCaagg                                                                                                                                                                                               >  2:462117/1‑236 (MQ=255)                                         cAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGGGCGGGTAATGTTGGTGCCTCCGGCGCTGCCAAACGGCGGGGTAGGGGCGACAATATGCGCCTTCGGCGGGGGTTCACCGTCGCT‑ACCCagcg                                                                                                                                                                                               >  1:625093/1‑237 (MQ=255)                                           ggTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGGTAGTGGCGCCAATATGCGCCTTTGGCCGGTTTTTACCGTGGCT‑ACCCCAAgtg                                                                                                                                                                                             >  1:662510/1‑239 (MQ=255)                                               acgacgAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCCGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACCATATGCGCCTTCGGCAGGGTTTTACCGCGGCTCCCCCAGAGTGtt                                                                                                                                                                                           >  1:566133/1‑238 (MQ=255)                                                 gacgaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAAGGGCGACAATATGCGCCTTCGGCCGGGTTTCACCGttctcttatacacatttccg                                                                                                                                                                                         >  1:681702/1‑219 (MQ=255)                                                   caacttcgttgcccTTGGCCGCCCTAAAACACCAC‑TTTTCGCGCACCCCCCACCGCCGCACGGC‑CGCCGGGATGCATCGTTTCGAAAATATGCTGCCTAAAGCACCACTTTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑CCCCAGAGTGTTGCCCAg                                                                                                                                                                                     <  2:229743/227‑1 (MQ=255)                                                      aTGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCCTAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTTGGGGGCACAATATGGGCCTTCGGCAGGTTTTCACCGTCGCC‑Accccccgcgtttccccgacg                                                                                                                                                                                  >  1:903842/1‑222 (MQ=255)                                                      aTGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGGGGCGTCAATATGCGCCTTTGGCAGGTTTTCACCGGTGCC‑TCCCCGAGGGTTGCCCAGGCg                                                                                                                                                                                  >  1:554946/1‑239 (MQ=255)                                                              gCCTTTGCCCCCTTCAAACAGCACGTTGCCCGGCA‑CCCCCCCCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTgcg                                                                                                                                                                           <  1:634091/238‑1 (MQ=255)                                                                              accAGCACGTTACCCGGCG‑CCACACCCGCCGCACGAC‑CGGCGGGATGCAGCGTTTCTATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGtt                                                                                                                                                          <  1:524693/237‑1 (MQ=255)                                                                               accacaccttttgcgccccccaccaccCCCGCAGGCC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTTGTGGTCGGATGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGtt                                                                                                                                                          <  2:848295/220‑1 (MQ=255)                                                                                   cACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGCCAATATGTGCCTTCGGCAGGGTTTCACCGTGCCC‑ACCCCGAGGGTTGGCCACACGGATTGGGCCCCCGGGGGCGGGGTTGCCt                                                                                                                                                      >  1:626682/1‑238 (MQ=255)                                                                                   cACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGCCAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTTGTCGGGGGTGCCt                                                                                                                                                      >  1:830623/1‑238 (MQ=255)                                                                                          tccccgccccaccccgcgcgcACGACGCGCTGGGATGCAGCGGTTGGATAATATGCTGCATAAAGCCAAACGGTTTGTTGCTGGTCGGGTGAACAAGGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAata                                                                                                                                                <  2:62163/221‑1 (MQ=255)                                                                                                     cgaccGACGCACGTC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcagg                                                                                                                                    <  1:429111/237‑1 (MQ=255)                                                                                                     cccccgccgCACGAC‑CGCCGGGATGCAGCGTTTCGATAATTTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcagg                                                                                                                                    <  1:709105/236‑1 (MQ=255)                                                                                                     ccccccgcccaaggc‑cgccgGGTTGGGGCGTTGTGGTAATATGCTGGATTACGCACATCTTGTTGTTGCGGGTCGGGTGAACAAGGGTGCGGGTAATGTTGGTGCCGGCGGCGCGGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcagg                                                                                                                                    <  2:554946/225‑1 (MQ=255)                                                                                                     ccaccGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGCCAAAATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGGGTTGCCCAGACGGATTGCGCCGCCGGGGTCGAGGTTGCCTTCAATATCTTgcagccgc                                                                                                                                    >  1:881494/1‑236 (MQ=255)                                                                                                         cgccgcACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGGCAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAAGGGGTTGCCCACACGGGTTTGGCCCCCGTGGTCGGGGGTGCCTTCAATATCGTGCAgcagggagg                                                                                                                                >  1:305722/1‑238 (MQ=255)                                                                                                            cgcACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCAc                                                                                                                                                                                                               <  1:747630/158‑1 (MQ=255)                                                                                                                   c‑cgccgGGATGCAGCGTTTCGATAATATGGTGCATAAAGCATAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGc                                                                                                                      <  1:951630/239‑1 (MQ=255)                                                                                                                            tggaaGCGTGTGGATAATTTGCGGACAAAAGTACAACTGTTGGTTGGGGGTGGGGTGAAAAAAGGGGGGGGGAATGGTGGGGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATAGGCGCCTTCGGAAGGTTTTCGCCGTCGCG‑ACCCAGAGGGTTGCCAAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAgacga                                                                                                               <  2:846296/234‑1 (MQ=255)                                                                                                                                  cGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCT‑ACCCAGAGGGTTGCCCCGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCCGTTCATCATGGCCAgacgcccgggg                                                                                                         >  1:689446/1‑231 (MQ=255)                                                                                                                                     ttCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGGGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCACTATTCTGCCGCAGGCAGTTCGACAGTGCCAGACGACGGGTGCCgg                                                                                                     >  1:869791/1‑239 (MQ=255)                                                                                                                                          taataTGGGGCATAAAGTCCAACCGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCAc                                                                                                 <  1:935272/238‑1 (MQ=255)                                                                                                                                          taataTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGa                                                                                                                                                                                <  1:31216/159‑1 (MQ=255)                                                                                                                                            atatGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATTCGCCTTCGGCAGGGTTTCAACGTCGCT‑CCCCAGAGCGTTGCCCCGACGGATTGGGCGCCCGTGGTCGGGGGTGCCTTCAATATCGTGCGGCAGGGAGGTCACTAGGGCCacggggaggggggcggggcccag                                                                                              >  2:1001020/1‑217 (MQ=255)                                                                                                                                                  ttcATAAAGCCCAACCGGTTGTTGGTGGTCGGGTGAACAAAGGGGGCGGTAATGTTGGTGCCTGCGGGGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCg                                                                                          <  1:1074234/235‑1 (MQ=255)                                                                                                                                                   gcATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGGGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGAGGACGGGGGCCGGGCACCaagtcgga                                                                                        >  1:573105/1‑231 (MQ=255)                                                                                                                                                    cATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGc                                                                                                                                                                                                         <  1:507268/125‑1 (MQ=255)                                                                                                                                                    cATAAAGCACAACTATTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGGGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCCGTTCATCAGTGCCCGACGACGGGTGCCcgggcccacgtcgatg                                                                                       >  1:98509/1‑225 (MQ=255)                                                                                                                                                       taaaccaactttttttttGCTGGGGGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGGGGGTTAGTGGCGAAAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGa                                                                                   <  2:573105/227‑1 (MQ=255)                                                                                                                                                       aaaGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGGTAGGGGCGACAATATGCGCCTTCGGCAGGGTTTTCCCGTCGGT‑ACCCAGAGGGTTGCCCCCGCGGAATGCGCCGCCGTGGGCGGGGTGTCCTTCAATAACCTGCACAAGGCAGCTCATCAGTGCCAGACGACGGGGGCCggggccccctttggggccgc                                                                                   >  2:181679/1‑222 (MQ=255)                                                                                                                                                              aaCTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCCGTTCGAGGCCGGTAAAcg                                                                             >  1:868416/1‑238 (MQ=255)                                                                                                                                                              aaCTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGGTAGTGGCGACCATATGCGCCTTCGGCCGGGTTTCACCGTCGCT‑ACCCCGAgt                                                                                                                                                                                              >  2:507268/1‑125 (MQ=255)                                                                                                                                                                  ggttgttgGTGGGCCGGGGAAAAAGGGTGCGGGTAATTGTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTGCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATGGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCgg                                                                        <  2:566133/237‑1 (MQ=255)                                                                                                                                                                     tgttgCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGGGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGGTGCCTTCAATATCGTGCAGCAGGCCGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCCAGGCCGATAAACGCCCGGTgg                                                                     >  1:724322/1‑239 (MQ=255)                                                                                                                                                                        tgCGGGTCGGGTGAACAAAGGTGGGGGTAATGTTGGTGCCTGCGGCGCTGCAAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATc                                                                  <  1:1068646/239‑1 (MQ=255)                                                                                                                                                                           tGGTCGGGTGAACAGAGGTGCGGGTAATGTTGGTGCCGGCGGCGTAGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑TCCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTgg                                                               <  1:1028519/239‑1 (MQ=255)                                                                                                                                                                                     aaCAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAgg                                                                                                                                                            <  1:642779/136‑1 (MQ=255)                                                                                                                                                                                         agggTGGGGGTAATGTTGGTCCCTGCGGCGGTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCGTGGGCAGGTTTTCCCCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCGGTGGTCGAGGTTGCCTTCCATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTc                                                 <  2:182154/237‑1 (MQ=255)                                                                                                                                                                                         aaGGTTCGGGTAATGTGGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGGTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTc                                                 <  1:198891/239‑1 (MQ=255)                                                                                                                                                                                                    aaTGGTGGTGCCTGCGGCGCTGCAAAACGGCGGGTGAGTGGCGACAATATGCGCCTTGGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATGGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAgg                                       <  1:40629/238‑1 (MQ=255)                                                                                                                                                                                                       gTTGGTGCCTGCGGCGCTGCTAACGGGCGGGTTAGGGGGGATAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGCGTGTTGCCCAGACGGAGTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGA‑CTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc                                  <  2:305722/239‑1 (MQ=255)                                                                                                                                                                                                        tttggtcttggggggccgCCAACCGGGGGGTTAGGGGCGACAATATGCGCCTTCGGCAGGTTTTCCCCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc                                  <  2:997495/222‑1 (MQ=255)                                                                                                                                                                                                        ttGGTGCCTGGGGCGCTGCCAAACGGGGGGTTAGGGGCTACAATATGGGGCTTCGGCAGGTTTTCACCGGCGAT‑ACCCGTGGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc                                  <  2:868416/239‑1 (MQ=255)                                                                                                                                                                                                        gtGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc                                  <  1:959800/238‑1 (MQ=255)                                                                                                                                                                                                           gTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGGACCAGTTCGAGGCCGATAAACGCGCGGTGGAACTGGAAATCCTGCCGGTCCCCATCAAGGTCGTCCAg                               >  1:997495/1‑239 (MQ=255)                                                                                                                                                                                                                      ccTGCCAAACGGCGGGTTAGGGGAGACAATTTGCCGCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGtt                     <  1:462117/237‑1 (MQ=255)                                                                                                                                                                                                                         gCCAAACGGCGGGTTAGGGGCGACAATATGCACCTTCGGGAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCg                 <  1:1056801/239‑1 (MQ=255)                                                                                                                                                                                                                          ccAACCGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCg                 <  1:181679/238‑1 (MQ=255)                                                                                                                                                                                                                            aagcggcggGTTAGTGGCGGAAATAGGCGCCTTCGGCAGGTTTTCACCGGCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACt              <  1:273099/236‑1 (MQ=255)                                                                                                                                                                                                                            aaaGGGCGGTTTAGGGGCGACAAATGGCGCCTTCGGGAGGTTTTCACCGTCGGT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGGCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACt              <  1:990396/239‑1 (MQ=255)                                                                                                                                                                                                                            aaaCGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACt              <  1:10397/239‑1 (MQ=255)                                                                                                                                                                                                                                     gTTAGGGGCGACAATATGCGCCTTCGGCAGGTTTGCCCCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc                                  <  1:495558/210‑1 (MQ=255)                                                                                                                                                                                                                                        aGTGGCGACAATTTGCCCCTTCGGCAGGTTTTCCCCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGt  <  1:405652/239‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              TCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCA‑CCACCACCGCCGCACGAC‑CGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGT  >  NC_000913/4582026‑4582493

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