Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,649	0	T	G	32.8%	47.4 / 11.6	67	G65G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base G (14/6);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTT  >  NC_000913/433405‑433891                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    aaGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                                                                                                                                                                                                                                               <  2:158366/251‑1 (MQ=255)             aaTTTGGCGACGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTcac                                                                                                                                                                                                                                   >  1:127305/1‑250 (MQ=255)             aaTTTGGCGAAGAGATCGCGCGCGTGGAGGACTAAGCCGTGCAGGACGAGTATTACAGGGCGCGGGCGCAAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGGAAACCACAGGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGAGACCGCCTAT‑GTcaca                                                                                                                                                                                                                                  <  2:27633/251‑1 (MQ=255)              aTTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGCGCGTCATTGTCAAAGATGGCGAAATTGTCGGCGAAGGTTACCACCAACGTGCGGGCGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCACCTAT‑GTCacac                                                                                                                                                                                                                                 >  1:476569/1‑251 (MQ=255)                 tccggatctgaccgcgcgCCTGGAGGCCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTAT‑GTCACACt                                                                                                                                                                                                                                <  2:458244/237‑1 (MQ=255)                       agagATCGCGCGCCTGGAGGACTAAGACGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGGTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGGGACAAAGCAAAAAGTGGGCCCGcccagtaacccccatgacccc                                                                                                                                                                                                                         >  2:353532/1‑232 (MQ=255)                         agaccgcgcgCCTGGAGGACTAAGCCGTGCAGGAGGAGTATTACATGGCGCGGGCGCTAAGGATGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGGCGGTTGCGTAATTGTCAAAGATGGCGAAATTGTGGGTGAAGGTTACCACCAACGTGGGGGTAAACCACAGGCGGAAGTACACGCGTTGCTTATGGCGGGTCAAAAAGCCAAGGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGt                                                                                                                                                                                                                      <  2:27096/247‑1 (MQ=255)                          gaTCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCCAAGGGGCGACCGCCCAT‑GTCACACt                                                                                                                                                                                                                                >  1:458244/1‑240 (MQ=255)                          gaTCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAAGGAGGGGGGGGTAATGGTAAAAGAGGGCGAAATGGTCGGTGAAGGTTACAACAAACGAGGGGGGGAACCACATGCGGAAGAACAGGCGTGGGGAAGGGCGGGGGAAAAAGCCAAAGGGGGGACCGcgtaaggcaaacacgaaacctgt                                                                                                                                                                                                                       >  1:196705/1‑228 (MQ=255)                            ccgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGGGGCCGCTAAAGCTGGCGCAACGAGGACGTTTTCCCACGCATCCCAACCCGATTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACTTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGc                                                                                                                                                                                                                   <  1:359020/250‑1 (MQ=255)                              gcgcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGCGGGCCCGCCTAT‑TTCACCCTACAACCCTGTaagc                                                                                                                                                                                                                  >  1:178338/1‑247 (MQ=255)                                     tGGAGGACTAGGCCGAGCAGGAAGCGTATTACATGGCGCGGGAGCTAAAGCTGGCGCACCGAGGACGTTTTACCACGCATCCCAACCCGCATGTCGGGTGCGTCATTGTCACAGATGGAGAAATTGACGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGACGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGAGGGACCCCCTCT‑GTCACCCCCGCACCCTggtgcccacatagg                                                                                                                                                                                                          >  2:209322/1‑238 (MQ=255)                                         ggaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAGGGGCGACCGCCCAT‑GTCACACTCGAACCCTGTAGCCAACATGGGCGt                                                                                                                                                                                                       >  2:358152/1‑250 (MQ=255)                                               aGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTAT‑GTCCCACTCGAACCCTGTAGCCATCATGGTCGTACCCCAc                                                                                                                                                                                                >  1:245864/1‑251 (MQ=255)                                                gCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCCCCTAT‑GTCACCATCCAACCCTGTAGCCATCCTgggcggtcgccccc                                                                                                                                                                                               >  2:189924/1‑239 (MQ=255)                                                gCCGTGCAGGACGAGTATTACAGGGCGGGGGCGCAAAAGCTGGCGCAACGAGGACTTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGACGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAcc                                                                                                                                                                                               <  1:269976/251‑1 (MQ=255)                                                       aGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCCAT‑GTCACACACGAACCCCATAGCCATCCTGGGCGTACGCCACCgcgcggt                                                                                                                                                                                        >  2:416311/1‑245 (MQ=255)                                                                     aTGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCATTGATGGTGACCTCTTCGGAGAAGTTTTCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGCCCGCCCATTGCCA‑ACTCGGACCCTGTAGCCCGCCTGGTCGGACCCCACCGTGCtggtgcgcccccctccc                                                                                                                                                                          >  2:407318/1‑236 (MQ=255)                                                                      gggCGCGGGAGCTAAAGCTGGCGCAACGAGGTCGTTTTCCCACGCATCCCAACCCGAATGTCGGGGGCTTCTTTGTCAAAGATGGCGAAATGGTCGGTGAAGTTTACCCCCAACGGGCGGGTGACCCACATGCCGATGTACGCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCTACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCTTGCTGTGACGCACTCATcgcc                                                                                                                                                                         <  2:290565/250‑1 (MQ=255)                                                                            gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCCAAGGTGCGACCGCCTCT‑GTCACACTCGAACCCTAGAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTcctctcccccggc                                                                                                                                                                   >  2:235564/1‑239 (MQ=255)                                                                             gccgcaaaaGCTGGCGCAACGAGGACGTTTTACCGCGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGAGGGGGAAAGTGTCGGTGAAGGTTGCCACCAGCGTGCGGGTGAACCACATGCGGAAGAACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGTTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCg                                                                                                                                                                  <  1:31453/245‑1 (MQ=255)                                                                                    aaGCTGGCGAAAGGAGGTCGTTTTCGCACGCACCCCAGCCCGAAGGTCGGGCGCGTCATTGTCAAAGATGGCGAAATGGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCTGATGCCGAAGTGCACGGGTTGCGGAGGGCGGGTTAAAAAGCCAAAGGGGCAGCCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCTTAgcgcg                                                                                                                                                           <  2:288607/251‑1 (MQ=255)                                                                                           cgcAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTTTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCCCTCCTCGCCCCTGGCGTAGCGCGCGTGGtt                                                                                                                                                    >  1:137562/1‑251 (MQ=255)                                                                                                         tttACCACGCATCCCAACCCGATTGGCGGGTGCGTCATTGTAAAAGAGGGCGAAATTGTGGGTGAAGGTTACCCCCAACGTGGGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAg                                                                                                                                       <  2:256101/250‑1 (MQ=255)                                                                                                           gaCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCCACGTGCGGGTGAACCACATGCCGGAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCGACGGTGCGGACGCCTAT‑TTCCCACTCGAACCCTGTAGACATCCTGGGCGTACGCCACCGGGCTGGGACGAACTCATCGCAGCTGGCGTAGAGCCCGGGGTTGCCTCGacacccgatc                                                                                                                                    >  1:441844/2‑242 (MQ=255)                                                                                                              aacgcaaacaagcgcgaaggccggatGCGTCATTAACAAAGGAGGCGAAGTTGTAGGTAAAGGGTACGACCAACGTGCGGGTGAACCGCAGGCCGAGGTACAGGCGATGACAAGGCCGGGTGAAAAAGCCAGAGGGGCGAGCGCGTGG‑GGCACACGGGAGCCCTGTAGCGAGCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTa                                                                                                                                 <  1:407318/226‑1 (MQ=255)                                                                                                                      ccAACCCGCATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGACAAAGGTGCGACCGCCTAT‑GTCCCACTCGGACCCTGTAGCCAACATGGTCGTACGCCCCCGTGCTTTGACGCCCTCATCCCCGCTGGCGTAGCGCGCGTGGTTGCCTCGAGGCAACATCCCAACCCGc                                                                                                                           >  2:321334/1‑249 (MQ=255)                                                                                                                       cAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCCCCGCTGGCGTAGCGCGCGTGGTTTCCTCGATGCCAGATCCTAACCCGCAgg                                                                                                                        >  1:214924/1‑251 (MQ=255)                                                                                                                               aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTTACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTTCCTCGATGCAAGATCCTAACCCGCAGGTCGCTggg                                                                                                                >  1:372200/1‑251 (MQ=255)                                                                                                                               aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGAACGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCCTCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGCTGCAAGATCCTAACCCGCCGGTCCCTggg                                                                                                                >  2:185694/1‑251 (MQ=255)                                                                                                                                 gTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCCAT‑GTCAAACTCGAACCCTGTCGCCATCCTGGTCGTACGCCACCGTGCTGTGCAGCACTCATCGCCGCTGGCGTAGCGCGCCTGGTTGCCTCGATGCCATATCCCAACCCCCAGGGCGCCGGGCg                                                                                                              >  2:153050/1‑251 (MQ=255)                                                                                                                                       tGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCCTGGTAGTACGCCCCCGTGCTGTGACGCCCTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAACTCCACACCCGCAGGGCCCCGGGCGTGGACt                                                                                                        >  2:320556/1‑251 (MQ=255)                                                                                                                                       ggCGTCTTTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGAAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACt                                                                                                        <  1:118001/250‑1 (MQ=255)                                                                                                                                          gTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCCAAGGGGCGGCCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTTGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGGGTGGGCtttt                                                                                                     >  1:80771/1‑250 (MQ=255)                                                                                                                                             aTTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAGGGGCGACCGCCTAT‑GTCACAATCGAAACCTGTGGCCATCATGGTCGTACGCCACCGTGCTGTGGCGCACTCATCGCCGCTGGCGTAGCACGCGTGGGTGGCTCGGTGACAGATCCTAAAACGGAGATCGCTGGGCgggggctgtaca                                                                                                   >  2:204125/1‑239 (MQ=255)                                                                                                                                              ttGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGAGAAAAAGCCAAAGGGGCGACCACATAT‑GTCACACTCGAACCCTGTAGCCATAATGGTCGTACGCCAACGTGCTGTGACCCACTAATCGCCGCTGGCGTAGCGCGCCTGGTTGCCTCGATGCAAGATCCTAACCCGCGGGTCGCTGGGGGAGGACCTTACCGt                                                                                                 >  1:196074/1‑251 (MQ=255)                                                                                                                                                 tCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTAACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTTGGTGCAAGATCCTAACCCGCAGGTCGCTGGGCGGGGGCTTTACCGTCTg                                                                                              >  2:201234/1‑251 (MQ=255)                                                                                                                                                  cAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGc                                                                                             >  1:4276/1‑251 (MQ=255)                                                                                                                                                        tGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCCTCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                                                       >  1:308154/1‑251 (MQ=255)                                                                                                                                                                 ggTCGGGGAAGGTTACCACCACCGTGCGGGTGAACCACAGGCCGAAGTACCCGCGTGGCGTAGGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTg                                                                              <  2:278815/250‑1 (MQ=255)                                                                                                                                                                  gtgggtgaatgttccacccaACGGGGGGGTGAACCACATGCCGAAGTACCCGCGTGGCGTATGGGGGGGAAAAAAGCCAAAGGGGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTg                                                                              <  2:155584/233‑1 (MQ=255)                                                                                                                                                                      gTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACACTGTAGCCATAATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCACGATCCTAACCCGCAGGTCGCTGGGCGTGGAATTTACCGTCTGCAACATGCTGACATTGACGTc                                                                         >  1:356350/1‑251 (MQ=255)                                                                                                                                                                           ggTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTAT‑TTCACACTCGAACCCTGTAGCCCTCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCCTTGACGTCAGccc                                                                    >  1:292862/1‑250 (MQ=255)                                                                                                                                                                           ggTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCGCCTAT‑GTAACACCCGAACCCTGGAGCCCTCATGGTCGTACGCCACCGTGCTGTGACGCACTCCTCCCCGCTGGCGGAGCGCGCGTGGGTGCCTCGATGCAAGATTCTAACCCGCAGGGCGCTGGGGGTGGGCTTTACCGTCTGCAACAGGCTgggcttggcgtccgccc                                                                    >  2:28813/1‑236 (MQ=255)                                                                                                                                                                               accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGCCCGCCTAT‑GTCACACTCGACCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCCCTCATCGCCGCTGGCGTAGCTCGCGGGGTTGCCTCGGTGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCCACAGGCTGGAATTGAGGTCCGcccccgc                                                                >  2:433677/1‑246 (MQ=255)                                                                                                                                                                                 ccccagcGTGCGGGTGACCCGCAGGCGGAAGTACACGCGTTGCGGATGGCGGGTGAAAAAGCCAAGGGGGCGACCGCCTAT‑GTCACGCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCGCTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGAGGCAAGGTCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCt                                                              <  1:219395/245‑1 (MQ=255)                                                                                                                                                                                  accaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGTGCGACCGCCTAT‑GTCACACTCGGACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAACGCGCGTGGTTGCCTCGATGCCAGGTCCTAACCCGCAGGTCGCTGGGCGTGGGCTTTACCGTCTGCCACAGGCTGGCATTGATGTCAGCCACGGCCtg                                                             >  2:458582/1‑251 (MQ=255)                                                                                                                                                                                    caACGTGGGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTAAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgat                                                           <  1:185694/251‑1 (MQ=255)                                                                                                                                                                                      aCGTGCGGGTGAACCACATGCCGAAGTACACGCGCTGCGTATGGCGGGGGGAAAAGCCAAAGGTGCGGCCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCCGGGCGGGGGCTTTACCGTCTGCAAAAGGCCGGCATTTACGGCCGGCCCGGCCtgatgg                                                         >  2:13000/1‑250 (MQ=255)                                                                                                                                                                                             ggTGACCCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGATCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGa                                                  <  1:321334/251‑1 (MQ=255)                                                                                                                                                                                              gTGAACCACCTGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGGGGACTTTAACGTCTGCAACAGGCTGGCCTTGACGGCAGCCACGGCCTGATGATGAgaga                                                  >  1:145694/1‑247 (MQ=255)                                                                                                                                                                                               tGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAGGAGCGACCGCCTAT‑GTCCCACTCGAACCCTGTAGCCATCAAGGTCGTACCCCACCGTGCTGTGACGCCCTCATCGCCGCTGGCGTAGAGCGCGTGGTTGCATCGATGCAAGATCCTAACCCGCAGGTCGCTGGGGGTGGACTTTACCGGCTGCAACAGGCTGGCATTGTCGTCAGCCCCGGCCTGTGGATGAggga                                                  >  2:115840/1‑246 (MQ=255)                                                                                                                                                                                                         gccgaagtccccgcgtgtcgttgggCGGTTGAAAAAGCCAAAGGCGCGCCCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTTACGCACTCATCTCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATTCAATATCCTAGCCCGCAGTTCGCGGGGCGTGGACTTTGCCGTCTGCGTCAGGCTGGCATTGACGTCAGCTACGTGCTGATGATGAGGGAAGCCGAGCaa                                        <  2:254791/226‑1 (MQ=255)                                                                                                                                                                                                         gCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGGACCCTGTAGCCATCATGGGCGGACGCCACCGTGGCATGACGCACGCAACGCCGATGGCGTAGCGCGCGGGGGTGCCGCTATGCAAGATCCTAACCCGCCGGGCGCTGGGCGTGGGCTTTACAGTCTGCAACAGGCTGGCCTTGTCGGCAGCCACGGCCTGATGGTGAGTGCAGCCGGGCagtt                                      >  2:139686/1‑248 (MQ=255)                                                                                                                                                                                                             aaGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAAt                                  <  1:358152/251‑1 (MQ=255)                                                                                                                                                                                                               gTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGCGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTTTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCCGGTCGGTGGGCGTGGACTTTACCGTCTGCCACAGGCTGGCCTTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGGGCAATTGAATa                                 >  2:343823/1‑250 (MQ=255)                                                                                                                                                                                                                          gcttagggCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCa                     <  2:127305/245‑1 (MQ=255)                                                                                                                                                                                                                           ggaagggcggtggAAAAAGCCAATGGTGAGACCGCCTAT‑GTCACACTCGAACCCTGGAGGCATCATGGTCGGACGCCTCCGTGCTGTGACGCACGCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGGAGGTCGCTGGGCGTGGACTTTAGCGTCTGCAGCAGGCTGTCATTGACGTAAGCCATGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCCTGCTCaa                    <  2:23127/239‑1 (MQ=255)                                                                                                                                                                                                                           cGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCaa                    <  1:173304/251‑1 (MQ=255)                                                                                                                                                                                                                               cccgatctaaaaaaaCCAAGGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCg                 <  2:501093/241‑1 (MQ=255)                                                                                                                                                                                                                                  cGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATg             >  1:448314/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggTGAAACAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGTACCCTGTAGCCATCATGGTCGTACGACACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCCGTTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGCTGCGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgcg           >  2:64666/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggTGAAAAAGCCAAAGGAGAGACCGACTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGAGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGAAATTGAATAAAGGCTTTCTCAAGCGGagg             >  2:299233/1‑247 (MQ=255)                                                                                                                                                                                                                                       gAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCg                 >  1:501093/1‑242 (MQ=255)                                                                                                                                                                                                                                            agcaaaaGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCt   <  1:189924/247‑1 (MQ=255)                                                                                                                                                                                                                                            aGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCt   <  1:437835/251‑1 (MQ=255)                                                                                                                                                                                                                                             gcaaaaGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCtt  <  1:64666/248‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    AAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTT  >  NC_000913/433405‑433891

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