Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	734,789	T→G	21.3%	pseudogene (657/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,789	0	T	G	21.3%	82.5 / 2.3	70	pseudogene (657/1521 nt)	rhsO	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (38/12);  new base G (14/6);  total (52/18)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCG  >  NC_000913/734568‑735022                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACGGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGCCGGCGAAGACGGTGGCAGGGTGGTGGTCCCGCCGGGCCTGTTGCCGATGCCCGACCGGCGGGCAA‑GGGAAATCCTGCAGGGCCGGGGGGGGGGGGAAAgcacccccgg                                                                                                                                                                                                                           >  2:907697/1‑231 (MQ=11) aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGCCCAAGACGGTGGCAGAGTGGGGTTCCCGCCGGGGCTGGTGCTGCCGCTCGACCGGCTGGAAA‑ATGGAAACCGGGCTGACCGGGGGGGGGCGGaaaagccccccga                                                                                                                                                                                                                           >  2:363263/1‑229 (MQ=11) aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGACATCCTGGCTGACCGGGGGGGGGAGGGAAGGcgccgcgg                                                                                                                                                                                                                           >  1:734316/1‑237 (MQ=9)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGGGGAAAGCCGCCCctggctgg                                                                                                                                                                                                                      >  1:237259/1‑239 (MQ=14)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggctg                                                                                                                                                                                                                       >  1:1035482/1‑238 (MQ=14)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGGGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GGGAAATCCTGGCTGACCGGGTGAGGTAGGAAAGCCGCCgctggccg                                                                                                                                                                                                                       >  1:637539/1‑236 (MQ=9)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGGTGTGCCCGCCGGGGCTGGGGCGGAGGGCCGGCCGGGGGGAAC‑GTGAAATGCTGTCTGACCGGGGGAGTGAGGAAAgccccggggggcggg                                                                                                                                                                                                                      >  2:607388/1‑227 (MQ=11)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGAAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGAACGGGCGAGTGGGCGATGACGCCgcgggctgt                                                                                                                                                                                                                      >  1:590968/1‑238 (MQ=11)           cccGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTAGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCTGGAGA‑GGGCAAACCTGGCCGCCCGGGGGAGGGGGGGAAACGGCCgctggggggcgcgg                                                                                                                                                                                                                 >  2:226024/1‑230 (MQ=11)                aTCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGCGGTGTTCCCGCCGGAGCTGGCGCAGATGCTCGCACGGCTGGAAA‑GTGAAATCATGGCTGACCGGGTGAGGGAGGAAAGCCACCCCTGGCTGGCATGGTGCg                                                                                                                                                                                                             >  1:633109/1‑238 (MQ=11)                aTCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCCGGACA‑GCGCAATCCTGGCTGACCGGGGGAGGGGGAAAAACCCCCACGGGCCGGCaacacgcgg                                                                                                                                                                                                            >  2:452431/1‑231 (MQ=11)                                  ccGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGGTGCTCGACCGGCTGGACA‑GTGAAATCCTGGCTTCCCGGGGGGAGGGGGAAAACCGCCGCGGGCGCGCCTCGGGGGGCCTGAgggggtggaaagg                                                                                                                                                                                          >  2:341258/1‑229 (MQ=11)                                    ggggAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGGGCAGATGc                                                                                                                                                                                        >  1:183995/1‑239 (MQ=14)                                            tCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGCAA‑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)                                                                                                                                                                                                         aaaTCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGACGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCAGCCTCa                    >  1:393660/1‑238 (MQ=255)                                                                                                                                                                                                            tCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCt                 >  1:278543/1‑238 (MQ=37)                                                                                                                                                                                                            tCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTg                >  1:462740/1‑239 (MQ=37)                                                                                                                                                                                                               tggttggcggggtgggtGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag             <  1:649008/230‑1 (MQ=255)                                                                                                                                                                                                               tGGCTGACGGGGGGGGGGGGGAAGCCCGCCGCTGGCTGCAATCGTGGGGCCTGACGGGGGAGAAGAGGCAAAACCAGA‑TGGCCCCGGTGTACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag             <  2:462740/239‑1 (MQ=11)                                                                                                                                                                                                               tGGCTGACCGGGTGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag             <  1:763357/239‑1 (MQ=37)                                                                                                                                                                                                                 gCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGGACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTTGCAGCACTGAAGGGGCAAAAGAATGGGGCGCAGAAACCGCTAAATGGGGCAACCGGCTAAATAAAAAAAAACCgaaaaaaacgcgaga           >  2:498667/1‑225 (MQ=11)                                                                                                                                                                                                                     aCCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAAGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCAGCAGCCCCagaagcgt       >  2:624457/1‑237 (MQ=255)                                                                                                                                                                                                                           gtagtgggAAAAGCGCCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTTCACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCg  <  1:907697/236‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           AACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA‑GTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCG  >  NC_000913/734568‑735022

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