Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,875	T→G	S153S (TCT→TCG)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,875	0	T	G	98.7%	256.3 / NA	76	S153S (TCT→TCG)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (41/34);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGACATCTACGTTTCCCCT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTA  >  NC_000913/3966645‑3967107                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                         tgagaacttttcGTTTCCCCTAAGCCAATCCCGCCGTTTCAACCTCCGCACTGGTGATCCCATCTCGGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                                                                                                                                                                       <  1:595951/224‑1 (MQ=255)  atgaCATCTACGTTTCCCCT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGaa                                                                                                                                                                                                                                    >  1:128925/1‑236 (MQ=255)  atgaCATCTACGTTTCCCCT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGAACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGaa                                                                                                                                                                                                                                    >  1:221837/1‑236 (MQ=255)      cATCTACGTTTCCCCT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGATAACGCCCGCAACAAAATCCTCTTTGAGGACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACtgatgac                                                                                                                                                                                                                                 >  2:439745/1‑234 (MQ=255)      cATCTACGTTTCCCCT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTGACCCCGCTGCACGCAAACTCCCGTCTGCGAATGGAACGTGGTAACGGTTGGACTgcacat                                                                                                                                                                                                                                 >  2:595951/1‑230 (MQ=255)       aTCTACGTTTCCCCT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGAt                                                                                                                                                                                                                                 >  2:715575/1‑234 (MQ=255)                 ccccT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTcgc                                                                                                                                                                                                                     >  2:596787/1‑236 (MQ=255)                 ccccT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTcg                                                                                                                                                                                                                      >  1:551624/1‑235 (MQ=255)                       aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACt                                                                                                                                                                                                                >  1:624392/1‑236 (MQ=255)                                  tcttttCAACCGCCGTATGGGTAATCCATTCTCGGGAAAAGTTCGCCCGCCGAAAGACGGTGAACGTTATCTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                                                                                                                                                                     <  2:428812/233‑1 (MQ=255)                                  ctttcacccccTGCGCAGTGGTGTCACCATTTCGGAAAAGATACGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                                                                                                                                                                     <  2:241703/227‑1 (MQ=255)                                                 cACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTg                                                                                                                                                                                      >  2:650035/1‑236 (MQ=255)                                                     ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGGAACGGTGCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGAATCACCTATCGGTCGTGGTCa                                                                                                                                                                                  >  2:434959/1‑236 (MQ=255)                                                      gTGATACCAGCTCTGGCCAGATTCGCCCGCCGAAAGAATGTCAACGCTATTTTGCTCTGCTGAACCTTAACGTAGTTAACTTCGACAACCCTGAAAACGCCAGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGAACTGGCATCACCTATCGGACGTGGTCAg                                                                                                                                                                                 >  2:707890/1‑236 (MQ=255)                                                        gATACCTTCTCTTGTAAATTCCGCCCGCCGAAAGAAGGTTAACGCTATTTTGCGCTGCTTAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCg                                                                                                                                                                               <  2:305281/236‑1 (MQ=255)                                                         aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTc                                                                                                                                                                                   >  2:685757/1‑231 (MQ=255)                                                         aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTc                                                                                                                                                                                   <  1:685757/231‑1 (MQ=255)                                                         aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCg                                                                                                                                                                               >  2:517155/1‑235 (MQ=255)                                                           aCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCTTGGTAAGCg                                                                                                                                                                               >  1:688154/1‑233 (MQ=255)                                                           aCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCg                                                                                                                                                                               >  1:245038/1‑233 (MQ=255)                                                             cATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGt                                                                                                                                                                           <  1:650035/235‑1 (MQ=255)                                                              aTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCt                                                                                                                                                                         <  1:684893/236‑1 (MQ=255)                                                              aTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCt                                                                                                                                                                         <  1:596787/236‑1 (MQ=255)                                                                ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGa                                                                                                                                                                       >  1:315999/1‑236 (MQ=255)                                                                ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGa                                                                                                                                                                       >  1:596764/1‑236 (MQ=255)                                                                ctacgatAAGATCCCCCCGCCGAAAGAAGGTTAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGa                                                                                                                                                                       <  2:624392/230‑1 (MQ=255)                                                                     gTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                                                                                                                                                                   >  2:675397/1‑235 (MQ=255)                                                                               cccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                                                                                                                                                         >  2:204513/1‑235 (MQ=255)                                                                                     gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGc                                                                                                                                                                                                       >  1:156891/1‑183 (MQ=255)                                                                                     gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGc                                                                                                                                                                                                       <  2:156891/183‑1 (MQ=255)                                                                                              tGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGt                                                                                                                                                                           <  1:668499/202‑1 (MQ=255)                                                                                              tGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGt                                                                                                                                                                           >  2:668499/1‑202 (MQ=255)                                                                                                      aTTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgctg                                                                                                                                 <  1:517155/236‑1 (MQ=255)                                                                                                      aTTTCGCGCTGCTGAAAGTTAACGAAGTTACCTTCGACAAACCTGAAACCGCCCGCAACAAAATCCTCTTTGAGACCTTACCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgctg                                                                                                                                 <  2:479535/236‑1 (MQ=255)                                                                                                        tttGCGCTGCTGACAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACACAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGAGCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCa                                                                                                                               <  1:275045/236‑1 (MQ=255)                                                                                                            cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgcgtctgaac                                                                                                                           >  2:469485/1‑228 (MQ=255)                                                                                                             gctgctGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                                                                                                                          <  1:615424/236‑1 (MQ=255)                                                                                                             ccgtgctaaaTTTGACGGATTTAACTTCGAAAACCCTGAAAACGCCCGTAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCCAGCAACCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                                                                                                                          <  2:302146/229‑1 (MQ=255)                                                                                                                         tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCCCCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa                                                                                                              <  2:551624/236‑1 (MQ=255)                                                                                                                         tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa                                                                                                              <  1:715575/236‑1 (MQ=255)                                                                                                                         tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa                                                                                                              <  1:392334/236‑1 (MQ=255)                                                                                                                         tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa                                                                                                              <  1:204513/236‑1 (MQ=255)                                                                                                                           aCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCCGAGCATTGCttacaacaac                                                                                                            >  1:320066/1‑227 (MQ=255)                                                                                                                            cGAAGTTAACTTCGACAAACCTGAAACCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCAt                                                                                                             <  1:367587/234‑1 (MQ=255)                                                                                                                                  tAACTTCCAACACCCTGAAAACGTCCGAAACAAAATCCTTTTTGAAAACTTAACCCCGTTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACa                                                                                                     <  2:596764/236‑1 (MQ=255)                                                                                                                                           cAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCgg                                                                                            <  1:675397/236‑1 (MQ=255)                                                                                                                                                cTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACAATGCTGCTGCAGCACATTGCTCAGCGAATTGCTTACAACTACCCggc                                                                                           >  2:89736/1‑231 (MQ=255)                                                                                                                                                     aaCCCCCGCAACAAATTCCTCTTTGAGATCTTACCCCCGCTGCACGCAACCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACGGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCt                                                                                  <  2:363434/236‑1 (MQ=255)                                                                                                                                                          ccGCAACAATATCCTCTTGGAGAACTTACCCCCGCTGCAGGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATgg                                                                             <  2:688154/236‑1 (MQ=255)                                                                                                                                                          ccGCAACAAATGCCTTTTTGAAAACTTAACCCCGCTGCACGCAACCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATgg                                                                             <  2:245038/236‑1 (MQ=255)                                                                                                                                                          ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGGTGGCTGAt                                                                               >  1:332224/1‑234 (MQ=255)                                                                                                                                                                        tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAg                                                                                                                                                                                 <  1:12834/122‑1 (MQ=255)                                                                                                                                                                        tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAg                                                                                                                                                                                 >  2:12834/1‑122 (MQ=255)                                                                                                                                                                        tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgc                                                                                                                                      <  2:71380/165‑1 (MQ=255)                                                                                                                                                                        tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgc                                                                                                                                      >  1:71380/1‑165 (MQ=255)                                                                                                                                                                        tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcga                                                                >  1:383131/1‑235 (MQ=255)                                                                                                                                                                            tGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGaa                                                            >  2:623316/1‑235 (MQ=255)                                                                                                                                                                            tGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGaa                                                            >  1:546484/1‑235 (MQ=255)                                                                                                                                                                             gaaaaCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACg                                                          <  2:132210/233‑1 (MQ=255)                                                                                                                                                                                aaCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTcc                                                       >  2:56329/1‑236 (MQ=255)                                                                                                                                                                                    tAACCCAGCTGCACGCAACCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGga                                                    <  1:577461/235‑1 (MQ=255)                                                                                                                                                                                        cccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGaagaag                                               <  1:56329/236‑1 (MQ=255)                                                                                                                                                                                            cTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGCACGTCCGGAAGAAGTAAc                                           >  1:360468/1‑236 (MQ=255)                                                                                                                                                                                                 cGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAt                                                                                          >  1:676433/1‑184 (MQ=255)                                                                                                                                                                                                 cGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAt                                                                                          <  2:676433/184‑1 (MQ=255)                                                                                                                                                                                                 cGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCgaga                                      >  1:459456/1‑236 (MQ=255)                                                                                                                                                                                                         ctcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACAACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGACCTCCGGGACGAAGGAACGGGGATGCa                                  >  2:221758/1‑232 (MQ=255)                                                                                                                                                                                                               cgCGTACGGAACGTTGTAACGGTTCGGCTGAAGTATTAACTGCTAGGTGACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCCCCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTa                        <  2:185331/235‑1 (MQ=255)                                                                                                                                                                                                                           gTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtgt                                                                                     <  2:755029/163‑1 (MQ=255)                                                                                                                                                                                                                           gTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtgt                                                                                     >  1:755029/1‑163 (MQ=255)                                                                                                                                                                                                                           gTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCGGGTAAAAGGGGAAgtt            >  1:111046/1‑236 (MQ=255)                                                                                                                                                                                                                           gTGGTAACCGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTTTTCAGCGTGGTCTTATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCCGAGCATTGCTTACAACCACCCGGATTGTGGGCTGATGGTTCTGCTCATCGACGAACGTCCGGAAGAAGAAACCTATATGCAGCAACAGGTaaacatgaacgct            >  2:90778/1‑226 (MQ=255)                                                                                                                                                                                                                             ggTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgtt            >  1:724462/1‑234 (MQ=255)                                                                                                                                                                                                                             ggTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgtt            >  1:30554/1‑234 (MQ=255)                                                                                                                                                                                                                                    gTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCtt     >  2:43758/1‑234 (MQ=255)                                                                                                                                                                                                                                     ttCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTa  >  1:400878/1‑236 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                         GATGACATCTACGTTTCCCCT‑AGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTA  >  NC_000913/3966645‑3967107

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