Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,457,392	C→T	E158E (GAG→GAA)	yjgA ←	putative ribosome biogenesis factor YjgA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,457,392	0	C	T	98.4%	222.0 / NA	64	E158E (GAG→GAA)	yjgA	putative ribosome biogenesis factor YjgA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/31);  minor base G (0/1);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCG  >  NC_000913/4457159‑4457634                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   ttCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCg                                                                                                                                                                                                                                   >  1:322826/1‑251 (MQ=255) ctCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCg                                                                                                                                                                                                                                   >  1:322844/2‑251 (MQ=255)   agctcTTAGCCAGATATAGTATGCTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGa                                                                                                                                                                                                                                 <  2:255390/247‑1 (MQ=255)    ggCCTTAGCCAGATATAGTATGCGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTTCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGAt                                                                                                                                                                                                                                <  1:309810/251‑1 (MQ=255)       cTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGa                                                                                                                                                                                                                                 >  1:255390/1‑247 (MQ=255)            ccAGATATAGCATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTATTTTTCTTTCTTCGCGTTGAGGATCAGAGtcc                                                                                                                                                                                                                        >  2:221129/1‑249 (MQ=255)             cAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCAGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACg                                                                                                                                                                                                                       >  2:77411/1‑251 (MQ=255)                  ataGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGtt                                                                                                                                                                                                                  >  1:148222/1‑251 (MQ=255)                     gTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACTCAGTTGCt                                                                                                                                                                                                               <  2:322844/251‑1 (MQ=255)                     gTATGCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACAGCAACAACGAAAATTACCCTGCGGTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCt                                                                                                                                                                                                               <  2:322826/251‑1 (MQ=255)                                 cGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGCTGCGACGATTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAACATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGc                                                                                                                                                                                                                                    >  2:187542/1‑218 (MQ=255)                                 cGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGc                                                                                                                                                                                                                                    <  1:187542/218‑1 (MQ=255)                                 cGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCAt                                                                                                                                                                                                   >  1:344250/1‑251 (MQ=255)                                  ggATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATc                                                                                                                                                                                                  >  2:24497/1‑251 (MQ=255)                                     tAAGGCGTTCACGCCGCATCCGGCATTTGAGCCTAGTGCCTGATGCGACGCTTGATGCGTCTTCTCAGACCTACAGGTGCACCGTATCCGGCAAACCATACTCAAATCAACAACGCAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCggg                                                                                                                                                                                               >  2:47041/1‑251 (MQ=255)                                           gTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGGTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAg                                                                                                                                                                                         <  1:217639/251‑1 (MQ=255)                                            ttCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCcac                                                                                                                                                                                           >  1:339958/1‑248 (MQ=255)                                               aCGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTAACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGAttt                                                                                                                                                                                     <  1:292411/251‑1 (MQ=255)                                                    gcATCCTGCATTTTAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCATAGTACGCAGTTGCTGACGATCGGCAGCGGGCCACAGATTTAATAc                                                                                                                                                                                <  2:277235/251‑1 (MQ=255)                                                     cATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCCAACCATACTCACATCAACAACGCAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATAcc                                                                                                                                                                               >  1:276586/1‑251 (MQ=255)                                                          ggCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGGGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGc                                                                                                                                                                          >  2:315698/1‑251 (MQ=255)                                                               ttGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGTGGCACCGTATCCTGCAAACCATACTCACATCAACAACGAAAATTCCCCTTCGGTCTCCGCTGACTCGCGCAGATACTGGAAAATCTGGCGTGCTGAGGTCGGCGGCGTAGTGCCTGCTTTTACATTCTGCGCGTTGCTGATCAGAGTAGTCAGTTGCTGCCGATCGGCATCGTGccctatgttgtattcgtcagcgatg                                                                                                                                                                      <  2:202163/250‑24 (MQ=255)                                                                  aGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAGATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGt                                                                                                                                                                  >  1:400361/1‑251 (MQ=255)                                                                   gCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCCAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGAAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGtc                                                                                                                                                                 >  2:442409/1‑251 (MQ=255)                                                                    cATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGtcc                                                                                                                                                                >  2:35568/1‑250 (MQ=255)                                                                          gcatGATGCGACGCTTGAGGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCAGCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCCCCt                                                                                                                                                          <  2:424912/248‑1 (MQ=255)                                                                           ccTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTGCCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTGATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCTGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCtt                                                                                                                                                         <  1:169851/251‑1 (MQ=255)                                                                               aTGCGACGCTTGATGCGTCTTATCAGGCCTACCGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAACATTACCCTTCGTTCTCCGCTAACTCGCGCAGATAATGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGCTCGGCATCTGGCCACAGATTTAATACCTCAGCGATGGCGTCATACCCTTgatc                                                                                                                                                     >  2:378596/1‑251 (MQ=255)                                                                               aTGCGACGCTGGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCCTCATCACCTTgatc                                                                                                                                                     <  2:173039/251‑1 (MQ=255)                                                                                tGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTATTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCt                                                                                                                                                          >  1:278002/1‑245 (MQ=255)                                                                                cgCGACGCTTGATGCTTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACAGCACCAACGAAAATTACCCTTCGCTCTCCGCTAACTCGCGCAGATACTGGAAAATCTTGCGCGCGGATTTCGGCGGCTTATGCCCTTCTGTTGCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGAGCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGGCATCCCCt                                                                                                                                                          <  2:278002/244‑1 (MQ=255)                                                                                       cTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTtctttttct                                                                                                                                                                                                                                                 >  2:125919/1‑151 (MQ=255)                                                                                       cTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTtctttttct                                                                                                                                                                                                                                                 <  1:125919/151‑1 (MQ=255)                                                                                         tgatgaGTCTTATCAGGCCTACAGGTGCACCGTATCCGTCAAACCATTCTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATCCTGGAAAATCTGGCGTGCGGATTCCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTAAGCGATGGCGTCATCATCTTGATCGATCAGAc                                                                                                                                             <  1:327613/249‑1 (MQ=255)                                                                                         tGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGAc                                                                                                                                             <  1:315698/249‑1 (MQ=255)                                                                                         tGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTGCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACg                                                                                                                                            <  2:386521/250‑1 (MQ=255)                                                                                             gCGTCTTATCAGGCCGACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGTCGTCATCACCTTGATCGATCAGACGATcg                                                                                                                                        <  2:155935/250‑1 (MQ=255)                                                                                                       aGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGATTACGCAGTTGCTGACGATCTGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTc                                                                                                                             <  1:326108/251‑1 (MQ=255)                                                                                                             aCAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTCCCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTctttttctttct                                                                                                                                                                                                                                             >  1:190230/1‑133 (MQ=255)                                                                                                             aCAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCTCAGATACTGGAAAATCTGGCTTGCGGATGTCGGCGGCTTATTCCCTTctttttctttct                                                                                                                                                                                                                                             <  2:190230/133‑1 (MQ=255)                                                                                                                        gTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGca                                                                                                            <  1:185690/251‑1 (MQ=255)                                                                                                                            ccGGCAAACCATACTCACATCAACAACGAAAATTACCCT‑CGTTCTCCCCCAACTCGCGCAGGTACTGGAAAATCTGGCGGGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTTGCGTTGCGGATCAGAGTAAGCAGTTGCTGACGGTAGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCCTCCCCTTTATCGGACAGACGATCGCGCAGGGTTTCCAGTTTATGGAACAGCaccacc                                                                                                       >  1:70272/1‑251 (MQ=255)                                                                                                                              ggCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACct                                                                                                      <  1:41597/251‑1 (MQ=255)                                                                                                                                aaaaCCATACAAACAACAACAACGAAAATTACCCTGCGTTCTCCGCTAACTCGCTCCGTTACTGGAAAATCTGGCGTGCGGATGTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCTTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCTCctgc                                                                                                    <  2:16607/250‑1 (MQ=255)                                                                                                                                  aaCCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACctgctg                                                                                                  <  1:269879/251‑1 (MQ=255)                                                                                                                                              aTCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGtt                                                                                      >  1:377334/1‑251 (MQ=255)                                                                                                                                                   caacGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCa                                                                                 <  1:35568/251‑1 (MQ=255)                                                                                                                                                        aaaaTTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCCGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACCGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCttt                                                                            >  2:449779/1‑250 (MQ=255)                                                                                                                                                                   tCGTTCTCCGCTACCTCGCGCAGATACTGTAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTGTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGcca                                                                  <  2:138696/250‑2 (MQ=255)                                                                                                                                                                   ccGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGAGTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGAGTTAATACCTCAGCGATGGCGTCAGCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGGGACGGTTTTTCAGCGTTTCCAGCGCCTg                                                                 <  2:65890/250‑1 (MQ=255)                                                                                                                                                                    cGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTg                                                                 >  1:128918/1‑250 (MQ=255)                                                                                                                                                                     gTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGAc                                                               >  2:120384/1‑251 (MQ=255)                                                                                                                                                                         tcCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAAt                                                           >  2:147841/1‑251 (MQ=255)                                                                                                                                                                             cTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGc                                                       >  1:283953/1‑251 (MQ=255)                                                                                                                                                                                   cgcgCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTa                                                   <  1:211574/249‑1 (MQ=255)                                                                                                                                                                                     cgcAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCTATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTc                                               <  1:304794/251‑1 (MQ=255)                                                                                                                                                                                        aGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTAGTCCCTTCTTTGTCTTGCTTCGCGTTGCGGAGCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTGTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTcgcg                                            <  2:398349/251‑1 (MQ=255)                                                                                                                                                                                           tACTGAAAAATCTGGCGTGCGAATTTCGCCGGCTTATTCCCTTATTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGACACTGATTTAATACCTCAGCGATGGCATCATCACCATGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACACCACCACCTGATGGTTGTGACGGATCTTCAGCTAATCCCGCTCCTGTCGACTAGGCCCGACGTCCCGCtc                                         >  2:138991/1‑250 (MQ=255)                                                                                                                                                                                                     tCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATc                               >  1:316494/1‑251 (MQ=255)                                                                                                                                                                                                      cTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCt                              >  2:86787/1‑251 (MQ=255)                                                                                                                                                                                                       tGGCGCGCGGATTTCGCCGGGTTATTCTCATCTTTTTATTTCTTAACGCTGCGGATATGGGAAAGCAGCTGCTGACGATTTGCATAGGGCCACGGATTTAATACCTCAGCGATGGCGTCATCACCTTGCACGATCCGACGATCGCGCAGGTTTTCCAGTTTCTGGAACAGCCCCACCTGCTGGCTGTGCCGGTTTTTCCGCTTGTCCAGCTGCTGACGCATAGGCTCTACGTCGCGCTGGCGCATCATCtt                             >  2:459483/1‑251 (MQ=255)                                                                                                                                                                                                                 aTTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGc                   <  1:24497/251‑1 (MQ=255)                                                                                                                                                                                                                       gcggcTTACTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATCCCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCCCCTGCTGGTTGTGACGGTGTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGccac             >  1:108646/1‑248 (MQ=255)                                                                                                                                                                                                                             tATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggc         <  1:147590/249‑1 (MQ=255)                                                                                                                                                                                                                             gaTTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCATAGTACGCAGTTGCTGACGATCGGCATCGGCCCACATATTTAATTCCGCAGCGATGGCGTCATTGCCTTGATCGATCAGACGATCGCGCAGGGTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAAGAGGCTCTACGTCGCGCTGTCGCAGCATCTTACCGATGAGCTGCAACTggcg        <  2:160296/249‑1 (MQ=255)                                                                                                                                                                                                                              aTTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggcggc      >  1:253914/1‑251 (MQ=255)                                                                                                                                                                                                                                 cccTTCTTTTTCTTTCTTCGCGGTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCTATGAGCTGCAACTggcggcgg    <  1:282291/250‑1 (MQ=255)                                                                                                                                                                                                                                  ccTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGgcggcggcg  >  1:41587/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCG  >  NC_000913/4457159‑4457634

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