Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	433,649	T→G	32.7%	G65G (GGT→GGG)	ribD →	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,649	0	T	G	32.7%	32.0 / 5.7	44	G65G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (10/18);  new base G (12/3);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 9.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGC  >  NC_000913/433419‑433883                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCAAA‑GAGGCGAccc                                                                                                                                                                                                                                       >  1:193026/1‑237 (MQ=255)               cgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTTTGGC‑GGGGGGGAAAACCAAA‑AGGGGGACCCCCCAT‑GTccaacccg                                                                                                                                                                                                                       >  1:882922/1‑232 (MQ=255)                     ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCCAA‑AGGGCGACCCCCCAT‑GTCACCCTCGAAccc                                                                                                                                                                                                                  >  1:67511/1‑238 (MQ=255)                              cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCAAA‑GGTGCGACCGCCCCT‑GTCAACCTCGAACCCTCTAGCcatc                                                                                                                                                                                                        >  1:41358/1‑239 (MQ=255)                               gatgCCGTGCGGGAGTAGTATGACATGGCGCGGGCTCAAAAGGGTGCCCACGGTGTACTTTGTACCCCCCATCCAACCACGAAGGTGGGTGGCTTCATTGTAAAAGAGGGCAAATTTGCCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATCCCGAAGTCCACCC‑TTTGCGTAGGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                                                       <  2:41358/236‑1 (MQ=255)                                 aGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCCAA‑AGGGCGACCGCCCAT‑GTCACACTCCAACCCCGTAGCCATCATg                                                                                                                                                                                                     >  1:214080/1‑239 (MQ=255)                                     gTGCAGGACGAGTTTTACAGGGCGGGGGCGATAAAGCGGGCGCAACGAGGACGTTTTCCCACGCATCCCACCCCGAATGTCGGGTGCTTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGC‑GTTGCGTAGGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCg                                                                                                                                                                                                 <  1:1177589/239‑1 (MQ=255)                                      tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑AGGGCGACCGCCTAT‑GTCACACTCGAACCCCGTAGCCCCCATcgtcgt                                                                                                                                                                                                >  1:1284265/1‑239 (MQ=255)                                      tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCAAA‑AGGGCGACCGCCCATGGTCAC‑CTCGCAACCCGTAGGCCACATGGTCGt                                                                                                                                                                                                >  1:1199913/1‑239 (MQ=255)                                                gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCCAT‑GTCACACTCGAACCCTGTAGCCATCATGGGCGTACGCCACCGt                                                                                                                                                                                      >  1:190422/1‑239 (MQ=255)                                                gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAACCAAA‑GGGGGGGCCCCCTCT‑TTCACACCCGCACCCCGGTGCCACCCTGGGCGTAcccccc                                                                                                                                                                                         >  1:576169/1‑234 (MQ=255)                                                 tATCACATGGCGCGGGCGCTAAAGCTGGCGCAAGGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTca                                                                                                                                                                                                                             <  1:599955/199‑1 (MQ=255)                                                      cgtgggcggggCGCTAAAGGTGGCGCAACGAGGACTTTTTCCCGCGCATCCCAACCCGAATGGCGGGTGCGTAATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCtgtg                                                                                                                                                                                <  1:863996/231‑1 (MQ=255)                                                         ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCAAA‑GGGGCGACCGCCTTT‑TTCACACTCGAACCCCGTAGCCACCATGGGCTTACGCCACCCGCCTGTGAc                                                                                                                                                                              >  1:924888/1‑238 (MQ=255)                                                           cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCAAA‑GGGGCGAACGCCTAT‑GTCACACTCGAACCCTTTAGCCAACATGGGCGTACCGCACCGGGcctgtgcgc                                                                                                                                                                            >  1:726619/1‑230 (MQ=255)                                                                            gcgcAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGGGAAAAAGCCAAA‑GGTGCGAACGCCCAG‑TTACCACCCCAACCCCGTAGCCATCAACGCCCTACGCCACCGTGCTGTGGCGCACCCATACCCGCCGGCgc                                                                                                                                                          >  1:775452/1‑238 (MQ=255)                                                                                  cGAGGGCTTTTTCCCCCGCACCCCAACCCGAATGGGGGGGGCGTCATTGAAAAAGATGGCGAAATTGGCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTAGGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcg                                                                                                                                                    <  1:253464/239‑1 (MQ=255)                                                                                         tttttACCACGCCTCCCAACCCGAATGTGGGGTGCGTCATTGTCAAAGATGGCGAAATGGTCGGTAAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGG‑GTTGCGTAGGGC‑GGGTAAAAAAGCCAAA‑GGGGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGt                                                                                                                                              <  1:1286785/237‑1 (MQ=255)                                                                                                  cGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGGGGCACCGCCTAT‑GGCACACTctga                                                                                                                                                                                                                     >  1:904652/1‑155 (MQ=255)                                                                                                    cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGGGCGACCGCCCAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCCCTGGCGTAGCGCGCGTGGGTGCCTCGAGGCa                                                                                                                                  >  1:868874/1‑239 (MQ=255)                                                                                                                 aTGGCGGGTGCGTCATTGTCAAAGGTGGCGAAATTGTCGGTAAAGGTTACCACCAACGTGGGGGTGAACCACATGCCGAAGTACCCGC‑GTTGCGTATGGC‑GGGTAAAAAAGCAAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCg                                                                                                                     <  1:449975/239‑1 (MQ=255)                                                                                                                       ggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCCCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGt                                                                                                                >  1:1154930/1‑238 (MQ=255)                                                                                                                           cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGCGGGGGGAAAAACCCAAA‑GG‑GCGGCCCCCCCT‑GTGCCCCTCGAAACCCGTAGCCATCATGGTCGTACGCCCCCCCGCTGTGGCGCCCCCCACCCCGCCGGCCTAGCGCCCGGGGGTGCCCCGCTGCCAAAACCCAACCCACAgggggggg                                                                                                           >  2:449975/1‑233 (MQ=255)                                                                                                                              cATGGTAAAAGGTGGCAAATTTGCAGGTGACGGTTACCCCAAACGTGGGGGTGAACCCCAGGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                        <  1:155321/239‑1 (MQ=255)                                                                                                                                tgtcaaaaGATGGCGAAATTGTGGGGTAAGGTTACCCCCAACGGGGGGGTAAACCACAGGCCGAAGTACACGC‑GTGGCGTATGGG‑GGGTAAAAAAGCAAAA‑GGTGGGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                      <  1:386828/234‑1 (MQ=255)                                                                                                                                          gggCAAAATTGGCGGTGAAGGTTACCCCAAGCGTGGGGGTAAACCCAATGCCGAAGT‑CCCGCGGTTGCGTAGGGC‑GGGTTAAAAAGCAAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                              <  1:400068/236‑1 (MQ=255)                                                                                                                                             cgattttttCGGTAGAGGTTCCCACATACGTGCGGGTAAACCACAGGCCGAGGTACCCGC‑TTGGGGTAGGGG‑GGTAGAAAAGCCCAAG‑GGTCGGACCGCCTAT‑GTCACCCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGt                                                                                          <  2:214080/230‑1 (MQ=255)                                                                                                                                                tttggtgggTGATGGTTCCCACCCGCGGGCGTGCAAACCAGATGCCGACGTACCCGT‑GTTGGGGTTGGT‑GGGTAAAAGACCCAAG‑GGTGCGACCGCCTAT‑GTCCCCCTCGAACCCTGTTGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTg                                                                                       <  2:332363/231‑1 (MQ=255)                                                                                                                                                                       aaCGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTTTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACCCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCTCGGTGCAAGATCCTAACCCCCCGGTCGCCGGGCGTGGACCTTTCCGGCTGCCACAGGCTGGCA‑TTGACGTCAg                                                               >  1:947268/1‑238 (MQ=255)                                                                                                                                                                         cgggCGGGTGACCCATGTGCCGAAGTCCCCGC‑GTTGCGTATGGC‑GGTGAAAAAAGCAAAG‑GGTGCGACCGCCTAT‑GTCACACGCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCa                                                            <  1:1251738/236‑1 (MQ=255)                                                                                                                                                                         cgggCGGGTAAACCCCATGCCGAAGTCCCCGC‑GTGGCGTATGGC‑GGGTGAAAAAGCCAAG‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCa                                                            <  1:34168/236‑1 (MQ=255)                                                                                                                                                                           tGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGGGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGGGCGGGCTTTTCCGCCTGCAACAGGCTGGCA‑TTGACGCCGGCCCCg                                                          >  1:369194/1‑239 (MQ=255)                                                                                                                                                                            gggggAGAACCCCATGCCGAATGAC‑CGT‑GTTGGGGTGGGG‑GGGGGAAAAAACCAAAGGGTGGGGCCGCCTTT‑GGCACACTCGAACCCTGTGGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACg                                                          <  2:486136/236‑1 (MQ=255)                                                                                                                                                                                               aaGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCCGGTCGCTGGGCGGGGACTTTACCGTCTGCCACAGGCTGGCA‑TTGACGTCAGCCCCGGCCCGATGGTGAGGGAAGcc                                      >  1:915124/1‑239 (MQ=255)                                                                                                                                                                                                attccccgc‑gTTGGGTAGGGC‑GGGTGAAAAACCCAAG‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                                     <  1:1295641/233‑1 (MQ=255)                                                                                                                                                                                                    acccgggTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGCCCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                   <  1:277060/232‑1 (MQ=255)                                                                                                                                                                                                           ttGCGTATGGC‑GGGTGGAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGGACCCTGTAGCCATCATGGTCGGACGCCCCCGGGCTGTGACGCACTCCTCGCCGCTGGCGGAGCGCGCCTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCGTCCGCAACAGGCCGGCA‑TTTGCGTCCGCCACCGCCTGATGTTAGGGGAAACCGGGCAATTGaa                           >  2:229860/1‑239 (MQ=255)                                                                                                                                                                                                             gCGTATGGG‑GGGGAAAAACGCAAAG‑GGGGCGACCGCCTAT‑GTCACCCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATa                         <  1:14619/239‑1 (MQ=255)                                                                                                                                                                                                              cGTAGGGC‑GGGTGAAAAAGCAAAG‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGt                                                                                                                                                                                                   <  1:58994/69‑1 (MQ=255)                                                                                                                                                                                                                 gtGGG‑GGGTGAAAAGCCAAAA‑GGTGTGCCCGCCTAT‑GGCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAg                      <  1:404270/237‑1 (MQ=255)                                                                                                                                                                                                                  gggC‑GGGAAAAAAAGCCAAG‑GGGGCGCCCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGc                    <  1:687243/238‑1 (MQ=255)                                                                                                                                                                                                                        ggggAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCCTGCTTTGCAGCCCTCATCGCCCCTGGCGTAGCGCGCCGGGTGGCCTCGATGCAAGATCCTAACCCGCAGGGCGCCGGGCGTGGCCTTTACCGTCTGCAACAGCCTGGCACTTCGCG‑CCGCCCCGGCCCCAAGATGAGTGAAACCGCGAAAATGaaaaaaagggtttc               >  2:687243/4‑231 (MQ=255)                                                                                                                                                                                                                                      a‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGCCTTTCTCAAGCGGATgggc  >  1:534743/1‑238 (MQ=255)                                                                                                                                                                                                                                      a‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGGGAGGCCGACCAATTGAATAAAGGCTTTCTCAAGCGGCTgcgc  >  1:204363/1‑238 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGC‑GTTGCGTATGGC‑GGGTGAAAAAGCCAAA‑GGTGCGACCGCCTAT‑GTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCA‑TTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGC  >  NC_000913/433419‑433883

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