Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,986	G→A	R73Q (CGA→CAA)	yahH →	putative uncharacterized protein YahH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,986	0	G	A	85.4%	88.2 / 5.5	41	R73Q (CGA→CAA)	yahH	putative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (14/21);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTC  >  NC_000913/339744‑340229                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    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)                                                                                                                                                                                                                                              atagcgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACCCCTTATCCGGCCTACGAATCGCGCGAGAATTGTAGGCCCGACAAGGGTAGAGCATCAGGCAGCTTGGCGTTGGTCATAACCCCATCACCCTCTGGAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATAATCAAAGATAACGccagccagac  >  1:178980/4‑249 (MQ=21)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    CCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTC  >  NC_000913/339744‑340229

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