Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,252,059	C→T	P332S (CCG→TCG)	dhaR →	DNA‑binding transcriptional dual regulator DhaR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,252,059	0	C	T	95.7%	123.0 / ‑5.9	46	P332S (CCG→TCG)	dhaR	DNA‑binding transcriptional dual regulator DhaR
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (29/15);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGAT‑CGCACGGACAATGAAAATGGCCGTCT  >  NC_000913/1251812‑1252291                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          aaaacacactaCGTCACCCGCCTTAGTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGc                                                                                                                                                                                                                                        <  2:442529/241‑1 (MQ=255)            cGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACAGCCGCTCAAACACGTAGCAGCAACCTTTGAAAGCCAGCACCAGTTTATGGATGCGGTGATAACCATTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGCTGCGGCAGTTGAGGACCAGTCCATTATGATCAGTCAGCA‑ATA‑ACTTCGCTCTTAAGCCACAAGACGAATCGCAAACccgcccc                                                                                                                                                                                                                             >  2:155870/1‑249 (MQ=255)             gTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACACGTCAATTAGGAAACGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCcgccgcc                                                                                                                                                                                                                            >  2:421067/1‑250 (MQ=255)                                    cAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAggcg                                                                                                                                                                                                     >  2:6062/1‑251 (MQ=255)                                        aaaaCAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGCTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGCAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGgcgc                                                                                                                                                                                                 >  1:85364/1‑251 (MQ=255)                                                   atccgccaaaaCACGTAGAAGCAACCTTTGAAAGCCAGCACCAGATTATTGATGCGGTGAAAACACTTAAACCGATACTAGAAACGCAGGGAACAAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCCGTCAATTAGGAAAAGTCAGCC‑CTA‑CCTTCGCTCATATGCCACAGGAACATTCGCCACCCCGCCGCTTTATTCATTTTTGTCTCAAGGCGGCGCGGAGTAACttt                                                                                                                                                                                      >  2:170481/9‑251 (MQ=255)                                                   aTCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGCTGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCCGTAGCttt                                                                                                                                                                                      >  2:210866/1‑251 (MQ=255)                                                   aTCCGCTCAAACACGTAAAAGCAACCTTTGAAAGCCAGAACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCTTTTTGGTCGCCAGGCGGCGCGCAGTAGCttt                                                                                                                                                                                      >  1:139748/1‑251 (MQ=255)                                                         tCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtc                                                                                                                                                                                >  1:229274/1‑251 (MQ=255)                                                                          aCCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGAGCGCAGTAGCTTTCCTGTCCTGCTTTGTGgaga                                                                                                                                                                  >  2:53675/1‑248 (MQ=255)                                                                                aaaaGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGt                                                                                                                                                         <  1:409428/250‑1 (MQ=255)                                                                                      cAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCAAAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCa                                                                                                                                                    >  1:309006/1‑250 (MQ=255)                                                                                          caatttttattgtgcggagttcaaCCCTAAAACCGAAAAACGAAACGTGGTGAACCATCTTTATTTTGTTGCCCCCTCCGTTGTAATATATGCGGCATTTGATGACACTTCAATTAGGAAAATTCCCCC‑ACA‑CCTTCGCTCATATCCCACAGGACGTGCCTCAATCCCGCCGCTTGATTCATTTTGGTCGCCAGGCTGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGCGCGTGGGCAAGGcc                                                                                                                                               <  2:83506/229‑2 (MQ=255)                                                                                            cAGTTTATGAATGGGGTGATACCCCTAAACCCGATAAAAGAAACGCAGGAAACCAGCTTTATTTGGTTGCTCCACCCTGTGGAAAAGAGGCGGCGTTTGATGCCCAGTAAATTAGAAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGAAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACt                                                                                                                                             <  2:194166/251‑1 (MQ=255)                                                                                                  gttGATGAGGTGATAACCCTTAAACCGATAATAGAACCGCAGGTAACCAGCTTTATTTTGTTTCTCCATCCTGTGGAACAGATGCGGCAGTTGATCACCAGTCAATTAGGAAACGTCAGCC‑ATA‑CCTTCGCACATATGCCACAGGACGATTCGCAAACCAGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCGGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCagtgtga                                                                                                                                         >  2:59353/2‑243 (MQ=255)                                                                                                     gATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCCTCCTGTGGAACAGATGCGGCAGTTGATGACCTGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGTGAAGAGGGCGTGGGCCAGGCACTGCTAAGTcc                                                                                                                                    >  2:193739/1‑250 (MQ=255)                                                                                                           gTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCCGTTGCTGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTAGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAGGAGGGCGTGGGCAAGGCCCTGCTAAGTCAGGCagt                                                                                                                              >  1:55954/1‑249 (MQ=255)                                                                                                             gATAACCCTTAAACCGATAATAGAAACGCAGGGACCCGTCTTTATTTTGTTGCTCCATCTTGTGGAACAGATGCGTCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATC‑CCTTCGCTCATATGCCCCAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAAtt                                                                                                                             <  1:18596/250‑1 (MQ=255)                                                                                                             gATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATACTGTGGAACAGATGCGGCAGTTGATGACCTGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTAGAGAAGAGGGAGTGGGCAAGGCACTGCTAAGTCAGGCAATTc                                                                                                                            >  2:166118/1‑251 (MQ=255)                                                                                                               tAACCCTTAAACCGATAATAGGAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGTAACAGATGCGGCAGTTGATGACCAGTCAATTCGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTAATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCGTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCCGGCAATTCAt                                                                                                                          >  2:231381/1‑251 (MQ=255)                                                                                                                                ataGAAAGGCAGGGAACCAGCTTTATTTGGTTGCTCCATCCTGTGGAACAGATGCGGCAGTGGATGACCAGTCAATTAGGAAAAGTCAGCA‑ATA‑CCTTCGCTCATAGGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGtgc                                                                                                         <  1:210866/251‑1 (MQ=255)                                                                                                                                           gggAACCAGCTTTATTTTGTTTCCCCATCCGGTGGAACAGAGGCGGCAGTTGATGACTAGTCAATGAGAAAAAGTCAGCC‑AGC‑CCTTCGCTCAAATGCCACAGGACGTTTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCtt                                                                                              <  2:407409/251‑1 (MQ=255)                                                                                                                                              aaCCAGATTTATTTTGTTGATCCATCCTGTGGAACAGATGCGGCAGTTGATGAACAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCCTATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGGGGGCGTGGGCACGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTtat                                                                                            >  2:467427/1‑250 (MQ=255)                                                                                                                                               aCCAGCTTTATTTTGTTGCTACATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCACTTAGGCAACGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTAGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCATTGTGGAGAAGAGGGCGTGGGCAAGGCGATGCTAAGTCAGGAAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTatat                                                                                          >  2:420726/1‑251 (MQ=255)                                                                                                                                                 cAGCTTTATTTTGTTGCTCCATCCTGTGGACCAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTACGGCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCg                                                                                        >  2:72945/1‑251 (MQ=255)                                                                                                                                                    cTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCg                                                                                     >  2:474011/1‑251 (MQ=255)                                                                                                                                                      ttATTTTTTTGCTCCATCCGGGGAAACAGATGCGGCAGTGGATGCCCAGGCAATTAGGAAAAGTCACCC‑ATC‑CCTTCGCTCATATGCCAAAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGt                                                                                    <  1:438690/250‑1 (MQ=255)                                                                                                                                                                    cATCCTGTGGAACAGATGCGGCAGTTGATGCCCAGTCAATTAGGTAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGGGAAGAGGGCGGGGGCACGGCTCTGCTAAGTCAGGCAATTCATAGTGAAAGCGAGCGTGCTGCCGGTCCTTATATCGACGGCAATtgggtggtgat                                                                     >  2:151241/1‑242 (MQ=255)                                                                                                                                                                          gtgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGtg                                                                                                          >  2:398092/1‑208 (MQ=255)                                                                                                                                                                          gtgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGCCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGtg                                                                                                          <  1:398092/208‑1 (MQ=255)                                                                                                                                                                             gAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAg                                                                                                                                                  <  1:486277/165‑1 (MQ=255)                                                                                                                                                                             gAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAg                                                                                                                                                  >  2:486277/1‑165 (MQ=255)                                                                                                                                                                                 aGATGCGGCAGTTAATGACAAGACACTTAGGAAAATTCAGCC‑ATC‑CCTTCGCTCATATACCACAGGCCGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTgcg                                                        <  1:39945/251‑1 (MQ=255)                                                                                                                                                                                       ggCATTTGATGACCAGTCAATTAGGAAAAGTCAGC‑‑ATAACCTTCGCTCATATGCCACAGGACGATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCg                                                  <  1:124424/251‑1 (MQ=255)                                                                                                                                                                                        gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCgg                                                 >  1:406576/1‑251 (MQ=255)                                                                                                                                                                                              gatgaCCAGACAATTAGGAAAAGCAAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGgaagaa                                            <  1:6062/250‑1 (MQ=255)                                                                                                                                                                                                    cAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATtg                                     >  1:424436/1‑251 (MQ=255)                                                                                                                                                                                                      gTCAATTAGAAAAAGTCCGCCCAT‑‑CCTTTCCTCATAGCCCACAGGACGATTCGACACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATtggt                                   <  2:204221/251‑1 (MQ=255)                                                                                                                                                                                                             aGGAAAAGTTAGCC‑TTC‑CCTTCTCTCATACTCCACAAGACAATTCCCCACCCCGCCGCGTGATTCATTTTGGTCTCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGTCTCTGCTAAGTCAGGCAATTCATAATGAAATCGAGCGTGCTGCAGTTCCTTATTTCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCGGTCGGAAGAATTTATTGGTGGCGAt                             <  1:66991/250‑1 (MQ=255)                                                                                                                                                                                                                 aaaGTCAGCC‑ATA‑CAGTCGCACATATGCCACAGGACGATTCGCAAACACGCCGCTTGATTCATTTTGGCCGCCAGGCGGGGCGCAGTAGCTTTCCTGTCCTGCTTAGTGGAGCAGAGGGCGCGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTCTATGGTGATGCTGCGCTGGCGGAGGAATTTATTGGTTGCGAT‑CGccc                       >  2:149331/1‑249 (MQ=255)                                                                                                                                                                                                                     tCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGAT‑CGCACGGAc                   >  1:153778/1‑251 (MQ=255)                                                                                                                                                                                                                        catccatgccctCGCTTATACGCCACAGGTCGAATCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCATGCTGTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGGGCTGCAGGTCCTTATATCGCAGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATGTATTGGTGGCGAT‑CCCACGGACa                  <  2:55954/240‑1 (MQ=255)                                                                                                                                                                                                                          c‑ATA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAagcgagcg                                                                                                            >  1:223482/1‑158 (MQ=255)                                                                                                                                                                                                                             tA‑CCTTCGCTCATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGCGAAGAGGGCGTGGGCAAGGCACTGCTAAGTAAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGGTGCTGCGCTGGCGGAAGGATTTATTGGTGGAGATTCG‑ACGGACCATGaa             >  2:103478/1‑250 (MQ=255)                                                                                                                                                                                                                                        cATATGCCACGGGAGGATCCGAATACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGAT‑CGCACGGACAATGAAAATGGCCGTCt  <  2:229274/251‑1 (MQ=255)                                                                                                                                                                                                                                        cATATGCCACAGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGAT‑CGCACGGACAATGAAAATGGCCGTCt  >  1:303569/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          TGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCC‑ATA‑CCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGAT‑CGCACGGACAATGAAAATGGCCGTCT  >  NC_000913/1251812‑1252291

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