Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	720,627	A→T	7.8%	L15F (TTA→TTT)	ybfK →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	720,627	0	A	T	7.8%	174.2 / 4.4	51	L15F (TTA→TTT)	ybfK	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (17/30);  new base T (2/2);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATAAAACTATTGCCGCCACATCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAA  >  NC_000913/720345‑720914                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                gATAAAACTATTGCCGCCACATCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACag                                                                                                                                                                                                                                                                                             <  2:781304/287‑1 (MQ=255)           ttGCCGCCACATCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGc                                                                                                                                                                                                                                                                                    >  1:850000/1‑286 (MQ=255)               cgccACATCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                                                                                                                                                                                                                                                              <  2:172719/288‑1 (MQ=255)                   acaTCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                                                                                                                                                                                                                                                              <  2:567420/284‑1 (MQ=255)                   acaTCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                                                                                                                                                                                                                                                              >  1:567420/1‑284 (MQ=255)                   aaagccatatTACGACTTTCATTAATTTAGAAACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATa                                                                                                                                                                                                                                                                         <  2:47084/283‑1 (MQ=255)                    caTCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATAc                                                                                                                                                                                                                                                                        <  2:604086/289‑1 (MQ=255)                                        ttAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACa                                                                                                                                                                                                                                                    <  1:481206/289‑1 (MQ=255)                                          aaTGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACaga                                                                                                                                                                                                                                                  <  2:554652/289‑1 (MQ=255)                                          aaTGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACaga                                                                                                                                                                                                                                                  <  2:165140/289‑1 (MQ=255)                                                 aTACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAg                                                                                                                                                                                                                                           >  1:279032/1‑289 (MQ=255)                                                  tACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAgg                                                                                                                                                                                                                                          <  1:101809/289‑1 (MQ=255)                                                     aTTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTTACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAg                                                                                                                                                                                                                                       >  1:343743/1‑289 (MQ=255)                                                       tctcGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAgg                                                                                                                                                                                                                                          <  1:126273/284‑1 (MQ=255)                                                       tctcGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAgg                                                                                                                                                                                                                                          >  2:126273/1‑284 (MQ=255)                                                             tGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTAc                                                                                                                                                                                                                               >  1:285401/1‑289 (MQ=255)                                                                   gTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATAc                                                                                                                                                                                                                                                                        >  2:332464/1‑242 (MQ=255)                                                                   gTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATAc                                                                                                                                                                                                                                                                        <  1:332464/242‑1 (MQ=255)                                                                            tCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCatat                                                                                                                                                                                                                                                                >  2:755213/1‑241 (MQ=255)                                                                            tCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCatat                                                                                                                                                                                                                                                                <  1:755213/241‑1 (MQ=255)                                                                              cGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTTACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGt                                                                                                                                                                                                              <  2:343743/289‑1 (MQ=255)                                                                                  cAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCgg                                                                                                                                                                                                           <  1:196694/288‑1 (MQ=255)                                                                                          aCTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                                                                                                                                                                                                                                                              >  2:544938/1‑213 (MQ=255)                                                                                          aCTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                                                                                                                                                                                                                                                              <  1:544938/213‑1 (MQ=255)                                                                                                cGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCAt                                                                                                                                                                                            >  1:821054/1‑289 (MQ=255)                                                                                                                 aCATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTTACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCg                                                                                                                                                                           >  1:540959/1‑289 (MQ=255)                                                                                                                  cATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGa                                                                                                                                                                          <  2:850000/289‑1 (MQ=255)                                                                                                                                         gTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGt                                                                                                                                                                                                              <  1:537798/230‑1 (MQ=255)                                                                                                                                         gTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGt                                                                                                                                                                                                              >  2:537798/1‑230 (MQ=255)                                                                                                                                             tAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCa                                                                                                                                                                                                                                                                   >  2:222214/1‑173 (MQ=255)                                                                                                                                             tAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCa                                                                                                                                                                                                                                                                   <  1:222214/173‑1 (MQ=255)                                                                                                                                                   gTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTgaa                                                                                                                                                                                        >  1:220730/1‑242 (MQ=255)                                                                                                                                                   gTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTgaa                                                                                                                                                                                        <  2:220730/242‑1 (MQ=255)                                                                                                                                                   gTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAg                                                                                                                                         >  1:512539/1‑289 (MQ=255)                                                                                                                                                       tCAATTTTTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGt                                                                                                                                     <  2:170328/289‑1 (MQ=255)                                                                                                                                                       tCAATTGCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGt                                                                                                                                     <  1:434831/289‑1 (MQ=255)                                                                                                                                                            ttttcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGcc                                                                                                                                <  1:301841/286‑1 (MQ=255)                                                                                                                                                            ttCTCGAATGAGAAATCGACCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGcc                                                                                                                                <  1:235254/289‑1 (MQ=255)                                                                                                                                                                  aaTGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTatat                                                                                                                          >  1:36440/1‑289 (MQ=255)                                                                                                                                                                                cgTTCCGGAAAAAAGCAAAGTTACCCCGCTGGTAAATCAATCTGGTGGATTGGTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGGAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCgg                                                                                                            <  2:279032/289‑1 (MQ=255)                                                                                                                                                                                cgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCgg                                                                                                            <  1:33319/289‑1 (MQ=255)                                                                                                                                                                                             aGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTTACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCa                                                                                               <  2:540959/289‑1 (MQ=255)                                                                                                                                                                                                                    taaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagagag                                                                                                                                                                                                                                               >  1:685808/1‑122 (MQ=255)                                                                                                                                                                                                                    taaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagagag                                                                                                                                                                                                                                               <  2:685808/122‑1 (MQ=255)                                                                                                                                                                                                                            gTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCtt                                                                                                                              <  1:465042/227‑1 (MQ=255)                                                                                                                                                                                                                            gTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCtt                                                                                                                              >  2:465042/1‑227 (MQ=255)                                                                                                                                                                                                                                 ttGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGa                                                             <  1:380554/287‑1 (MQ=255)                                                                                                                                                                                                                                                          cTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTg                                  <  2:512539/289‑1 (MQ=255)                                                                                                                                                                                                                                                          cTTACCAGTCAGGTGGCGGATAATACAAGGGGTAACAGGGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTg                                  <  2:821054/289‑1 (MQ=255)                                                                                                                                                                                                                                                            tACAAGTCAGGTGGCTGATAATACAAGGGTTAAAAGAGAATGATGACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTAGGCTGAAGTGCTCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATTGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGaa                                <  2:579598/289‑1 (MQ=255)                                                                                                                                                                                                                                                                        ggCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGTTATTCTATCCAGGAAATACGGGTGAATTGTGAAGGTGa                    >  2:481887/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                           aaCAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATCCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCaa  >  2:665407/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                GATAAAACTATTGCCGCCACATCGATATTACGACTTTCATTAATGTAGATACATTCTCGCTGCGTGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAA  >  NC_000913/720345‑720914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑