Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	734,765	T→G	23.1%	pseudogene (633/1521 nt)	rhsO →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	734,765	0	T	G	23.1%	75.4 / 2.9	73	pseudogene (633/1521 nt)	rhsO	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (39/13);  new base G (16/4);  total (56/17)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATG  >  NC_000913/734531‑734988                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 gTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCCTCCGCAGTCCGGCGGCCGAGAACGGGGCCGTGCGGGGGTCCCCCCGGGCGCGGGGCAGACGGCCGAACGGCgggaaa‑gggaaa                                                                                                                                                                                                                              >  2:56948/1‑239 (MQ=11)      gCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCCGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGGTGTTCCCGCCGGGGCGGGGGCAAGGGGGCGCCCGGCTGGGAA‑GGGAAAAccgg                                                                                                                                                                                                                         >  2:841849/1‑237 (MQ=11)       cTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTg                                                                                                                                                                                                                         >  1:84614/1‑238 (MQ=31)       cTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGGGGGTCCCGCCGGGGCTGGGGCAGAGGCCCGACCGGCTGGAAA‑GTGGAATCCTg                                                                                                                                                                                                                         >  1:542050/1‑238 (MQ=11)          ggACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCCCAGCCTGGCGGCAACCCCTCAGCAGTCCGGCGGCGAAGACGGTGGCCGGGGGGGGTTCCCCCCGGGGCCGGTGCGGGGGCTCAGCCGGCTGGGAA‑AGGAAAACCCgggca                                                                                                                                                                                                                     >  2:1076441/1‑236 (MQ=11)             gggcccccGCCTCCCCACGATACAGACCGAAAACCCCCGCACCCAGACGTTTTA‑CCGGCGGGGAGGCTCACCCCCACTCACCAGGGTTGAAACCGCCACCGGTGAGCGGGCGAAACCGCAGCGCCGCAGCCTGGCGGATCCCCTTCAGCAGCCGGGCGGCGAAGACGGTGGGAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGAc                                                                                                                                                                                                                   <  2:994396/232‑1 (MQ=11)                cGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGGCCGGCGGCGAAGACGGTGGCAGGGTGGGGTTCCCGCCCGTGCTGGTGCCGATGCTCGACCGGCTGGAAA‑GTGAAAACCTGGCGGAACggg                                                                                                                                                                                                               >  2:645456/1‑239 (MQ=11)                     gCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGCCGGTGGCAGGGTGGTGTTCCCGCCGGGGCTGGGGCAGATGCTCGACCGGCTGGAAA‑GGGAAACCCTGTCTGCCCgggggggg                                                                                                                                                                                                          >  2:747947/1‑234 (MQ=11)                        tGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCggggggggggg                                                                                                                                                                                                       >  1:434282/1‑231 (MQ=14)                                    aGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTTAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACGGGCTGGAAA‑GTGAAAACCCGGCTGACCGGgtgcggtgggcaaggccgcgc                                                                                                                                                                                           >  1:742885/1‑221 (MQ=11)                                      aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACGGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGCCGGCGAAGACGGTGGCAGGGTGGTGGTCCCGCCGGGCCTGTTGCCGATGCCCGACCGGCGGGCAA‑GGGAAATCCTGCAGGGCCGGGGGGGGGGGGAAAgcacccccgg                                                                                                                                                                                         >  2:907697/1‑231 (MQ=11)                                      aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGCCCAAGACGGTGGCAGAGTGGGGTTCCCGCCGGGGCTGGTGCTGCCGCTCGACCGGCTGGAAA‑ATGGAAACCGGGCTGACCGGGGGGGGGCGGaaaagccccccga                                                                                                                                                                                         >  2:363263/1‑229 (MQ=11)                                      aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGACGACGGTGGCACTGTGGTGTTCCCGCCGGTGCCGGGGCCGATGCCCGACCGGCCTGAAA‑AGGAAACCCCGGCGGCCCggggaggtgggggaacgcgcccccg                                                                                                                                                                                         >  2:787705/1‑217 (MQ=11)                                      aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGACATCCTGGCTGACCGGGGGGGGGAGGGAAGGcgccgcgg                                                                                                                                                                                         >  1:734316/1‑237 (MQ=9)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGGGGAAAGCCGCCCctggctgg                                                                                                                                                                                    >  1:237259/1‑239 (MQ=14)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggctg                                                                                                                                                                                     >  1:1035482/1‑238 (MQ=14)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGGGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGTGAGGTAGGAAAGCCGCCgctggccg                                                                                                                                                                                     >  1:637539/1‑236 (MQ=9)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGGTGTGCCCGCCGGGGCTGGGGCGGAGGGCCGGCCGGGGGGAAC‑GTGAAATGCTGTCTGACCGGGGGAGTGAGGAAAgccccggggggcggg                                                                                                                                                                                    >  2:607388/1‑227 (MQ=11)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGAAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGAACGGGCGAGTGGGCGATGACGCCgcgggctgt                                                                                                                                                                                    >  1:590968/1‑238 (MQ=11)                                                cccGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTAGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCTGGAGA‑GGGCAAACCTGGCCGCCCGGGGGAGGGGGGGAAACGGCCgctggggggcgcgg                                                                                                                                                                               >  2:226024/1‑230 (MQ=11)                                                     aTCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGCGGTGTTCCCGCCGGAGCTGGCGCAGATGCTCGCACGGCTGGAAA‑GTGAAATCATGGCTGACCGGGTGAGGGAGGAAAGCCACCCCTGGCTGGCATGGTGCg                                                                                                                                                                           >  1:633109/1‑238 (MQ=11)                                                     aTCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCCGGACA‑GCGCAATCCTGGCTGACCGGGGGAGGGGGAAAAACCCCCACGGGCCGGCaacacgcgg                                                                                                                                                                          >  2:452431/1‑231 (MQ=11)                                                                       ccGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGGTGCTCGACCGGCTGGACA‑GTGAAATCCTGGCTTCCCGGGGGGAGGGGGAAAACCGCCGCGGGCGCGCCTCGGGGGGCCTGAgggggtggaaagg                                                                                                                                                        >  2:341258/1‑229 (MQ=11)                                                                         ggggAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGGGCAGATGc                                                                                                                                                      >  1:183995/1‑239 (MQ=14)                                                                                 tCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGCAA‑GGGAACTCCCGGCTGACCGGGGGGGGGAGGAAAGCCCCCCCTCGCTGGCACCGGGCGGCCTTGCGGGGGAGCCGGAGCAAAAACaa                                                                                                                                              >  2:593539/1‑238 (MQ=11)                                                                                                  ggTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGGGAGTGAGGCAAGCCACCCCTGGCTGGCATCGTGCGGCCTGACGGGGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGAac                                                                                                                              >  1:1026881/1‑238 (MQ=255)                                                                                                  ggTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGGGAGGGAGGGAAGCCGCCGCTGGCCGGCACCGTGCGGCCTGACGGGGGGGCAGATGAAAAAACAGAGGGGACCCgggtcaa                                                                                                                              >  2:912701/1‑234 (MQ=11)                                                                                                    ttGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGCGGTacac                                                                                                                            >  1:949982/1‑238 (MQ=16)                                                                                                     tGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGTGGGTGAGGCAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGGGGAGCAGATGCAAAACCAGA‑GGGACCCgggggaaccc                                                                                                                           >  1:772741/1‑230 (MQ=255)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGGAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCCTCGTGCGGCCTGACGGGGGGGCAGATGCAAAACCAGA‑AGGAACCGGGGCacaccgcg                                                                                                                        >  2:259228/1‑235 (MQ=11)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGTGCAGATGCAAAACCAGA‑GGGGCCCGGGGTACACGCCg                                                                                                                        >  1:36919/1‑239 (MQ=16)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCg                                                                                                                        >  1:747123/1‑239 (MQ=33)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGCGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTCCGGCCTGACGGTGGGGCAGATGCAAAACCAAA‑TGGACCCGGGTTACACCCCg                                                                                                                        >  1:588383/1‑239 (MQ=11)                                                                                                            gCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCGGGCATCGTCCGGCCTGACGGCGGAGCAGCTGCAAAACCAGA‑TGGACCCGgggtcaccaccgccgcg                                                                                                                   >  1:59289/1‑226 (MQ=11)                                                                                                            gCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGTGAGGGGGGAAAGCCGCCGCTGGCGGGCACCGTGGGGCCTGACGGTGGGGCAGAGGAAAAACCAAA‑TGGACCCGGGGGAGACACCgggggg                                                                                                                   >  2:287799/1‑235 (MQ=11)                                                                                                             ccACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGCGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAACTCCGGGCTGACCGGGTGCGGGCGGAAAGCCGCCGCTGCCTGCGATGTGGCGGCCGGACGCGGGGGCAGATGCAAAACCGAA‑GGGACCCGgggtcagccccggcggca                                                                                                                  >  2:605957/1‑225 (MQ=11)                                                                                                              cACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGGGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACCCCGACGCGaa                                                                                                                 >  1:646339/1‑239 (MQ=17)                                                                                                              cACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGCGGAGCAGATGCAAACCCAGA‑TGGACCCGGTGTACACGCCGACGCGaa                                                                                                                 >  1:785227/1‑239 (MQ=17)                                                                                                                           gCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGGCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCGGACCGGGGGAGTGAGGAACGCCCCCGCGGGGCGGCATCGTGCGGCCTGGCGGTGGAGCAGAGGCAAAACCAGA‑AGGGCCCGGGGGACACGCCGGCGCGAAAAAGCCaacggca                                                                                                    >  2:1095692/1‑233 (MQ=11)                                                                                                                                               ccTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGGGGGGCCTGACGGGGGAGCAGATGACAAACCAGA‑TGGACCCGGGGTACACGCCGGGGCGAGAAACCAACCTgacaaccgggcacacacgcggcc                                                                                >  2:391693/1‑217 (MQ=11)                                                                                                                                                         aCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgc                                                                      >  1:1105460/1‑239 (MQ=31)                                                                                                                                                                     tCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGAc                                                                                                                                       >  1:328552/1‑162 (MQ=31)                                                                                                                                                                     tCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCCCTGGCGCTTGTCAGCAc                                                          >  1:732065/1‑239 (MQ=18)                                                                                                                                                                     tCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCAGGCGCGACAAATCCACCTGTACCACTGCGACCATCGCGGCCCGCCGCTGGCGCTTGTcccgc                                                           >  1:862809/1‑234 (MQ=16)                                                                                                                                                                              gagggcGGGGGCAGTGGGGTGTTCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                                                  <  1:391693/233‑1 (MQ=16)                                                                                                                                                                              agaggtggtggcgggggtgttttCCCGCCGGTGCGGGTGCAGATGCTCGCCCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                                                  <  1:17781/217‑1 (MQ=11)                                                                                                                                                                              aaagcggggggtgggtgggtttccccGCCGGTGCTGGTGCAGATGCCCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                                                  <  1:841849/215‑1 (MQ=11)                                                                                                                                                                               aaGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCCCTGCGACCATCGCGGCCTGCCCCTGGCGCTTGTCAGCACGGAACGGGCa                                                >  1:912498/1‑239 (MQ=18)                                                                                                                                                                                      tGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGGAGATGCTAAACCCGA‑TGGCCCCGGGGTCCACGCCGGCGGGAAAAAGCCACCGGTACCACGGCACCCACCGCGGCCTGCCGGCGGAGCTTGGAACGCCCGAAGGGGCaaacggag                                         >  2:327232/1‑233 (MQ=11)                                                                                                                                                                                          ggggtggtgTCCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGAAG‑GTGAAACCCTGGCTGCCCGGGTGAGTGAGGAAACCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGt                                      <  1:259228/235‑1 (MQ=255)                                                                                                                                                                                            tgttgtGTCCCCCCGGGTTGTGGTGCAGTTGCTCGACAGGCTGGAAG‑GTGAAATCCCGGCTGACCGGATGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCCTCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                   <  2:434282/239‑1 (MQ=11)                                                                                                                                                                                            tgtgGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGGGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCACCACGGAAGGGGCAACAGAATGGGgc                                    >  1:637605/1‑238 (MQ=255)                                                                                                                                                                                             gtgattaaggggcaggTGCGGGGGCAGAGGCTCGACCGGCGGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                   <  1:401927/224‑1 (MQ=11)                                                                                                                                                                                             gggttgcccccGGCGGTGGTGGGGGAGGTGCCCGACGGGGTGGAAA‑GTAAAATCCTGGCGGCCGGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                   <  2:911409/230‑1 (MQ=11)                                                                                                                                                                                             ggggtgTTCCCGCGGGGTGTGGGGGAGATGCCCGGCCGGGTGAAAG‑GTGAAATCCGGGCTGACCGGGGGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                   <  1:593539/236‑1 (MQ=11)                                                                                                                                                                                              tggtgTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGgcgcg                                 >  1:836612/1‑238 (MQ=37)                                                                                                                                                                                                gtttccccGCCGGTGCTGGTGCAGATGCTCGCCCGGCTGGAAA‑GTGAAAGCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAg                                <  1:1076441/233‑1 (MQ=255)                                                                                                                                                                                                  gtcccggccgGTGCGGGGGGAGATGCTCGACCGGGTGGAAA‑GGGAAATCCTGGCTGACCGGGGGGGGGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGaa                              <  1:226024/232‑1 (MQ=11)                                                                                                                                                                                                  gtTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCCGA‑TGGACCCGGGGTCCACGCCGGCGCGAAAAATCCACCAGGACCACGGCGACCATCGGGGCCTGCCCCCGGGGCTTTTCAACACGGAACGGGCACAAGACTGGGCGGCAGaaa                             >  2:813905/1‑238 (MQ=11)                                                                                                                                                                                                  gtTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCCCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGGCGCAGAAt                             >  1:990779/1‑239 (MQ=37)                                                                                                                                                                                                     cccGCCGGGGTGGGGCCAGATGCCCGACGGGTGTGAAA‑GGGAAAGCCGGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACg                          <  1:481366/239‑1 (MQ=11)                                                                                                                                                                                                      ccgccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGACTGGTGCGCCGAATACGa                         >  1:922185/1‑239 (MQ=37)                                                                                                                                                                                                      ccgccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGCCCATCGCGGCCCGCCGCGGCCCCTTGTCAGCACGGAAAGGGCACCAGAATCGGGCGAAGAAAACGa                         >  2:292720/1‑239 (MQ=11)                                                                                                                                                                                                        gccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGCAGGGGCAACAGAATGTGGCGCAGAATACGATg                       >  1:625992/1‑239 (MQ=255)                                                                                                                                                                                                           ggTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGAAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAAt                    >  1:365277/1‑239 (MQ=37)                                                                                                                                                                                                             tGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACCCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCGCCCACTGGGGCTTGTCCGCCCGGAAGGGGAAACAAAAAGGGGGGCaaaaaaaaggaaaaga                  >  2:612099/1‑227 (MQ=11)                                                                                                                                                                                                                 ggggCAGATCCCCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGTGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCa              <  1:1095692/236‑1 (MQ=255)                                                                                                                                                                                                                    cacaaggccagACCGGGTGGAAA‑GGGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAAc            <  1:758404/228‑1 (MQ=255)                                                                                                                                                                                                                     gcgggggcgcacccgggggAAAGGTGAAACCCGGGGCGGCCGGGTGGGTGGGGAAAGCCGCGGCGGGCGGGCAGCGGGCGGCCGGACGGTGGAGAAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAAc            <  2:1105460/221‑1 (MQ=11)                                                                                                                                                                                                                       gATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCTTCGTGCGGCCTGACGGGGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGCCCTGCCCCTGGCCCGTGTCACCACCGAACGGGCCACAAAATGGTGCGCAAAATACGa                         >  2:401927/1‑222 (MQ=11)                                                                                                                                                                                                                       gATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCCTCGCGGCCTGCCCCTGGCGCTTGGCAGCACGGCAGGGGCAACAACACGGGGCGCGCAAAACGATGAATGGGGCAACcagc        >  2:763357/1‑239 (MQ=11)                                                                                                                                                                                                                         tGCTCGACCGGCTGGAAA‑GTGAAACCCTGGCTGACCGGGGGAGTGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAACCCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACctgc        <  1:1014149/237‑1 (MQ=255)                                                                                                                                                                                                                            tggccgggCGGGAAA‑GTGAAATCCGGGGTGGCCGGGGGGGGGGGGAAAGCCCCCGCGGGCTGGCATCGTGCGGCCTGACGGTGGAGAAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCtgaat   <  2:59289/233‑1 (MQ=11)                                                                                                                                                                                                                             cGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGGACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCGGCCGCTGGCGCTTGTCAGCCCGGAAGGGGCAACAGAAAGGTGCGCAGAAAACAATAAATGGGGCAACCCGCtgaatg  >  2:649008/1‑239 (MQ=11)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 GTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATG  >  NC_000913/734531‑734988

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