Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,742	T→C	H380H (CAT→CAC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,742	0	T	C	97.5%	108.6 / NA	40	H380H (CAT→CAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/16);  minor base A (1/0);  total (24/16)
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

CACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCC  >  NC_000913/4620499‑4620982                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 cACCAACTTCGCTGACGTCGACTCTTCCTGGGGCCACCGACGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGaa                                                                                                                                                                                                                                            <  1:206977/250‑1 (MQ=255) cACCAACTTAGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAAACGGGCTCACTGGGTCACCGTGaa                                                                                                                                                                                                                                            >  1:89039/1‑250 (MQ=255)   ccAACTTCGTTGACTTCGACGCGTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGaaa                                                                                                                                                                                                                                           <  2:20637/249‑1 (MQ=255)      acctCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGGCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTGACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTc                                                                                                                                                                                                                                      <  1:138504/248‑1 (MQ=255)      aCTTCGTTGAAGTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTc                                                                                                                                                                                                                                      <  2:94385/251‑1 (MQ=255)                ccgcgacTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGGTATGCTGCGGGTCTGGAGCTGATCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATGCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGAGGCACCGGGAAACCTGCTCGggtgatg                                                                                                                                                                                                                            >  2:246063‑M2/4‑244 (MQ=255)                     aCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCAGCGGGTCTGGAGCGGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTACGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGACACCTTCGCGggtgatgctgcc                                                                                                                                                                                                                       <  2:59558‑M2/251‑13 (MQ=255)                     aCTCTTCCTGGGGCCACAGTCGCGACGTCGCCGGCTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCAGGTACTGGTTTACGGCCCGAAAGTAAAAACGGGCTCACTGGGTCACCGTGAAACCTCCGCGggtgatgaggcc                                                                                                                                                                                                                       <  2:185004‑M2/251‑13 (MQ=255)                           ccTGGGGCCCCCGTCGCGACGTCGCCGGTTATTCTGCGGGTCTGGAGCTGTTCGACCGCAGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCCCGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGGTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgcaaactt                                                                                                                                                                                                                  >  2:180835‑M2/1‑233 (MQ=255)                           ccTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaactt                                                                                                                                                                                                                  >  1:261859‑M2/1‑233 (MQ=255)                                  ccACCGTCGCGACGTCGCCGGTTATGCTGCGGGGCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgattt                                                                                                                                                                                                          >  2:194916‑M2/1‑226 (MQ=255)                                                   ccGGATATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgttt                                                                                                                                                                                         <  1:244858‑M2/251‑43 (MQ=255)                                                    cGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttt                                                                                                                                                                                        <  1:115088‑M2/251‑44 (MQ=255)                                                      gTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttacggatttagtgtatgatggtgtttttg                                                                                                                                                                                      >  2:251146‑M2/1‑206 (MQ=255)                                                         aTGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagg                                                                                                                                                                                   >  1:5093‑M2/1‑203 (MQ=255)                                                           gctgcgggtCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTAATGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAAAACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtg                                                                                                                                                                                 <  1:255320‑M2/248‑51 (MQ=255)                                                              gcgGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgct                                                                                                                                                                               >  1:122444‑M2/1‑198 (MQ=255)                                                                           tGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggagatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgt                                                                                                                                                                 <  1:266329‑M2/251‑67 (MQ=255)                                                                               cGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGACAGTAAAACCGGGCTCACTGGGTCACCGTGAACCCTTCGCGggtgatgctgccaacttgctgatttcgtgtatgatggtgtttttgaggtgctcccgtggcttctgtttct                                                                                                                                                             <  1:64533‑M2/251‑71 (MQ=255)                                                                                        tCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgt                                                                                                                                                    <  1:4862‑M2/251‑80 (MQ=255)                                                                                                       gTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagc                                                                                                                                      >  1:268417‑M2/1‑157 (MQ=255)                                                                                                                gcgcGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggg                                                                                                                            >  1:81605‑M2/1‑148 (MQ=255)                                                                                                                       gacgacATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACCCGCGTGAACACCTTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggagatgctgccaacttactgatttagtgtatgatggtgtttttgag                                                                                                                                                                                    <  2:37648‑M2/188‑48 (MQ=255)                                                                                                                       gacgacATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACGGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgag                                                                                                                                                                                    >  1:37648‑M2/1‑141 (MQ=255)                                                                                                                        acgacAGCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctccggttcagctactgacggggtggtgcg                                                                                                                    >  1:97809‑M2/1‑140 (MQ=255)                                                                                                                         cgacATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgt                                                                                                                   >  1:250661‑M2/1‑139 (MQ=255)                                                                                                                                            gCTGACCACGGTTGTGATCCGGCCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcgtctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgcc                                                                                                 >  1:162745‑M2/1‑120 (MQ=255)                                                                                                                                            gCTGACCACGGTTGTGAGCCGACCTGGACCGGTACTGACCACACGCGTGAACACATCCCGGTACGGGTTTACGGCCCGAAAGAAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgaggctgccaacttactgatttattgtatgatggtttttttgaggtgctccagtggcttttgtttcaattagctgtccctcctgttcagctactgacgggggggtgcgtaacggcaaaagcaccgccg                                                                                                >  2:125134‑M2/1‑120 (MQ=255)                                                                                                                                               gACCACGGTTGTGATCCGACCTGGACCGGTCCTGACCACACGCGTGAACACATTCCGGTACTGGTTAACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTTAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctccggttcagctactgacggggtggtgcgtaacggcaaaagcaccgccgga                                                                                              >  2:226929‑M2/1‑117 (MQ=255)                                                                                                                                                            aTCCGACCTGGACCGGTACGGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatct                                                                                >  1:37201‑M2/1‑104 (MQ=255)                                                                                                                                                                gACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaccttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggggcgtaacggcaaaagcaccgccggacatcagcgctatctctgc                                                                            <  2:40405‑M2/251‑152 (MQ=255)                                                                                                                                                                     ggACCGGTACTGACCACACGCGTGAACCCATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacgtggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctca                                                                       >  2:101656‑M2/1‑95 (MQ=255)                                                                                                                                                                     ggACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctc                                                                        <  1:42641‑M2/250‑156 (MQ=255)                                                                                                                                                                                           tGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACATTCGCGggtgatgctgccaacttactggtttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctaatgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactg                                                 <  1:241091‑M2/251‑179 (MQ=255)                                                                                                                                                                                            gAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgtgcagctacggacggggtggtgggtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgc                                                >  1:41970‑M2/1‑72 (MQ=255)                                                                                                                                                                                                    tCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcac                                         >  1:169248‑M2/1‑64 (MQ=255)                                                                                                                                                                                                             ggtttacggcccgaaaGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacaccagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccg                                 >  2:168883‑M2/8‑55 (MQ=255)                                                                                                                                                                                                                ttacggcccgAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgagggtgttttggaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgcgggacatcagcgctatctcggctctcactgccggaaaacatggcaactgcagttcacttacaccgcttc                             >  2:262986‑M2/5‑52 (MQ=255)                                                                                                                                                                                                                           aaGTAAAACCGGGCTCACTGGGTCTACGTGAAACCTTCGCGggtgaggctgccaacttactgacttagtgtatgctggtgtttttgaggcgctccagtggcttctgtttctatcagctctccttcctgttcagctactgtcggggtggtgcgtaacggcaaaatcaccgccggacatcagcgttatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtac                 >  2:64533‑M2/1‑41 (MQ=255)                                                                                                                                                                                                                                          cACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcat  >  1:263221‑M2/1‑26 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 CACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCC  >  NC_000913/4620499‑4620982

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