Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,941	C→T	T58I (ACT→ATT)	yahH →	putative uncharacterized protein YahH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,941	0	C	T	95.0%	90.3 / ‑6.3	40	T58I (ACT→ATT)	yahH	putative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (19/19);  total (20/20)
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

TCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTC  >  NC_000913/339700‑340182                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  tCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAAGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATg                                                                                                                                                                                                                                          >  1:397292/1‑251 (MQ=255)  cGAGAAAGCCATTCTTGGATGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑CTATCCGGCCTACAAAAGGCGCCCGAATCTGGAGGCCTGATAAACGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATgc                                                                                                                                                                                                                                         >  2:366130/1‑251 (MQ=18)   gagaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGCATGGCGCAAGAATCTGTAGGCCTGATAAGGGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATgcg                                                                                                                                                                                                                                        >  2:77089/1‑251 (MQ=255)     gaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCCAGGAACTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCCCCATTGCAGGATgcggc                                                                                                                                                                                                                                      >  1:87526/1‑251 (MQ=18)     gaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATgcggc                                                                                                                                                                                                                                      <  1:93039/251‑1 (MQ=255)       aaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGt                                                                                                                                                                                                                                    <  1:391917/251‑1 (MQ=255)         gCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTg                                                                                                                                                                                                                                   >  1:24609/1‑250 (MQ=255)          ccATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTgaa                                                                                                                                                                                                                                 <  1:308378/251‑3 (MQ=255)              tCTTGGCTGCTGCGACAATTCCGGCGACTGACATACATGCTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGACAGGCGCCCGAATATGTAGGCCTGATAAGCGTGGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTTGCCCCC‑TTATTCCGCCTCCGAATGGCGCACGAATCTGCAGGCCGGCTGAACGTAGCGCAT‑CAGGCATTTGTGCCCATTGACGGATGCAGCGTGGACGcc                                                                                                                                                                                                                             >  2:37508/1‑251 (MQ=2)                     tGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTGTCACCATTGCCGGTTGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCCAGAATCCTTAGGCCTGATAAGCGTAGCGCCT‑CCGGCATTTGTCACCACTGCCGGATGCGGCGTGGaacccctttccc                                                                                                                                                                                                                      >  2:47332/1‑240 (MQ=17)                            aaCGGTACGGCGTCTGGCCTACATGTTCATGCCGGATGCGGCGTGATCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATGTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CATGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACg                                                                                                                                                                                                               <  1:282909/251‑1 (MQ=14)                             aCGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACg                                                                                                                                                                                                               >  1:97133/1‑250 (MQ=255)                               gTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGt                                                                                                                                                                                                            <  1:260671/251‑1 (MQ=255)                                  aCGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGAGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACCCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAACGTAGCGCAT‑CAGGAATTTGTCACCATTGCCGGATGCGGCGTGAACGACTTATCCGGCCTACGAGTGgc                                                                                                                                                                                                         >  2:442757/1‑251 (MQ=14)                                             ccTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAttgt                                                                                                                                                                                              <  1:192464/251‑4 (MQ=255)                                               tACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGACTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGACAGCC‑TTATCCGTCCTACGAATGGCGCAAGAATTTGTAGGCCTGATAAGCGGAGCGCATTCAG‑CATTTGTCACCATTGACTGATGCGGCGGGGACCCCTTCTCCggtccacgcgtggggggagcatgtt                                                                                                                                                                                              >  1:235567/1‑226 (MQ=14)                                                          gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAACGTAGCGCCT‑CAGGCATTTGTCCCCCTTTCCGGATTCGGCGGTGACCCCTTTTCCCGCCTCCGTGTGGCGCGGGAAT‑TTTAGGCCTGATaacc                                                                                                                                                                                >  1:236759/1‑249 (MQ=11)                                                                   ggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCAt                                                                                                                                                                       <  2:441825/251‑1 (MQ=9)                                                                      gtgtaCGCCTTATCCGGCCTTAGAATGGCGCAAGAATCTGTAGGCCCGAAAAGCGTTGCGGATCCGGCATTTGTCACCCTTGGCGGATGCGGCGGGAACGCCCTTA‑CCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                                                                    <  2:445827/247‑1 (MQ=11)                                                                            gCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCAGGATGCGGCGTGAACACA‑TTATACGGCCTACGAATGGCGCAAGCATCTGTAGGCCTGATAAGCGTAGAGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGAGTGAACGCCTTATACGGCCTACGCGGGGCGCGAGAGT‑TGTAAGCCTGATAAGCAGAGCGCATCAGGCaga                                                                                                                                                               >  1:267241/1‑248 (MQ=11)                                                                                      gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTCCGAGTGGCGGGAGGAT‑TGTAGGCCGGATAAGCGTAccgcctccggcggactggcgttggg                                                                                                                                                    >  1:191240/1‑230 (MQ=11)                                                                                                                   cTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCg                                                                                                                                                                                                                                              >  2:84582/1‑133 (MQ=255)                                                                                                                   cTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCg                                                                                                                                                                                                                                              <  1:84582/133‑1 (MQ=255)                                                                                                                                                            aTGCGGCGTGCACGCC‑TTCTCAGTCTTAGGAATGTAGAAAGAATCTGGAGGCGGGATAAGCGGAGCTGGT‑CAGGCAGTTGTCTTCACTGCCGGATGGGGCGTGAACGCCTGATCCGGCCTACGAATGGCGCAAGAATCGGTAGGCCGGATAAGCGTAGCGCAGCAGGCAGTTGGCTCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGAGTGGCGCGATAATTGTAGGTCTGATAAGCGTAgc                                                                               <  2:364572/251‑1 (MQ=2)                                                                                                                                                                              tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTg                                                                <  1:47332/249‑1 (MQ=2)                                                                                                                                                                              tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTg                                                                <  2:267701/249‑1 (MQ=2)                                                                                                                                                                                    gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGGGAATTGTAGGCCTGATAAGCGTAGAGCATCAGGCAGTCTGGCGTTGGt                                                        >  1:46632/1‑251 (MQ=11)                                                                                                                                                                                     ccTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTc                                                       <  1:325921/251‑1 (MQ=255)                                                                                                                                                                                       tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGGGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCAt                                                     >  1:199900/1‑251 (MQ=255)                                                                                                                                                                                         cGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                   <  1:89022/251‑1 (MQ=255)                                                                                                                                                                                                     gaagcTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCct                                       <  2:339191/247‑1 (MQ=255)                                                                                                                                                                                                       aTCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCGTGGCGGGAGAATTGTAGGGCTGATAAGCGGTGCGCATCAGGCAGTCTGGCGTTGGTCATAAccccctcccccccc                                     >  1:407299/1‑241 (MQ=11)                                                                                                                                                                                                           gTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTa                                  <  1:263158/250‑1 (MQ=255)                                                                                                                                                                                                               gCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGa                             <  2:222080/251‑1 (MQ=255)                                                                                                                                                                                                                         cGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCCGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaac                    >  1:254881/1‑250 (MQ=255)                                                                                                                                                                                                                               gcAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGAGTGGCGCGAGTATTGTTGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGGGTTGGTCATAACCCCATCACCCCCTGTAGCGGACATaacaactttaa               >  1:35059/1‑244 (MQ=14)                                                                                                                                                                                                                                cAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAaga            <  1:445370/251‑1 (MQ=255)                                                                                                                                                                                                                                    cAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGGTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCt         <  1:434886/251‑1 (MQ=17)                                                                                                                                                                                                                                        cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCCCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgtg     <  1:73791/251‑1 (MQ=255)                                                                                                                                                                                                                                        cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGGGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgtg     >  1:348467/1‑251 (MQ=255)                                                                                                                                                                                                                                           ttgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGt   <  2:50824/246‑1 (MQ=255)                                                                                                                                                                                                                                            tgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTc  <  1:173061/247‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  TCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTC  >  NC_000913/339700‑340182

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: ‑