Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,290,120	T→A	11.3%	Y281F (TAT→TTT)	acs ←	acetate‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,290,120	0	T	A	11.3%	210.0 / 3.0	75	Y281F (TAT→TTT)	acs	acetate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (59/7);  new base A (2/7);  total (61/14)
Fisher's exact test for biased strand distribution p-value = 5.15e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGTTTTTCAGCGCGTCA  >  NC_007779/4289888‑4290353                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          gggTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGGGCC‑CCAGGAGATATCACCCGGATGGTAATCAAAGACATATTTAAA‑GGTCAGCCCCGCGTACACCAAAAAAccg                                                                                                                                                                                                                                      >  1:2116632/1‑238 (MQ=255)      cGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGgta                                                                                                                                                                                                                                >  1:367320/1‑239 (MQ=255)                   atTGACCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCC‑CCAGTAGATATCCCCCGGATGAGAATCAAAAACCTATTTTAA‑GGGCAGCTCCCCGGACACCAGATAACCCCCGGTCGTAAgcagaccc                                                                                                                                                                                                                    >  2:1241997/1‑234 (MQ=255)                   atTGACCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGcaca                                                                                                                                                                                                                    >  1:551200/1‑238 (MQ=255)                   atTGACCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGGGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCCGCGCCGCGTCCACCAGATAACCGCCGGTAGTATGCAGcaca                                                                                                                                                                                                                    >  1:1507115/1‑238 (MQ=255)                       aCCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGGAGTATGCAGCACACCttt                                                                                                                                                                                                               >  1:1910879/1‑239 (MQ=255)                            aTGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCC‑CCAGTAGATATCCACCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGCACACCAGATAACCGCCGGGAGTATGCAGCACACCTTTTGGCt                                                                                                                                                                                                          >  1:1323361/1‑239 (MQ=255)                            aTGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGAAATATTTAAA‑GGTCAGCGCCGCGAACACCAGATAACCGCCGGTAGTATGCAGCAAACCTTTGGGCt                                                                                                                                                                                                          >  1:930282/1‑239 (MQ=255)                            aTGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGGAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCt                                                                                                                                                                                                          >  1:1477854/1‑239 (MQ=255)                             tGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTgg                                                                                                                                                                                                            >  1:1228549/1‑236 (MQ=255)                                   cccccccccagcgcCAGAC‑‑GGCAGGCGGCCGGCCTGGTGGGTCCTCCTTCAAACGTCGGGGTGG‑CCCGCCGCAGGGCCGGGGGGCCATACAGCAAGTAGTTGTGCCCGGTCC‑CCCACCCCACGTCGG‑GGGGCAGAGAGTATACTCCCCCGGGAAGATAAAAAAAGACATATTAAAA‑GTCCCGCGCCGCGAACCCCAGATATCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGt                                                                                                                                                                                                   <  2:1670139/229‑1 (MQ=255)                                        caccTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGaa                                                                                                                                                                                               >  1:849112/1‑238 (MQ=255)                                         ccccccgcCACACGGGGCGGGC‑GGCGGCCGTTTGGGTACCCTTAAAAACATCGGGGGGGCCGCACCGCG‑GCCCAGGGGCCCGTACAGCAAGTACCTGTGTCCGGCCC‑CCCACCCCACACCGGCGGTGCC‑CCGGTAGATCCCCCCCGGATGATAATCAAAGA‑ATTTTTAAAGGGTCAGCGCCGCTTACACAAGATACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAAc                                                                                                                                                                                              <  2:997846/234‑1 (MQ=255)                                         accTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAAc                                                                                                                                                                                              >  1:1491098/1‑238 (MQ=255)                                                  aTAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGGAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCCCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCCTTTGGCTTACCGGGAGAACCGGAGGGGTa                                                                                                                                                                                    >  1:736750/1‑239 (MQ=255)                                                   tAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGGAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGGGTag                                                                                                                                                                                   >  1:604076/1‑239 (MQ=255)                                                           ggcc‑GTCGGCCGTTTGGCTCCCCCTTAAACCGTCCGCGGGCTCCCCCGGCA‑GGCAGGCGGCCCGTACAGCATTTATGTGTGCCCGCCCA‑CCCACCCCACCTCGGCTGGGCC‑CCAGTCGATCTCCCCCGGATGAAAATAAAAGACATATTTAAG‑GGCCGGGTCCGCGTACACCATATACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATa                                                                                                                                                                             <  2:1495484/235‑1 (MQ=255)                                                            ggC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGTGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCCGCGCCGCGTAAACCAAAAAACCGCCGGAAGTATTCAGCACCACTTTTTGCTTACCGGGAGAACCGGAGGGGGAGAGaaaaaa                                                                                                                                                                           >  2:430708/1‑239 (MQ=255)                                                              c‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGGAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGCGGAGAGAATAAACa                                                                                                                                                                         >  1:1430300/1‑239 (MQ=255)                                                                gTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCC‑CCCGGAGAAAACACCCGGGATATAATCAAAGACATAATTAAA‑GGGCAGCGCCCCGTACCACAGATAAACCCCGGGAGTAGTCAGCACCACTTTTGGGTTACCCGTAGAACCGCGGGGGGAGAGAAAAAAccg                                                                                                                                                                        >  2:1218326/1‑237 (MQ=255)                                                                              tACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGt                                                                                                                                                          >  1:611888/1‑239 (MQ=255)                                                                                       aaaCATCGGCGTGGTCGCCCCGCA‑GGCCGGGGGGCCGTACAGCAAGTAACGTTGCCCGGTCC‑CCCACCCCACATCGGGGGTGCC‑CCAGTAGATACCCCCCGGATGATAATAAAAGACATATTAAAA‑GGTCAGCGCCGCGTACACCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCtt                                                                                                                                                 <  1:430708/239‑1 (MQ=255)                                                                                             gtgCGGGGCCGCCCCGCG‑GGCCGGCGGGCGGAACAGAAAG‑AACCGGGTCCCGGCACCCCACCCCACGACGGCGGTGAA‑CAGGGAAACACCCCCGGGATAATAAAAAAGAAAATATTAAAG‑GGCCGCCGCCCCTTACACAAGAAACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGcc                                                                                                                                            <  2:1883471/236‑1 (MQ=255)                                                                                             cAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACCTCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGCAGTATGCAGCACACCTTTTTGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGAGCTTCGGCg                                                                                                                                                           >  1:2057561/1‑223 (MQ=255)                                                                                             cAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATgcagca                                                                                                                                                                                                                      >  1:1755634/1‑164 (MQ=255)                                                                                                     ccGCACCGCA‑GGCCAGCGCGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGGACC‑CAAGGAGATATACCCCGCCTGAACATCCACCACAAACTTAAC‑AGTCCGCCCCCAGGACACCACATAACCGCCGGGAGTATGCAGCCCACCCTTTGGCCTACCGGTAGAACCGGAGGTGTGGCGAATACACAGCGGAGATTCGCCGTTCATCACTTCGGCGTAGTGcgtg                                                                                                                                   >  2:867447/2‑236 (MQ=255)                                                                                                         aCCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGGAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGGGTAGAGAATAAACAGCGGATCTTCTGCGTTCATCTCTTCCGCCCGGGGCTggt                                                                                                                                  >  1:1883471/1‑234 (MQ=255)                                                                                                             cA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCCTCGGCGTTCATCCCTTCCGCCTGGTGCTGAtcgctc                                                                                                                             >  1:231335/1‑237 (MQ=255)                                                                                                                      cGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGAGGATCATCCAAGCCATATTTAAA‑GGCCAACGCCGCGTACACCAGATAACCCCCGGCAGTAAGCAGCACACCTTTTGGCGTACCGGTAGAACCCGAAGGGGAGAGAAAAAACAGCAGACGTTCGGCGTTTCTCCCCTCCGGCTGGTGGTGATTGCTCGGCTGCTCa                                                                                                                    >  2:382001/1‑238 (MQ=255)                                                                                                                          ccGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAg                                                                                                               >  1:1749790/1‑239 (MQ=255)                                                                                                                                aGCAATAAACGGTGCCCGGCCA‑CCCGGCCACCTTCGGCGGTGCA‑CCAGTAAAAACCCCCGGGATAATAATAAAAGAAATTTTAAAG‑GGTCGGCGCGGCGTACACCAAAAACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTg                                                                                                         <  1:1489824/239‑1 (MQ=255)                                                                                                                                  aaattAACTGTGTCCGGTCC‑CCCAGCCCACATCGGGGGTGCA‑CCAGTAGATATCCCCGGGATGATAATAAAAGACATATTAAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGt                                                                                                                                                                                                   <  1:2155213/143‑1 (MQ=255)                                                                                                                                   aattaccTGTGCCCGGCCC‑CCCACCCCACATCGGCGGTCCA‑CCAGTAGATACCCCCCGGATAATAATAAAAGACATATTAAAA‑GGCCGGCGCCGCGTACACCAAAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTg                                                                                                         <  1:1241997/230‑1 (MQ=255)                                                                                                                                    aGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGcca                                                                                                      >  1:1651577/1‑238 (MQ=255)                                                                                                                                       aaCTGTGCCCGCCCA‑CCCAGCCCACATCGGCGGTGCC‑CCAGTAGATACCACCCGGATGAAAATAAAAGAAATATTTAAA‑GGCCAGCCCCGCGTACACCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCcaccac                                                                                                  <  1:457629/239‑1 (MQ=255)                                                                                                                                            tgtCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATaaa                                                                                                                                                                           >  1:564395/1‑161 (MQ=255)                                                                                                                                            tgtCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGGGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGt                                                                                              >  1:1462832/1‑238 (MQ=255)                                                                                                                                            tgtCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGGATGCAGCACACCTTTTGGCTTACCGGTAGAACCCGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCCCAggc                                                                                              >  1:2111644/1‑237 (MQ=255)                                                                                                                                            tgtCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGGAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGGGCTGATCGCTCGCTTGCTCAACCAGGGCGTGCCACCACAGGt                                                                                              >  1:1030692/1‑238 (MQ=255)                                                                                                                                            tgtCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAAATATAACCCGGATGATAATCAAAGACATATTTAAA‑AGGCAGCGCCGCGTACACCAGATAAACCCCGGTAGACTGCAGCACACCTTTTGGCTTACCGGTAGAACCgggg                                                                                                                                                                                         >  2:2155213/1‑145 (MQ=255)                                                                                                                                                           gCCCACACCGGCGTGCCA‑CACGAAGATACCCCCCGAATAAAAATCAAAGATATTTAAAAC‑GCCCAGCGCCGCATACACAAGAAACCCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGcc                                                                               <  2:604076/239‑1 (MQ=255)                                                                                                                                                             ccACATGGGCGGTGCA‑CGAGTGAATATCACCGGGAAAAAAATAAAAAACATATTAAAG‑GGCCCGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCAACTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGc                                                                                <  1:1218326/236‑1 (MQ=255)                                                                                                                                                             ccACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCCCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCCGCCAg                                                                             >  1:1769349/1‑239 (MQ=255)                                                                                                                                                             ccACACCGTCGTTGCA‑CAGGTATATACCCCCCGGATGATAATCAAAGCAATATTTAAA‑GGTCGGCCCGGCGTACACCGAATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGc                                                                                <  1:916553/236‑1 (MQ=255)                                                                                                                                                                caTCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCa                                                                          >  1:1686789/1‑239 (MQ=255)                                                                                                                                                                 aTCGGGGGTGCA‑CCAGTAGATACCACCCGGATAAAAATAAAAGAAATATTAAAA‑GGTCGGCGCCGCGTACACCAGAAAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGAtcgctcg                                                                                                                            <  1:24485/188‑1 (MQ=255)                                                                                                                                                                      cggTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGAACGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAAtttt                                                                     >  1:660376/1‑238 (MQ=255)                                                                                                                                                                         tGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACAACAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGGGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTccc                                                                  >  1:93834/1‑238 (MQ=255)                                                                                                                                                                                 gTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTAc                                                          >  1:2060334/1‑239 (MQ=255)                                                                                                                                                                                    gATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCt                                                       >  1:1514335/1‑239 (MQ=255)                                                                                                                                                                                       atCACCCGGATGATAATCAAAGACATATTTCAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTAGTGCCACCACAGGTCGCGCCCTTCGTGCCAGTCAATTTTCCCGCCAGTACGCTTCa                                                    >  1:1732274/1‑239 (MQ=255)                                                                                                                                                                                       atCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTCCGCTTCa                                                    >  1:2167815/1‑239 (MQ=255)                                                                                                                                                                                       atCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACATTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCAt                                                                                                                                                      >  1:777470/1‑141 (MQ=255)                                                                                                                                                                                         cACCCGGATGATAATCAAAGCCATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACAGCCGGTCGCATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTAGTGCCTCCACAGGTCGCGCCCTTCCTGCCAGTCCATTTTCCCGCCAGTACGCTTCAGt                                                  >  1:1182054/1‑239 (MQ=255)                                                                                                                                                                                         cACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTc                                                                                    >  1:133704/1‑205 (MQ=255)                                                                                                                                                                                         cACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGt                                                  >  1:1750491/1‑239 (MQ=255)                                                                                                                                                                                         cACCCGGATGATAATCAAAGAAATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCAAAGGTCGCGCCCTTCCTGCCCGTCAATTTTCACGCCAGTACGCTTCCGt                                                  >  2:1234011/1‑239 (MQ=255)                                                                                                                                                                                          aCCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGAGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACATTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCAAGTCAATTTTCCCGCCAGTACGCTTCAGt                                                  >  1:554840/1‑238 (MQ=255)                                                                                                                                                                                            ccGGATGATAATCCAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTAGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAAACAGGTCGTTCCACCACAGGTCGCCCCCTTCCTGCCAGTCAATTTTCCCGCCCGTCACCTTCAGTac                                                >  2:255049/1‑238 (MQ=255)                                                                                                                                                                                            ccGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTa                                                 >  1:622040/1‑237 (MQ=255)                                                                                                                                                                                                           aGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTAc                                 >  1:886539/1‑238 (MQ=255)                                                                                                                                                                                                               atatTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACg                                >  1:1680049/1‑235 (MQ=255)                                                                                                                                                                                                                 atTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCCCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACAACAGGTCCCGCCCTTCCTGCCAGTCCATTTTCCCGCCAGAACGCCTCAGTACCACCACATACTCTACGCTGGt                           >  2:2115434/1‑238 (MQ=255)                                                                                                                                                                                                                   ttAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGCTCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCCCCCTTCCTGCCAGCCAATTTTCCAGCCAGTACGCTTACGTACCATCAAATGCTATACGCCTGTGAc                        >  2:1256711/1‑239 (MQ=255)                                                                                                                                                                                                                          gTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGc                                                                                >  1:864704/1‑177 (MQ=255)                                                                                                                                                                                                                            cAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGt                >  1:1941776/1‑239 (MQ=255)                                                                                                                                                                                                                            cAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGt                >  1:166034/1‑239 (MQ=255)                                                                                                                                                                                                                             aGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTg                                                                                                         <  1:219353/149‑1 (MQ=255)                                                                                                                                                                                                                             aGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTg                                                                                                         >  2:219353/1‑149 (MQ=255)                                                                                                                                                                                                                             aGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTGAACGGTAGAACGGGAGGTGTAGCGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCGCAACCAGGGCGCGCCGCCACAGCTCGCGCCCGTCCGGCCAGTCAATTTTCCCGCCTGACCGCTTCAGTCCTACCACATGCGCTACGCTGGTGCTGTGCGGGtt               >  2:1488298/1‑239 (MQ=255)                                                                                                                                                                                                                                    cgTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGTTTTTCAgc        >  2:1011541/1‑239 (MQ=255)                                                                                                                                                                                                                                    cgTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGTTTTTCAgc        >  1:1786681/1‑239 (MQ=255)                                                                                                                                                                                                                                    cgTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGACGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGTTTTTCAgc        >  1:1049374/1‑239 (MQ=255)                                                                                                                                                                                                                                         acCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGTTTTTCAGCGCGtc   >  1:250822/1‑239 (MQ=255)                                                                                                                                                                                                                                           cAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGTTTTTCAGCGCGtca  >  1:653478/1‑238 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          GGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATAC‑GGGCAGGC‑GTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCA‑GGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCA‑CCCAGCCCACATCGGCGGTGCA‑CCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAA‑GGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTCTTCCGCCTGGTGCTGATCGCTCGCTTGCTCAACCAGGTCGTGCCACCACAGGTCGCGCCCTTCCTGCCAGTCAATTTTCCCGCCAGTACGCTTCAGTACCACCACATGCTCTACGCTGGTGACGTTCGGGTTTTTCAGCGCGTCA  >  NC_007779/4289888‑4290353

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