Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,577,196	T→G	M133L (ATG→CTG)	gntK ←	gluconate kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,577,196	0	T	G	98.6%	238.8 / ‑3.2	71	M133L (ATG→CTG)	gntK	gluconate kinase 2
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (37/33);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

AAGCCAGGAAAGCGTGCATACGCGCCTTCATGACTAAAAACAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACATTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGT  >  NC_000913/3576968‑3577428                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         aaGCCAGGAAAGCGTGCATACGCGCCTTCATGACTAAAAACAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCCCCAACCGTTGGGtt                                                                                                                                                                                                                                   >  2:85359/1‑236 (MQ=255)      aGGAAAGCGTGCATACGCGCCTTCATGACTAAAAACAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTaaa                                                                                                                                                                                                                              >  1:635525/1‑236 (MQ=255)      aGGAAAGCGTGCATACGCGCCTTCATGACTAAAAACAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTaaa                                                                                                                                                                                                                              >  1:345196/1‑236 (MQ=255)       ggAAAGCGTGCATACGCGCCTTCATGACTAAAAACAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCG‑‑ATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGaa                                                                                                                                                                                                                           >  1:21450/1‑236 (MQ=255)                                gACTAAAAACAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggc                                                                                                                                                                                                    <  1:446373/236‑1 (MQ=255)                                    aaaaaCAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCt                                                                                                                                                                                                >  2:26679/1‑236 (MQ=255)                                            cagcagTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTTCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAAt                                                                                                                                                                                          <  1:380629/234‑1 (MQ=255)                                            cagcagTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCGGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATca                                                                                                                                                                                        >  2:404007/1‑236 (MQ=255)                                            cagcagTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAAt                                                                                                                                                                                          >  2:380629/1‑234 (MQ=255)                                                agTAAAACTGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACCACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAAt                                                                                                                                                                                          <  2:135533/230‑1 (MQ=255)                                                          aCCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCt                                                                                                                                                                          >  1:698693/1‑236 (MQ=255)                                                            ccTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCttt                                                                                                                                                                        >  1:326452/1‑236 (MQ=255)                                                                      ttAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACGTCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGCAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATg                                                                                                                                                              >  2:633004/1‑236 (MQ=255)                                                                          aaCAAGCGTTAATGTAGTCACTTCTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAag                                                                                                                                                          <  1:287257/236‑1 (MQ=255)                                                                               gCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACCCCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAgagaga                                                                                                                                                     <  1:179582/236‑1 (MQ=255)                                                                                cGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGAt                                                                                                                                                    >  2:411250/1‑236 (MQ=255)                                                                                cGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGAt                                                                                                                                                    >  1:556668/1‑236 (MQ=255)                                                                                   tAATGTAGTCACTACTTATTTGCCTTTTTGAATAACCTCAATGGTGCTTGCCGCAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTGTCGTCCTCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAgagaga                                                                                                                                                     <  2:553623/232‑1 (MQ=255)                                                                                      tgtatgtacttctTATTTGCCGTTTTGAATAACCCAAATGGAGATTGCCACAACACCTTCCAGCGGTTGATCGATACCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGAt                                                                                                                                              <  2:309303/229‑1 (MQ=255)                                                                                               tactTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAATTAGATGAAAGGGCGATTCGGATTACATTCAc                                                                                                                                     >  2:200862/1‑236 (MQ=255)                                                                                                actTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCCAATAGATGACAGAGAGATTCGGATTACCTTCAc                                                                                                                                     >  2:210896/1‑235 (MQ=255)                                                                                                 ctTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAACTCAATGAAAGAGAGATTCGGCTTAACttgaggt                                                                                                                                   >  2:406830/1‑231 (MQ=255)                                                                                                    aTTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACGCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCACATAGATGAAAGAGAGATTCGGATTACCGTCAagcagc                                                                                                                                >  2:347486/1‑236 (MQ=255)                                                                                                     tttGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCaaa                                                                                                                                                                    <  2:370821/199‑1 (MQ=255)                                                                                                     tttGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCaaa                                                                                                                                                                    >  1:370821/1‑199 (MQ=255)                                                                                                     tttGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACgcagca                                                                                                                               <  2:326452/236‑1 (MQ=255)                                                                                                        gCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACATTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGt                                                                                                                            <  1:540817/236‑1 (MQ=255)                                                                                                         ccTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTc                                                                                                                           >  1:42558/1‑236 (MQ=255)                                                                                                         ccTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTc                                                                                                                           >  1:72803/1‑236 (MQ=255)                                                                                                              tttAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGAt                                                                                                                      <  1:26679/236‑1 (MQ=255)                                                                                                                   taaCCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAgtgt                                                                                                                 <  2:39895/236‑1 (MQ=255)                                                                                                                             ggTGCTTGCCACAGCTCCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCTCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATg                                                                                                       <  1:443344/236‑1 (MQ=255)                                                                                                                                gCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAg                                                                                                    <  1:105595/236‑1 (MQ=255)                                                                                                                                 cTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCaa                                                                                                                                                                     >  2:399788/1‑170 (MQ=255)                                                                                                                                 cTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCaa                                                                                                                                                                     <  1:399788/170‑1 (MQ=255)                                                                                                                                 cTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATg                                                                                                       >  1:404486/1‑232 (MQ=255)                                                                                                                                 cTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATg                                                                                                       <  2:404486/232‑1 (MQ=255)                                                                                                                                 cTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGa                                                                                                   >  1:510166/1‑236 (MQ=255)                                                                                                                                  ttGCCACAACACCTTCCAGCGGTTGATCG‑‑ATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGaa                                                                                                  <  2:21450/234‑1 (MQ=255)                                                                                                                                    gCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACa                                                                                                >  2:41323/1‑236 (MQ=255)                                                                                                                                         aacaCCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGa                                                                                           >  1:223029/1‑236 (MQ=255)                                                                                                                                           cacCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATc                                                                                         >  1:508837/1‑236 (MQ=255)                                                                                                                                                  cAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGaca                                                                                  <  1:102014/236‑1 (MQ=255)                                                                                                                                                   aGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGacac                                                                                 <  2:510166/236‑1 (MQ=255)                                                                                                                                                         tGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTttat                                                                            <  1:404007/235‑1 (MQ=255)                                                                                                                                                          gATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTttatt                                                                           >  2:650873/1‑235 (MQ=255)                                                                                                                                                             ttgccatcaccaccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTg                                                                       <  1:234196/229‑1 (MQ=255)                                                                                                                                                               atatCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTgcg                                                                     <  1:411250/236‑1 (MQ=255)                                                                                                                                                                tatCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTgcg                                                                     <  1:650873/235‑1 (MQ=255)                                                                                                                                                                tatCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTgcg                                                                     <  2:345196/235‑1 (MQ=255)                                                                                                                                                                tatCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTgcg                                                                     <  2:223029/235‑1 (MQ=255)                                                                                                                                                                  tCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGaca                                                                                  <  1:420246/220‑1 (MQ=255)                                                                                                                                                                        accaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGc                                                            <  2:556668/236‑1 (MQ=255)                                                                                                                                                                         ccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCATCAATTCGCGATAGAGATTTTTCAATGCAGAACAGACGATCAGCGACACTttatt                                                                           >  2:420246/1‑220 (MQ=255)                                                                                                                                                                             tACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAAc                                                        >  1:50199/1‑235 (MQ=255)                                                                                                                                                                               cagcGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAAcgcg                                                     <  1:200862/233‑1 (MQ=255)                                                                                                                                                                                aTCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACCCTTTATTAGTGCGCTGCATAGCAAAcgcg                                                     >  1:88103/1‑235 (MQ=255)                                                                                                                                                                                 tCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAACGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGCCGATCAGCGACACTTTATTAGTGAGCTGCATAGCAAACgcggc                                                   >  1:336625/1‑236 (MQ=255)                                                                                                                                                                                    gTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAGTGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGCCGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGcgt                                                 >  2:466742/1‑235 (MQ=255)                                                                                                                                                                                       tCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAgagaga                                                                                                                                                     >  1:724585/1‑132 (MQ=255)                                                                                                                                                                                       tCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAgagaga                                                                                                                                                     <  2:724585/132‑1 (MQ=255)                                                                                                                                                                                       tCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGcgtcgt                                              >  1:736645/1‑235 (MQ=255)                                                                                                                                                                                       tCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGAGTTTTTTCAATGCAGAACAGACGATTAGCGACACTTTATTAGTGCGCTGCATAGCAAACACGCagtcgt                                              >  2:610163/1‑235 (MQ=255)                                                                                                                                                                                         gTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACGTTATTAGTGCGCTGCATAGCAAAGCCGGCGTCGTTc                                            >  2:637013/1‑235 (MQ=255)                                                                                                                                                                                               aCCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGcc                                      <  2:736645/235‑1 (MQ=255)                                                                                                                                                                                                      tCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCa                              <  2:635525/236‑1 (MQ=255)                                                                                                                                                                                                               gTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCg                     >  1:110282/1‑236 (MQ=255)                                                                                                                                                                                                                 ttCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGAt                   <  1:347486/236‑1 (MQ=255)                                                                                                                                                                                                                          cGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCa                              <  2:363836/216‑1 (MQ=255)                                                                                                                                                                                                                          cGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCa                              >  1:363836/1‑216 (MQ=255)                                                                                                                                                                                                                                  aCAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGt  <  2:72803/236‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         AAGCCAGGAAAGCGTGCATACGCGCCTTCATGACTAAAAACAGCAGCAGTAAAACAGACCCTACTGCTGTTAAAACAAGCGTTAATGTAGTCACTACTTATTTGCCTTTTTTAATAACCTCAATGGTGCTTGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACATTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGT  >  NC_000913/3576968‑3577428

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