Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,366,007	A→C	14.1%	G26G (GGT→GGG)	rhmD ←	L‑rhamnonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,366,007	0	A	C	14.1%	232.2 / 2.1	73	G26G (GGT→GGG)	rhmD	L‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (56/3);  new base C (2/12);  total (58/15)
Fisher's exact test for biased strand distribution p-value = 2.05e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAG  >  NC_007779/2365773‑2366245                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               acaTTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCAccg                                                                                                                                                                                                                                             >  1:1336205/1‑238 (MQ=255)   aTTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCAccgcc                                                                                                                                                                                                                                           >  1:691648/1‑238 (MQ=255)   aTTTACACTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCAccgcc                                                                                                                                                                                                                                           >  1:1652760/1‑238 (MQ=255)     ttACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCCCGTAATCGCGGTATTTACTCATCGGGGGGGCAAT‑ATGATCGTCTATCCAG‑TGTTTCGCCCCCTGGTC‑GTGATAATCACCgccgcc                                                                                                                                                                                                                                        >  1:117032/1‑239 (MQ=255)             aaTGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCACAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCtt                                                                                                                                                                                                                                >  1:1596148/1‑239 (MQ=255)               tGAAACGGTTAAGATGTTTTTCGACAATAAAACAGTCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGGCTGCTCGGAATCCCGGGATTTACTCCTCGGGGGGGCAAT‑ATTATCGTCTATCCAG‑TGATTCGCCCCCTGGCT‑GGGGTAATCCACCCCGCCAGCCCCtttt                                                                                                                                                                                                                              >  2:321735/1‑239 (MQ=255)                   aCGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTct                                                                                                                                                                                                                           >  1:993965/1‑238 (MQ=255)                       ttAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTct                                                                                                                                                                                                                           >  1:1500797/1‑234 (MQ=255)                       ttAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTct                                                                                                                                                                                                                           >  1:546362/1‑234 (MQ=255)                       ttAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGGGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCCGGTC‑GGTAGAATCACCCCCCCCAGCGCCTTTTTttgctt                                                                                                                                                                                                                       >  2:415778/1‑237 (MQ=255)                       ttAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCCCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCCCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTct                                                                                                                                                                                                                           >  1:157823/1‑234 (MQ=255)                               tttttCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCg                                                                                                                                                                                                               >  1:1000884/1‑238 (MQ=255)                                 ttacGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAGTACGG‑TCTGGCCGCTTTTCGC‑CTCGGCTTCAACAACCAAAGTGCCAACAA‑CCTTTATGCCAAAT‑GGCTGGCTGGAATGCTCCCAATCGCGGTGTTTAAACAAACGGGGGGCAAA‑TTTATCCCTCATCCAG‑AGAATCGTCCCCCGGTT‑CCGACAATCCACGTCCCCCGCGCTCTTTTTTGCTGTTGCAcccgca                                                                                                                                                                                                            >  2:698881/4‑235 (MQ=255)                                                ctcccAATTCCTCGGCTAGCAAAACGGGCGATCCGGTTCTGGCC‑TTTTCTTCGTTGTATTTCACCACCCAACGTCCAAAGAA‑CTTTATTCCCAAAT‑GATTGGGTTGATTGCTCGTAATCGCGTTTTTTACTCAGGGGGTTGGAAAT‑ATGTTCGTTTATCCAG‑T‑ATTCCCCCCCTTGTCGGTGATAACCCCCGCCGCCGGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGcgcg                                                                                                                                                                                              <  2:1305133/237‑1 (MQ=255)                                                  cccATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                                                                                                                                                                            >  1:1735047/1‑238 (MQ=255)                                                  cccATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                                                                                                                                                                            >  1:1062056/1‑238 (MQ=255)                                                    cATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                                                                                                                                                                            >  1:1752002/1‑236 (MQ=255)                                                    cATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                                                                                                                                                                            >  1:1019961/1‑236 (MQ=255)                                                        cccccGCCTGTCGAAACGGCGATTCCGG‑TTGGCCCGTTTCCTGC‑TTCGCCTTCACCACCCAACGTCCCAAAAC‑CTTTAATCCAAAAG‑GACTGGCGTGCTTTCTCGTACTCCGGGTTTTTCCCCTCGGGGGTGGAAAT‑ATGTTCGTTTTTCCGG‑TGTTCCCCCCCCTGTTG‑GTAATAACCCCCGCCGCCACCGTTTTTTTTTGTTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGttt                                                                                                                                                                                     <  2:1494005/236‑1 (MQ=255)                                                                 gTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg                                                                                                                                                                            >  1:15413/1‑239 (MQ=255)                                                                       aCGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAg                                                                                                                                                                                                                                      >  2:376884/1‑175 (MQ=255)                                                                       aCGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAg                                                                                                                                                                                                                                      <  1:376884/175‑1 (MQ=255)                                                                             attCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTCCCAAAAA‑CGTAAATGCAAAAT‑GACTGGCGTGACTGCTCGAAATCGGGGTATTTACTCATCGGGGTGGAAAT‑ATGATCGTCTATCCAG‑TGATTCCCCCCCTGGTC‑GTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat                                                                                                                                                                <  1:546618/237‑1 (MQ=255)                                                                              ctcggg‑TCTGGCGTTTTTTTGC‑TACGACTACAACACCCAACGTGCCAAGAA‑TTTAAATGCAAAAG‑GCGGGGGGTGACTGCTCTTAATGGCGGTATTTCCTCACGGGGGTGGAAAT‑ATGTTCGTCTATCCAT‑TGATCCCCCCCCGGGTC‑GTAAAAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat                                                                                                                                                                <  2:1752002/234‑1 (MQ=255)                                                                                 cGG‑TCTGGCCCTTTTCTTC‑CTCACCCTCAAACACCAACGCGCCAAAAACCGTTAATACAAAAT‑GCCGGGTGTGCTTCCTCGTACGCGGGTTTTTTACTCAGGGGGGTAGAAAT‑TTTATCTTCTACCCTG‑TCTTCCCCCCCCGGGCC‑GGAATAACCCCCGCCCCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGtt                                                                                                                                                             <  2:1769824/239‑1 (MQ=255)                                                                                    gttgggCCTTTTTCTGT‑GTCCACTTAACCAACAACCGTGCAAAGAA‑CTTTAAT‑CCAAATGGACGGGGTTGATTCCTCGAAACCGGGGT‑TTTTCTCCATGGGGGGGCAATAATGTTTTTCTTTCCAG‑TGTTCCCCCCCTTGGTC‑GTGATAACCCCCGCCGCAAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTtct                                                                                                                                                           <  2:1154937/234‑1 (MQ=255)                                                                                     tCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTcc                                                                                                                                                         >  1:530453/1‑239 (MQ=255)                                                                                      cTGGCCTTTTTCGGC‑TTCGACTTCAACAACCAACGTGCAAAGAC‑CGTTAATGCCAAAT‑GACTGGCGTGCCTGCTCGTAATCGCGTTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCa                                                                                                                                                        <  1:1216190/239‑1 (MQ=255)                                                                                         gCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat                                                                                                                                                                >  1:616455/1‑228 (MQ=255)                                                                                         gCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatg                                                                                                                                                                  >  1:528158/1‑226 (MQ=255)                                                                                              tttCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATg                                                                                                                                                >  1:568511/1‑239 (MQ=255)                                                                                                               aacaacCAACTCGCCAAGAA‑CTTAAATGCCAAAG‑GACGGGCGTGCCTGTTAGTAACCGCTGTTTTTCCTCATCGGGGTGGAAAT‑ATGATCTTCTACCCGG‑TGTTCCGCCCCCTTGTC‑GTGATAAACCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCaaa                                                                                                                                <  2:1000884/239‑1 (MQ=255)                                                                                                                    ccaacGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCCGGCCCGAACCTGCTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAggg                                                                                                                             >  2:1637749/1‑237 (MQ=255)                                                                                                                        cGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑GTGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATATTCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACttt                                                                                                                        >  1:1312180/1‑238 (MQ=255)                                                                                                                        cGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTACTTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGCCttt                                                                                                                        >  1:555828/1‑238 (MQ=255)                                                                                                                         gTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAg                                                                                                                                              >  1:918369/1‑215 (MQ=255)                                                                                                                            ccAAGAA‑CTTAAATCCCAAAT‑GACTGGCGGGACTGCTCGTAACCGCGGTTTTTACTCGCGGGGGTGGAAAT‑ATGATCGTCTACCCAG‑TGATCCCCCCCCTGGTC‑GTAAAAATCCCCGCCGCCAGCCCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCAt                                                                                                                   <  1:698881/239‑1 (MQ=255)                                                                                                                              aaGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATc                                                                                                                  >  2:926718/1‑238 (MQ=255)                                                                                                                                               aaT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAata                                                                                                 >  1:1109334/1‑239 (MQ=255)                                                                                                                                                         cGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTaaa                                                                                                                                                                                                       >  1:1550025/1‑128 (MQ=255)                                                                                                                                                                  tCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTaaa                                                                                                                                                                                                       >  1:697459/1‑119 (MQ=255)                                                                                                                                                                           cgGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGccc                                                                      >  1:1654973/1‑239 (MQ=255)                                                                                                                                                                             ttATTTCCTCACCGGGGTGGAAAT‑ATGTTCGTCTA‑CCAGTTGATTCGCCCCCTGGTC‑GTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCg                                                                     <  1:1763143/237‑1 (MQ=255)                                                                                                                                                                                       atgggggTGGAAAT‑AGGATCGTCTATCCAG‑TAATTCGCCCCCGGGTC‑GTAATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCa                                                          <  1:67955/236‑1 (MQ=255)                                                                                                                                                                                        tCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTACTTTTTGGt                                                                                                                                                                           >  1:254505/1‑125 (MQ=255)                                                                                                                                                                                        tCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCAc                                                         >  1:1287917/1‑239 (MQ=255)                                                                                                                                                                                        tCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCAc                                                         >  2:1256263/1‑239 (MQ=255)                                                                                                                                                                                        tCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCAc                                                         >  1:562989/1‑239 (MQ=255)                                                                                                                                                                                         cGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACg                                                        >  1:350528/1‑239 (MQ=255)                                                                                                                                                                                            ggTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAAt                                                            >  1:1762418/1‑232 (MQ=255)                                                                                                                                                                                               ggCAAT‑ATGATCGTCTATACAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTACATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCAGCCAGATAATCACGGTAt                                                    >  2:500004/1‑237 (MQ=255)                                                                                                                                                                                                     ttttttCGTCTACCCTG‑TGATCCGCCCCCGGGTC‑GTGATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCAt                                             <  1:1210534/234‑1 (MQ=255)                                                                                                                                                                                                       tGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgc                                                                                                                                                                                                 >  1:848966/1‑89 (MQ=255)                                                                                                                                                                                                        gATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa                                                                                                                                                                          <  1:1130270/111‑1 (MQ=255)                                                                                                                                                                                                        gATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa                                                                                                                                                                          >  2:1130270/1‑111 (MQ=255)                                                                                                                                                                                                           cctCTATCCTA‑TGATTCGCCCCCACGTC‑GTGATAATCACCGCCACCACCGCCTTTTTCTGCTGTCGCTCCGCCAGTAAACCAGGCGCGAACCAGTTTACTTTTTGGAAGGGTCTTGACGGAATACATTGTTATGAGGCTTGTAAGTCAAAGGGACTTACCCAGCCAATCAGAGGTGAAATATCACTGACGCCGGCCATGGCCTTGCCCGCCAGATGATCATGTCATTCTTCCTGGAttt                                       >  2:908950/3‑239 (MQ=255)                                                                                                                                                                                                             tCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAt                                         >  1:1168702/1‑235 (MQ=255)                                                                                                                                                                                                               tATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCa                                                                                                                                                        <  1:1112193/122‑1 (MQ=255)                                                                                                                                                                                                               tATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCa                                                                                                                                                        >  2:1112193/1‑122 (MQ=255)                                                                                                                                                                                                               tATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAat                                                                                                  >  1:1481430/1‑176 (MQ=255)                                                                                                                                                                                                               tATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTaa                                   >  2:1279052/1‑239 (MQ=255)                                                                                                                                                                                                                aTCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTaa                                   >  2:31639/1‑238 (MQ=255)                                                                                                                                                                                                                 tCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGc                                 >  1:1449590/1‑239 (MQ=255)                                                                                                                                                                                                                       tGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGtt                                                                                                                                                                                      >  1:288065/1‑85 (MQ=255)                                                                                                                                                                                                                       tGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAAGAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTAttc                                                  >  1:1692232/1‑217 (MQ=255)                                                                                                                                                                                                                       tGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATcc                                                                                              >  1:21842/1‑173 (MQ=255)                                                                                                                                                                                                                       tGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCa                                                          >  1:72694/1‑209 (MQ=255)                                                                                                                                                                                                                         aTTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtt                                                                                                                                                                                 >  1:483356/1‑88 (MQ=255)                                                                                                                                                                                                                         aTTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAcc                          >  1:1699044/1‑239 (MQ=255)                                                                                                                                                                                                                            cGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTa                            >  1:1056749/1‑234 (MQ=255)                                                                                                                                                                                                                            cGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGa                        >  1:1697982/1‑238 (MQ=255)                                                                                                                                                                                                                                  cTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTAttctt                                                >  1:1219839/1‑208 (MQ=255)                                                                                                                                                                                                                                                tCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTa    >  1:1484485/1‑239 (MQ=255)                                                                                                                                                                                                                                                  aCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAg  <  1:31639/239‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               ACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGG‑TCTGGCCGTTTTCTGC‑TTCGACTTCAACAACCAACGTGCCAAGAA‑CGTTAATGCCAAAT‑GACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAAT‑ATGATCGTCTATCCAG‑TGATTCGCCCCCTGGTC‑GTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAG  >  NC_007779/2365773‑2366245

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