Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,580	C→A	E134* (GAA→TAA)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,580	0	C	A	94.1%	146.5 / ‑1.9	52	E134* (GAA→TAA)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (24/24);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCA  >  NC_000913/789339‑789822                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    gaGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGGGCGCCGTACAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGCCCGCGACTTACAGTGAAgg                                                                                                                                                                                                                                            >  2:456721/2‑250 (MQ=255)  aGTTACTGTGGATTGTGACGACAGCCACCCCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAAGTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCTCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCg                                                                                                                                                                                                                                           <  2:504862/251‑1 (MQ=255)                 ggccgaccgctcacccgttGGGCATGGAAACTGCGTTGGTCCCAAGGGGGCGGAAATCGAACAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGAACGTTAAGCGCGATTTGTGCGTCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTGCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAACCg                                                                                                                                                                                                                            <  2:182781/234‑1 (MQ=255)                         gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGCAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTCACCCGGCTCCCTg                                                                                                                                                                                                                    >  2:324538/1‑251 (MQ=255)                          ccACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTACGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGATCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGc                                                                                                                                                                                                                   >  1:411900/1‑251 (MQ=255)                           cacaCCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCg                                                                                                                                                                                                                  <  1:115742/251‑1 (MQ=255)                                            aaaCTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATc                                                                                                                                                                                                 >  1:12381/1‑251 (MQ=255)                                             aaCTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGACCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATcc                                                                                                                                                                                                >  2:169290/1‑250 (MQ=255)                                                 gCTTTAGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGAGAAAGCTGCTGCAATACGGTTCCGACCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCAcc                                                                                                                                                                                              <  2:445349/249‑1 (MQ=255)                                                   tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGt                                                                                                                                                                                          >  1:490131/1‑251 (MQ=255)                                                       gTCCTCAGTGAGTGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCAc                                                                                                                                                                                     <  1:285876/251‑1 (MQ=255)                                                          cccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCAcgcc                                                                                                                                                                                  <  1:268876/251‑1 (MQ=255)                                                              gggAGCGGCAATAGAGCAGCAAGGCATGAGCTTTCTTGCCGAGCGCGTAAATTTGCTGATCCATGATCCCGCAGTGTCAGCCGACAAACTGGTTTGCTGCTACCTGACCGTGAAGTGCGATTTGTTCGCCGTCCAGCGGCAGAGGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGATAAGCGGAAGAACTTAACCCGGCACCCTGCGGGACATTGCCGCTCATCCCCATGT‑CCAGGccgccg                                                                                                                                                                              <  2:215396/249‑1 (MQ=255)                                                               tGAGCGGCAATCGATCAGCACGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGccgccg                                                                                                                                                                              >  2:140452/1‑250 (MQ=255)                                                                gtgtggCAATCGATCAGCAGCTCATGATTTTTCTTGCCGAGCGAGGAAATTACCTGCGCCATGATCCCGCGGTTCCAGCCTACAACCTGGTTTTCTGCTTCCTGACCGTTAAGTGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACATTGAAGCGGAATAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCCCCATGT‑CCACGCCGCCGa                                                                                                                                                                             <  2:349926/246‑1 (MQ=255)                                                                    ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCtg                                                                                                                                                                        <  2:462907/251‑1 (MQ=255)                                                                            aTCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTGCAGCCTACAAACTGGTTTTCTGCTTCCTGACCGGTAAGCGCGATTTGTGCGCCGTCAAGCGGCAGGTGCTAAAGCTGCTGCAATACGGTTCCAACCGCGACTTACAGTGACGCAGAAGAACtt                                                                                                                                                                                                                                  >  1:38767/1‑186 (MQ=255)                                                                                         gATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCCGCGGCAGCTGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGACGAACTTAACCCGGCACCCTGCGGCACCTTGCCGCTGATCACCATGT‑CCACGCCGCCCAAGCTGTTGTTACGCAGTTGCAAATgg                                                                                                                                                   >  2:335169/1‑250 (MQ=255)                                                                                              ttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTGCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCTCCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGttcca                                                                                                                                               <  2:376459/250‑4 (MQ=255)                                                                                                     cGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgccc                                                                                                                                       >  1:2753/1‑250 (MQ=255)                                                                                                         cgcgGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                                                                                                                                   <  1:289008/251‑1 (MQ=255)                                                                                                              aaaTTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAg                                                                                                                              <  1:351853/251‑1 (MQ=255)                                                                                                              aaaTTAGATGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGAGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGATCCGAACGCGACTTACAGTGAAGCGGAAGAACATAACACGGCACCCTGCGGCACATTGCCGCTGATCACCATGTACCA‑GCCGCCGAAGCTGTTGTTACGAAGGTACAGAGGTTTCACCACGCCACGAACga                                                                                                                                >  2:442614/1‑248 (MQ=255)                                                                                                                     cTGATCCATGATCCCGCAGTTACAGCCTACAACCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCCGATGTTTCACCACGCCACGAACGTAGTTAGccc                                                                                                                       >  2:493019/1‑251 (MQ=255)                                                                                                                      tGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGccc                                                                                                                       >  1:400242/1‑250 (MQ=255)                                                                                                                       gATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAc                                                                                                                                                                                                                               >  1:322999/1‑146 (MQ=255)                                                                                                                       gATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAc                                                                                                                                                                                                                               <  2:322999/146‑1 (MQ=255)                                                                                                                          ccATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCGTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTg                                                                                                                   <  1:414517/250‑1 (MQ=255)                                                                                                                           cATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCAATTGAt                                                                                                                 >  2:167223/1‑251 (MQ=255)                                                                                                                              gATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGt                                                                                                              >  1:224478/1‑251 (MQ=255)                                                                                                                              gATCCCGCAGTGACAGCCTACAAACTGGTTTTCTGCTTCATGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGt                                                                                                              <  1:221035/251‑1 (MQ=255)                                                                                                                                       gTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtttg                                                                                                     >  2:434205/1‑248 (MQ=255)                                                                                                                                         tACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGc                                                                                                    <  1:140566/250‑1 (MQ=255)                                                                                                                                                      cgggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCAt                                                                                      >  2:426772/3‑251 (MQ=255)                                                                                                                                                                    ccTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcg                                                                        <  2:143999/251‑1 (MQ=255)                                                                                                                                                                     cTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                                                                                                                                   <  2:156302/191‑1 (MQ=255)                                                                                                                                                                     cTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGaa                                                                                                                                   >  1:156302/1‑191 (MQ=255)                                                                                                                                                                     cTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgt                                                                       >  1:431024/1‑251 (MQ=255)                                                                                                                                                                       gACCGGTAAGCGCGATTGGGGCGCCGGCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGGTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtc                                                                      <  2:169496/250‑1 (MQ=255)                                                                                                                                                                              aaGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGAt                                                              <  1:387417/251‑1 (MQ=255)                                                                                                                                                                                    gATTTGGGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATCCGGTTCCGACCGCGACTTACAGGGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCAt                                                        <  2:148044/251‑1 (MQ=255)                                                                                                                                                                                     aTTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCAt                                                        >  1:430641/1‑250 (MQ=255)                                                                                                                                                                                           gcgcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGc                                                 >  1:133664/1‑251 (MQ=255)                                                                                                                                                                                                   cAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTACCCCGGCACCCGGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCa                                         <  1:169290/251‑1 (MQ=255)                                                                                                                                                                                                     gCGGCAGATGATAAACCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcg                                       <  1:426772/251‑1 (MQ=255)                                                                                                                                                                                                           gatgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACt                                 >  2:283878/1‑251 (MQ=255)                                                                                                                                                                                                               atAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTACCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTa                              <  1:314591/250‑1 (MQ=255)                                                                                                                                                                                                                     gtgctgcTATCCGGTTTCGACCTCGACTTCCAGTGAATCGGAAGCACTTATCCCGGCCCCCTGCCACACATTTCCGTTGATCACCATGG‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACTCCACGAACGTAGTTATCCCATTGATAGTTTTCAGGTGCGACAATGGGCGCATCGAGGGAAATCTCGTCGAGTTGAATTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCAt                       <  2:361389/250‑1 (MQ=255)                                                                                                                                                                                                                               aCGGTTCCGCCCGCGACTTACTGTGAAGCGGAAGAACTTAACCCGGCACCCTGTGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTCGTTACTCAGTTGCAGATTTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCa                                                                                                         <  1:358916/159‑1 (MQ=255)                                                                                                                                                                                                                               aCGGTTCCGACCGCGCCTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCAGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCCCCACGCCACGAACGTCGTTAACCCATTGATAGTTTTcc                                                                                                         >  2:358916/1‑158 (MQ=255)                                                                                                                                                                                                                                 ggTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCAGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACACCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGCTTTTCATAATCGACTGCCATCACGCGCACTTTACGATCATCGCGTGGTGcac           >  1:275192/1‑251 (MQ=255)                                                                                                                                                                                                                                          cgcgACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCAGCGCGTGGTGCACAACTGATCa  <  2:490131/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    AAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCA  >  NC_000913/789339‑789822

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