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

AACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCG  >  NC_000913/734568‑735022                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACGGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGCCGGCGAAGACGGTGGCAGGGTGGTGGTCCCGCCGGGCCTGTTGCCGATGCCCGACCGGCGGGCAAGGGAAATCCTGCAGGGCCGGGGGGGGGGGGAAAgcacccccgg                                                                                                                                                                                                                           >  2:907697/1‑231 (MQ=11) aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGCCCAAGACGGTGGCAGAGTGGGGTTCCCGCCGGGGCTGGTGCTGCCGCTCGACCGGCTGGAAAATGGAAACCGGGCTGACCGGGGGGGGGCGGaaaagccccccga                                                                                                                                                                                                                           >  2:363263/1‑229 (MQ=11) aaCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGACATCCTGGCTGACCGGGGGGGGGAGGGAAGGcgccgcgg                                                                                                                                                                                                                           >  1:734316/1‑237 (MQ=9)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGGGGAAAGCCGCCCctggctgg                                                                                                                                                                                                                      >  1:237259/1‑239 (MQ=14)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggctg                                                                                                                                                                                                                       >  1:1035482/1‑238 (MQ=14)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGGGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGTGAGGTAGGAAAGCCGCCgctggccg                                                                                                                                                                                                                       >  1:637539/1‑236 (MQ=9)      cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGAAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAACGGGCGAGTGGGCGATGACGCCgcgggctgt                                                                                                                                                                                                                      >  1:590968/1‑238 (MQ=11)           cccGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTAGCAGTGTGGGGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCTGGAGAGGGCAAACCTGGCCGCCCGGGGGAGGGGGGGAAACGGCCgctggggggcgcgg                                                                                                                                                                                                                 >  2:226024/1‑230 (MQ=11)                aTCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGCAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGCGGTGTTCCCGCCGGAGCTGGCGCAGATGCTCGCACGGCTGGAAAGTGAAATCATGGCTGACCGGGTGAGGGAGGAAAGCCACCCCTGGCTGGCATGGTGCg                                                                                                                                                                                                             >  1:633109/1‑238 (MQ=11)                aTCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCCCGACCGGCCGGACAGCGCAATCCTGGCTGACCGGGGGAGGGGGAAAAACCCCCACGGGCCGGCaacacgcgg                                                                                                                                                                                                            >  2:452431/1‑231 (MQ=11)                                  ccGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGGTGCTCGACCGGCTGGACAGTGAAATCCTGGCTTCCCGGGGGGAGGGGGAAAACCGCCGCGGGCGCGCCTCGGGGGGCCTGAgggggtggaaagg                                                                                                                                                                                          >  2:341258/1‑229 (MQ=11)                                    ggggAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGGGCAGATGc                                                                                                                                                                                        >  1:183995/1‑239 (MQ=14)                                            tCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGCAAGGGAACTCCCGGCTGACCGGGGGGGGGAGGAAAGCCCCCCCTCGCTGGCACCGGGCGGCCTTGCGGGGGAGCCGGAGCAAAAACaa                                                                                                                                                                                >  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)                                                                                                                                                                                                        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)                                                                                                                                                                                                                    aCCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAAGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCAGCAGCCCCagaagcgt       >  2:624457/1‑237 (MQ=255)                                                                                                                                                                                                                          gtagtgggAAAAGCGCCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑TGGACCCGGTGTTCACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCg  <  1:907697/236‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           AACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGA‑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 < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑