Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	882,452	T→G	T95P (ACC→CCC)	deoR ←	DNA‑binding transcriptional repressor DeoR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	882,452	0	T	G	97.6%	136.8 / NA	42	T95P (ACC→CCC)	deoR	DNA‑binding transcriptional repressor DeoR
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (21/20);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCA  >  NC_000913/882205‑882697                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        gccacccccgccGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGGTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTCCTCTTACAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTGCATTAGCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTg                                                                                                                                                                                                                                                      <  1:407307/245‑1 (MQ=255)   cccccccgccGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGAt                                                                                                                                                                                                                                                    <  2:390474/247‑1 (MQ=255)                       aaGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCCGCGTCGc                                                                                                                                                                                                                                 >  1:676048/1‑250 (MQ=255)                           gATATTCGGTCAAAAATTATTCAGCGTTTGCTGACAAGCGATGGGTTTGAAAATGGCGTTGCGGGCGTGAATTTCACCACCGCTAAGATACGAGCGGCAATGGGGTTTCTCTGTCAGCTCCAGAAAGTTATTTAGCGAAGAACAAACGGCGGTATAAGGGATTTCATTATCAAACGCTGCAATAATCCGCGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGcc                                                                                                                                                                                                                             <  2:541785/250‑1 (MQ=255)                                                           gAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTccacca                                                                                                                                                                                              <  1:75360/249‑1 (MQ=255)                                                            aaaaTCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGc                                                                                                                                                                                           <  1:94435/251‑1 (MQ=255)                                                              aaTCGATGGGTTTGAAAATGGCGTTGCTGGCGTGCAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGAAAAAGGGGTTTCATTATCAATCGGTTCAATAATCCCCGGCGTGGTGGGGCCACCATCAAAAAAGAGGGGCTGATCAGGTTTTACCAGCGTCGCCGCCAGTTTTTCCGCCCGGGGTTTTTTTTCCACCAggagg                                                                                                                                                                                         >  2:295986/1‑248 (MQ=255)                                                                 cGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAATTTTGCCGCCCGGCGTTTTTCTTCCCCCAAGCggggt                                                                                                                                                                                      >  2:187352/1‑249 (MQ=255)                                                                                 ggCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAAc                                                                                                                                                                       <  1:274419/250‑1 (MQ=255)                                                                                    gTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGt                                                                                                                                                                   <  1:222685/251‑1 (MQ=255)                                                                                       gcgggCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAAGCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAt                                                                                                                                                                <  2:281806/248‑1 (MQ=255)                                                                                          ggCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGtt                                                                                                                                                                                                                                              <  2:512041/170‑1 (MQ=255)                                                                                          ggCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGtt                                                                                                                                                                                                                                              >  1:512041/1‑170 (MQ=255)                                                                                             gTGAAATTCACCGCCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCCGGCGGGATTTTTGATCGCTTAACAGGGAAtggctgg                                                                                                                                                          >  1:76417/1‑251 (MQ=255)                                                                                                       ccaccGCAAAGAAACGCGCGGCCATGGGGTTTCTCTTTAAGCGCCAGAAAGGAATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCAACAGTCAAAAAAGAGGGGCTAATGGGGGTATACCAGCGCCGCAGACAGTTTTGCCGCCCGGCGTTTTTCTTCCAACAGGCGGGAGTTTTGCTCGCTATACAGGAAAtggatggccctcgccgg                                                                                                                                                >  2:105156/1‑242 (MQ=255)                                                                                                          ccGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTAGTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTc                                                                                                                                            <  2:76417/251‑1 (MQ=255)                                                                                                           cGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTcc                                                                                                                                           <  1:486631/251‑1 (MQ=255)                                                                                                           cGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATGATCAATCGCTTCAATTATCCACGGCGTGGTGGTGCCACAGGCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTGTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGCAATGGCTGGCACTGCGCGG‑TTcc                                                                                                                                           <  2:77289/251‑1 (MQ=255)                                                                                                                         cgGCAATGGGGTTACTCTTTCAGCTCCAGTAAGGTAGTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAAGCCACGTCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGTAGCCAGTTTTGTCGCACGGCGTTTTTCTTCTACCATGCGGTATTTTTGTTCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAgcc                                                                                                                             <  1:310218/251‑1 (MQ=255)                                                                                                                                ggggTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGGAATGGCTGGGAATGCGCGG‑TTCCAGAACAATATAGCCGCCGAGc                                                                                                                      >  1:566407/1‑251 (MQ=255)                                                                                                                                 gggTTTCTCTTTCAGCGCCAGAAAGGAATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCCCCGTCAAAAAAAAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGATTTGCCGGCCGGCGTTTTTCTTCCACCAGGCGGGGTTTTTGATCGCTTAAAAAGTAAATGCTTGCACTGCGCGGGTTCC‑CAACAATATAAccgccgggg                                                                                                                      >  2:176697/1‑247 (MQ=255)                                                                                                                                        tctTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGGGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCACGGCGTTTTTCTTCCCCCAGGCGGGCTTTTTGATCGCTTAACAGGTAATGGTTGGCCCTTCGCGG‑GTCCAGAAAAATATATCCGACGACCAAAacagc                                                                                                              >  2:200969/1‑248 (MQ=255)                                                                                                                                               gcgcCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCa                                                                                                       <  1:475203/251‑1 (MQ=255)                                                                                                                                                cgcCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCCAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGGTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACCGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAACCGCCGAGCAAAACGACGGGCGCAc                                                                                                      >  1:348069/1‑251 (MQ=255)                                                                                                                                                          gTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCCTAACAGGTAATGGCTGGCACTGCGCGG‑GTTCAGAACAATATAGACGCCGAGCAAAACGACGGGCGCACTGTGg‑ttgtt                                                                                           >  1:210416/1‑251 (MQ=255)                                                                                                                                                              ttAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCCGGCGGGATTTTTGATCGCTTAACCGGGAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCCCCGGGCAAAacgacg                                                                                                             >  1:660462/1‑230 (MQ=255)                                                                                                                                                              ttAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCGTCAATAATCCACGGCGTGGTGGTGACACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTGTGAGCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAacgacg                                                                                                             <  2:660462/230‑1 (MQ=255)                                                                                                                                                                  cGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCCGGCGGGGTTTTTGTTCGCTTAACAGGTAATGGCTGGCACTGCCCGG‑TTCCCGAAAAATATATCCGCCGAGCAAAAAGACGGGCGCAC‑GTGGGTTGTTCAGATccc                                                                                   >  2:44608/1‑249 (MQ=255)                                                                                                                                                                       aaCAAACGGCGGAAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGGCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGGTTTTTGATCGCTTAAAAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGagacg                                                                               >  1:387036/1‑250 (MQ=255)                                                                                                                                                                           aaCGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACtgt                                                                                                   <  1:133267/227‑1 (MQ=255)                                                                                                                                                                                     aaaGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGAt                                                                                                                                                                               >  1:343225/1‑142 (MQ=255)                                                                                                                                                                                     aaaGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGGGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGAt                                                                                                                                                                               <  2:343225/142‑1 (MQ=255)                                                                                                                                                                                           aTTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑GTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATAGCGACGAATCGTCATCTCCGCAAcccc                                                          >  1:88582/1‑251 (MQ=255)                                                                                                                                                                                               cATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGACGAGCAAAACGACCGGCGCACTGGGG‑TTGTTCAGAACGCAACGAATCGTAATCTCCGAAACACCa                                                         >  2:273098/1‑248 (MQ=255)                                                                                                                                                                                                aTTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCa                                                     <  1:187352/251‑1 (MQ=255)                                                                                                                                                                                                      aaTCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGggcgg                                               >  1:18066/1‑251 (MQ=255)                                                                                                                                                                                                                 taatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCttt                                     >  1:33668/1‑250 (MQ=255)                                                                                                                                                                                                                            gtggtggtgCCAAAGTCAAAAAAGAGGGGCTGTTCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACtt                         <  2:210416/251‑1 (MQ=255)                                                                                                                                                                                                                                            cAAAAAAGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGAACGGCGCACTGTGG‑TTGTTCAGATCGCGACGAACCGTCATCTCCGAAACCCAAAGCAGGGCGGCGGCGTCTTTACGATGTAACTTTTCGCTCCGTTTTAAt         >  2:466219/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaaGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGc                 <  2:376641/239‑1 (MQ=255)                                                                                                                                                                                                                                                aaaGAGGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGc                 >  1:376641/1‑239 (MQ=255)                                                                                                                                                                                                                                                   gagGGGCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTgca  <  1:279765/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        TGCACGCCCGCCGCAGAATAAAAAGCGATATCCGGGCAAAAATTATTCAGCGTTTGCTGAAAATCGATGGGTTTGAAAATGGCGTTGCTGGCGTGAAATTCACCACCGCAAAGAAACGCGCGGCAATGGGGTTTCTCTTTCAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGG‑TTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGG‑TTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCA  >  NC_000913/882205‑882697

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