Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,979	G→A	V71M (GTG→ATG)	yahH →	putative uncharacterized protein YahH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,979	0	G	A	81.0%	79.3 / 10.2	42	V71M (GTG→ATG)	yahH	putative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (13/21);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAG  >  NC_000913/339730‑340223                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     gTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGt                                                                                                                                                                                                                                                     <  1:260671/251‑1 (MQ=255)    aCGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGAGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACCCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAACGTAGCGCATCAGGAATTTGTCACCATTGCCGGATGCGGCGTGAACGACTTATCCGGC‑CTACGAGTGgc                                                                                                                                                                                                                                                  >  2:442757/1‑251 (MQ=14)               ccTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAAttgt                                                                                                                                                                                                                                       <  1:192464/251‑4 (MQ=255)                            gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAACGTAGCGCCTCAGGCATTTGTCCCCCTTTCCGGATTCGGCGGTGACCCCTTTTCCCGC‑CTCCGTGTGGCGCGGGAAT‑TTTAGGCCTGATaacc                                                                                                                                                                                                                         >  1:236759/1‑249 (MQ=11)                                     ggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCAt                                                                                                                                                                                                                <  2:441825/251‑1 (MQ=9)                                        gtgtaCGCCTTATCCGGCCTTAGAATGGCGCAAGAATCTGTAGGCCCGAAAAGCGTTGCGGATCCGGCATTTGTCACCCTTGGCGGATGCGGCGGGAACGCCCTTA‑CCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                                                                                                             <  2:445827/247‑1 (MQ=11)                                              gCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCAGGATGCGGCGTGAACACA‑TTATACGGCCTACGAATGGCGCAAGCATCTGTAGGCCTGATAAGCGTAGAGCATCAGGCATTTGTCACCATTGCCGGATGCGGAGTGAACGCCTTATACGGC‑CTACGCGGGGCGCGAGAGT‑TGTAAGCCTGATAAGCAGAGCGCATCAGGCaga                                                                                                                                                                                                        >  1:267241/1‑248 (MQ=11)                                                        gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTCCGAGTGGCGGGAGGAT‑TGTAGGCCGGATAAGCGTAccgcctccggcggactggcgttggg                                                                                                                                                                                             >  1:191240/1‑230 (MQ=11)                                                                                                                              aTGCGGCGTGCACGCC‑TTCTCAGTCTTAGGAATGTAGAAAGAATCTGGAGGCGGGATAAGCGGAGCTGGTCAGGCAGTTGTCTTCACTGCCGGATGGGGCGTGAACGCCTGATCCGGC‑CTACGAATGGCGCAAGAATCGGTAGGCCGGATAAGCGTAGCGCAGCAGGCAGTTGGCTCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGAGTGGCGCGATAATTGTAGGTCTGATAAGCGTAgc                                                                                                                        <  2:364572/251‑1 (MQ=2)                                                                                                                                                tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTg                                                                                                         <  2:267701/249‑1 (MQ=2)                                                                                                                                                tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTg                                                                                                         <  1:47332/249‑1 (MQ=2)                                                                                                                                                      gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGGGAATTGTAGGCCTGATAAGCGTAGAGCATCAGGCAGTCTGGCGTTGGt                                                                                                 >  1:46632/1‑251 (MQ=11)                                                                                                                                                       ccTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTc                                                                                                <  1:325921/251‑1 (MQ=255)                                                                                                                                                         tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGGGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCAt                                                                                              >  1:199900/1‑251 (MQ=255)                                                                                                                                                           cGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                                                            <  1:89022/251‑1 (MQ=255)                                                                                                                                                                       gaagcTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCct                                                                                <  2:339191/247‑1 (MQ=255)                                                                                                                                                                         aTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCGTGGCGGGAGAATTGTAGGGCTGATAAGCGGTGCGCATCAGGCAGTCTGGCGTTGGTCATAAccccctcccccccc                                                                              >  1:407299/1‑241 (MQ=11)                                                                                                                                                                             gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTa                                                                           <  1:263158/250‑1 (MQ=255)                                                                                                                                                                                 gCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGa                                                                      <  2:222080/251‑1 (MQ=255)                                                                                                                                                                                           cGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCCGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaac                                                             >  1:254881/1‑250 (MQ=255)                                                                                                                                                                                                 gcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGAGTGGCGCGAGTATTGTTGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGGGTTGGTCATAACCCCATCACCCCCTGTAGCGGACATaacaactttaa                                                        >  1:35059/1‑244 (MQ=14)                                                                                                                                                                                                  cATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAaga                                                     <  1:445370/251‑1 (MQ=255)                                                                                                                                                                                                     cAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTGCGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGGTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCt                                                  <  1:434886/251‑1 (MQ=17)                                                                                                                                                                                                         cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCCCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgtg                                              <  1:73791/251‑1 (MQ=255)                                                                                                                                                                                                         cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGGGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgtg                                              >  1:348467/1‑251 (MQ=255)                                                                                                                                                                                                            ttgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGt                                            <  2:50824/246‑1 (MQ=255)                                                                                                                                                                                                             tgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTc                                           <  1:173061/247‑1 (MQ=255)                                                                                                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACGAGTGGCGCGAGAATTGTAGGCCCGATAAGCGTAGCGCCTCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGAAATAACAACATGAAAGAGGTTGTGGTCGGTGCCATtggg                               >  2:134932/1‑250 (MQ=18)                                                                                                                                                                                                                         cGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCTCCATTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCGACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTGTGGCGTTGGTCATAAGCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTGGTGGTCGTTGCCAAtggtg                              <  2:31426/251‑1 (MQ=18)                                                                                                                                                                                                                         cGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAAAAACATGAAAGAGCTTGTGGTCGTTGCCATttgtg                              >  2:205054/1‑251 (MQ=255)                                                                                                                                                                                                                           gATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGc                            <  2:302776/251‑1 (MQ=255)                                                                                                                                                                                                                                ggcgTGAACTCCTTATCCGGCTCT‑CGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAGGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGGGGCGCGAGAATTGTATGCCTGTTAAGCGTAGCGCAGCAGGCAGTCTGGCGTTGGTCATAGCCCCATCCCCCGCTGGAGCGGACATAACAACATGAAAGAGCGTCTGGGCTTTGCCATTGGTGGCAACAg                       <  2:438238/251‑1 (MQ=17)                                                                                                                                                                                                                                   gtgaaCGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATACGGCCTAAGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCCTCAGGCAGACTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACCTAATACCATGAAAGAGCTTGTGATAGTTGCCATTGGGGGCAAAAGCAt                    >  1:185622/5‑251 (MQ=18)                                                                                                                                                                                                                                    tgaaCGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTCGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCCTCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACCTAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCAtt                   >  2:55779/4‑251 (MQ=255)                                                                                                                                                                                                                                     gaaCGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCAtt                   <  2:97133/248‑1 (MQ=255)                                                                                                                                                                                                                                      aaCGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCGGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATATCCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTAt                 <  1:300036/250‑1 (MQ=255)                                                                                                                                                                                                                                              aTCCGGC‑CTACGAATGGCGCAAGAAGCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGGTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAATAGCTTGTGGTCGTTGCCAGTGGTGGCAACAGCATTATCGAAGAt          <  2:438208/250‑1 (MQ=21)                                                                                                                                                                                                                                                  ggC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAAAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAAACGCATTATCAAAGATAACgc     >  1:25483/1‑251 (MQ=255)                                                                                                                                                                                                                                                    c‑cTACGAATGGCGCAAGCATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCGTGGCGCGCGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGAGGGCATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGAAACAGCATTATCAAAGATAACgccc   >  2:216335/1‑250 (MQ=25)                                                                                                                                                                                                                                                    c‑cTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgcca   <  2:165117/251‑1 (MQ=255)                                                                                                                                                                                                                                                    c‑cTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgcc    >  1:15386/1‑250 (MQ=255)                                                                                                                                                                                                                                                        aCGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGTATTATCAAAGATAACgccag  <  2:255086/249‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     GTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAG  >  NC_000913/339730‑340223

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