Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	734,765	T→G	25.9%	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	25.9%	63.3 / 3.1	61	pseudogene (633/1521 nt)	rhsO	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (33/9);  new base G (14/4);  total (48/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATG  >  NC_000913/734531‑734988                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 gTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCCTCCGCAGTCCGGCGGCCGAGAACGGGGCCGTGCGGGGGTCCCCCCGGGCGCGGGGCAGACGGCCGAACGGCgggaaagggaaa                                                                                                                                                                                                                              >  2:56948/1‑239 (MQ=11)      gCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCCGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGGTGTTCCCGCCGGGGCGGGGGCAAGGGGGCGCCCGGCTGGGAAGGGAAAAccgg                                                                                                                                                                                                                         >  2:841849/1‑237 (MQ=11)       cTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                                                                                                                                                                                         >  1:84614/1‑238 (MQ=31)       cTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGGGGGTCCCGCCGGGGCTGGGGCAGAGGCCCGACCGGCTGGAAAGTGGAATCCTg                                                                                                                                                                                                                         >  1:542050/1‑238 (MQ=11)          ggACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCCCAGCCTGGCGGCAACCCCTCAGCAGTCCGGCGGCGAAGACGGTGGCCGGGGGGGGTTCCCCCCGGGGCCGGTGCGGGGGCTCAGCCGGCTGGGAAAGGAAAACCCgggca                                                                                                                                                                                                                     >  2:1076441/1‑236 (MQ=11)             gggcccccGCCTCCCCACGATACAGACCGAAAACCCCCGCACCCAGACGTTTTA‑CCGGCGGGGAGGCTCACCCCCACTCACCAGGGTTGAAACCGCCACCGGTGAGCGGGCGAAACCGCAGCGCCGCAGCCTGGCGGATCCCCTTCAGCAGCCGGGCGGCGAAGACGGTGGGAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAc                                                                                                                                                                                                                   <  2:994396/232‑1 (MQ=11)                cGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGGCCGGCGGCGAAGACGGTGGCAGGGTGGGGTTCCCGCCCGTGCTGGTGCCGATGCTCGACCGGCTGGAAAGTGAAAACCTGGCGGAACggg                                                                                                                                                                                                               >  2:645456/1‑239 (MQ=11)                     gCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGCCGGTGGCAGGGTGGTGTTCCCGCCGGGGCTGGGGCAGATGCTCGACCGGCTGGAAAGGGAAACCCTGTCTGCCCgggggggg                                                                                                                                                                                                          >  2:747947/1‑234 (MQ=11)                        tGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCggggggggggg                                                                                                                                                                                                       >  1:434282/1‑231 (MQ=14)                                      aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACGGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGCCGGCGAAGACGGTGGCAGGGTGGTGGTCCCGCCGGGCCTGTTGCCGATGCCCGACCGGCGGGCAAGGGAAATCCTGCAGGGCCGGGGGGGGGGGGAAAgcacccccgg                                                                                                                                                                                         >  2:907697/1‑231 (MQ=11)                                      aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGCCCAAGACGGTGGCAGAGTGGGGTTCCCGCCGGGGCTGGTGCTGCCGCTCGACCGGCTGGAAAATGGAAACCGGGCTGACCGGGGGGGGGCGGaaaagccccccga                                                                                                                                                                                         >  2:363263/1‑229 (MQ=11)                                      aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGACATCCTGGCTGACCGGGGGGGGGAGGGAAGGcgccgcgg                                                                                                                                                                                         >  1:734316/1‑237 (MQ=9)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGGGGAAAGCCGCCCctggctgg                                                                                                                                                                                    >  1:237259/1‑239 (MQ=14)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggctg                                                                                                                                                                                     >  1:1035482/1‑238 (MQ=14)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGGGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGTGAGGTAGGAAAGCCGCCgctggccg                                                                                                                                                                                     >  1:637539/1‑236 (MQ=9)                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGAAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAACGGGCGAGTGGGCGATGACGCCgcgggctgt                                                                                                                                                                                    >  1:590968/1‑238 (MQ=11)                                                cccGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTAGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCTGGAGAGGGCAAACCTGGCCGCCCGGGGGAGGGGGGGAAACGGCCgctggggggcgcgg                                                                                                                                                                               >  2:226024/1‑230 (MQ=11)                                                     aTCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGCGGTGTTCCCGCCGGAGCTGGCGCAGATGCTCGCACGGCTGGAAAGTGAAATCATGGCTGACCGGGTGAGGGAGGAAAGCCACCCCTGGCTGGCATGGTGCg                                                                                                                                                                           >  1:633109/1‑238 (MQ=11)                                                     aTCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCCGGACAGCGCAATCCTGGCTGACCGGGGGAGGGGGAAAAACCCCCACGGGCCGGCaacacgcgg                                                                                                                                                                          >  2:452431/1‑231 (MQ=11)                                                                       ccGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGGTGCTCGACCGGCTGGACAGTGAAATCCTGGCTTCCCGGGGGGAGGGGGAAAACCGCCGCGGGCGCGCCTCGGGGGGCCTGAgggggtggaaagg                                                                                                                                                        >  2:341258/1‑229 (MQ=11)                                                                         ggggAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGGGCAGATGc                                                                                                                                                      >  1:183995/1‑239 (MQ=14)                                                                                 tCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGCAAGGGAACTCCCGGCTGACCGGGGGGGGGAGGAAAGCCCCCCCTCGCTGGCACCGGGCGGCCTTGCGGGGGAGCCGGAGCAAAAACaa                                                                                                                                              >  2:593539/1‑238 (MQ=11)                                                                                                  ggTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGCAAGCCACCCCTGGCTGGCATCGTGCGGCCTGACGGGGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGAac                                                                                                                              >  1:1026881/1‑238 (MQ=255)                                                                                                  ggTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGGGAGGGAGGGAAGCCGCCGCTGGCCGGCACCGTGCGGCCTGACGGGGGGGCAGATGAAAAAACAGAGGGGACCCgggtcaa                                                                                                                              >  2:912701/1‑234 (MQ=11)                                                                                                    ttGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGCGGTacac                                                                                                                            >  1:949982/1‑238 (MQ=16)                                                                                                     tGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGTGGGTGAGGCAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGGGGAGCAGATGCAAAACCAGA‑GGGACCCgggggaaccc                                                                                                                           >  1:772741/1‑230 (MQ=255)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGGAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCCTCGTGCGGCCTGACGGGGGGGCAGATGCAAAACCAGA‑AGGAACCGGGGCacaccgcg                                                                                                                        >  2:259228/1‑235 (MQ=11)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGTGCAGATGCAAAACCAGA‑GGGGCCCGGGGTACACGCCg                                                                                                                        >  1:36919/1‑239 (MQ=16)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCg                                                                                                                        >  1:747123/1‑239 (MQ=33)                                                                                                       aaaCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGCGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTCCGGCCTGACGGTGGGGCAGATGCAAAACCAAA‑TGGACCCGGGTTACACCCCg                                                                                                                        >  1:588383/1‑239 (MQ=11)                                                                                                            gCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGTGAGGGGGGAAAGCCGCCGCTGGCGGGCACCGTGGGGCCTGACGGTGGGGCAGAGGAAAAACCAAA‑TGGACCCGGGGGAGACACCgggggg                                                                                                                   >  2:287799/1‑235 (MQ=11)                                                                                                              cACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGGGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACCCCGACGCGaa                                                                                                                 >  1:646339/1‑239 (MQ=17)                                                                                                              cACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGCGGAGCAGATGCAAACCCAGA‑TGGACCCGGTGTACACGCCGACGCGaa                                                                                                                 >  1:785227/1‑239 (MQ=17)                                                                                                                           gCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGGCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCGGACCGGGGGAGTGAGGAACGCCCCCGCGGGGCGGCATCGTGCGGCCTGGCGGTGGAGCAGAGGCAAAACCAGA‑AGGGCCCGGGGGACACGCCGGCGCGAAAAAGCCaacggca                                                                                                    >  2:1095692/1‑233 (MQ=11)                                                                                                                                                         aCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgc                                                                      >  1:1105460/1‑239 (MQ=31)                                                                                                                                                                     tCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGAc                                                                                                                                       >  1:328552/1‑162 (MQ=31)                                                                                                                                                                     tCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCCCTGGCGCTTGTCAGCAc                                                          >  1:732065/1‑239 (MQ=18)                                                                                                                                                                     tCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCAGGCGCGACAAATCCACCTGTACCACTGCGACCATCGCGGCCCGCCGCTGGCGCTTGTcccgc                                                           >  1:862809/1‑234 (MQ=16)                                                                                                                                                                              gagggcGGGGGCAGTGGGGTGTTCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                                                  <  1:391693/233‑1 (MQ=16)                                                                                                                                                                               aaGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCCCTGCGACCATCGCGGCCTGCCCCTGGCGCTTGTCAGCACGGAACGGGCa                                                >  1:912498/1‑239 (MQ=18)                                                                                                                                                                                          ggggtggtgTCCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGAAGGTGAAACCCTGGCTGCCCGGGTGAGTGAGGAAACCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGt                                      <  1:259228/235‑1 (MQ=255)                                                                                                                                                                                            tgttgtGTCCCCCCGGGTTGTGGTGCAGTTGCTCGACAGGCTGGAAGGTGAAATCCCGGCTGACCGGATGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCCTCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                   <  2:434282/239‑1 (MQ=11)                                                                                                                                                                                            tgtgGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGGGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCACCACGGAAGGGGCAACAGAATGGGgc                                    >  1:637605/1‑238 (MQ=255)                                                                                                                                                                                             gggttgcccccGGCGGTGGTGGGGGAGGTGCCCGACGGGGTGGAAAGTAAAATCCTGGCGGCCGGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                   <  2:911409/230‑1 (MQ=11)                                                                                                                                                                                             ggggtgTTCCCGCGGGGTGTGGGGGAGATGCCCGGCCGGGTGAAAGGTGAAATCCGGGCTGACCGGGGGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                   <  1:593539/236‑1 (MQ=11)                                                                                                                                                                                              tggtgTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGgcgcg                                 >  1:836612/1‑238 (MQ=37)                                                                                                                                                                                                gtttccccGCCGGTGCTGGTGCAGATGCTCGCCCGGCTGGAAAGTGAAAGCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAg                                <  1:1076441/233‑1 (MQ=255)                                                                                                                                                                                                  gtcccggccgGTGCGGGGGGAGATGCTCGACCGGGTGGAAAGGGAAATCCTGGCTGACCGGGGGGGGGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGaa                              <  1:226024/232‑1 (MQ=11)                                                                                                                                                                                                  gtTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCCGA‑TGGACCCGGGGTCCACGCCGGCGCGAAAAATCCACCAGGACCACGGCGACCATCGGGGCCTGCCCCCGGGGCTTTTCAACACGGAACGGGCACAAGACTGGGCGGCAGaaa                             >  2:813905/1‑238 (MQ=11)                                                                                                                                                                                                  gtTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCCCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGGCGCAGAAt                             >  1:990779/1‑239 (MQ=37)                                                                                                                                                                                                     cccGCCGGGGTGGGGCCAGATGCCCGACGGGTGTGAAAGGGAAAGCCGGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACg                          <  1:481366/239‑1 (MQ=11)                                                                                                                                                                                                      ccgccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGACTGGTGCGCCGAATACGa                         >  1:922185/1‑239 (MQ=37)                                                                                                                                                                                                      ccgccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGCCCATCGCGGCCCGCCGCGGCCCCTTGTCAGCACGGAAAGGGCACCAGAATCGGGCGAAGAAAACGa                         >  2:292720/1‑239 (MQ=11)                                                                                                                                                                                                        gccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGCAGGGGCAACAGAATGTGGCGCAGAATACGATg                       >  1:625992/1‑239 (MQ=255)                                                                                                                                                                                                           ggTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGAAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAAt                    >  1:365277/1‑239 (MQ=37)                                                                                                                                                                                                                 ggggCAGATCCCCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGTGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCa              <  1:1095692/236‑1 (MQ=255)                                                                                                                                                                                                                    cacaaggccagACCGGGTGGAAAGGGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAAc            <  1:758404/228‑1 (MQ=255)                                                                                                                                                                                                                       gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCTTCGTGCGGCCTGACGGGGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGCCCTGCCCCTGGCCCGTGTCACCACCGAACGGGCCACAAAATGGTGCGCAAAATACGa                         >  2:401927/1‑222 (MQ=11)                                                                                                                                                                                                                       gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCCTCGCGGCCTGCCCCTGGCGCTTGGCAGCACGGCAGGGGCAACAACACGGGGCGCGCAAAACGATGAATGGGGCAACcagc        >  2:763357/1‑239 (MQ=11)                                                                                                                                                                                                                         tGCTCGACCGGCTGGAAAGTGAAACCCTGGCTGACCGGGGGAGTGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAACCCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACctgc        <  1:1014149/237‑1 (MQ=255)                                                                                                                                                                                                                            tggccgggCGGGAAAGTGAAATCCGGGGTGGCCGGGGGGGGGGGGAAAGCCCCCGCGGGCTGGCATCGTGCGGCCTGACGGTGGAGAAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCtgaat   <  2:59289/233‑1 (MQ=11)                                                                                                                                                                                                                             cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGGACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCGGCCGCTGGCGCTTGTCAGCCCGGAAGGGGCAACAGAAAGGTGCGCAGAAAACAATAAATGGGGCAACCCGCtgaatg  >  2:649008/1‑239 (MQ=11)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 GTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCA‑CGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑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 < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑