Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,532	T→C	L119L (TTA→CTA)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,532	0	T	C	100.0%	187.0 / NA	56	L119L (TTA→CTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/32);  total (24/32)

ACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATG  >  NC_000913/4638287‑4638774                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aCCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTaaaa                                                                                                                                                                                                                                               <  1:622637/251‑1 (MQ=255)        ggATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGAtt                                                                                                                                                                                                                                        <  2:232798/250‑1 (MQ=255)          aTGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTc                                                                                                                                                                                                                                      <  1:637969/251‑1 (MQ=255)           tgtgGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCg                                                                                                                                                                                                                                     >  1:406772/1‑251 (MQ=255)              ggAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATg                                                                                                                                                                                                                                  <  1:184833/251‑1 (MQ=255)               gcaGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGGTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTAGCCATTGTTGGTTACCTGAGTCATTGATGGTTGATGTCAATCAGAACGTGGAAGTACGGAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGAGTTAACGCTAAAATCcgatgtcgatgt                                                                                                                                                                                                                                 <  2:333249/249‑1 (MQ=255)                  gaggcgtttCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTc                                                                                                                                                                                                                              <  2:420196/244‑1 (MQ=255)                     gcgattacccgaaGTACCAGCGGGCCGGAAAAACCCGTTGGGCCGCTTATCGCTATTTCTGTGAACGAGCTTTATTCTGTTCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACTCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTcgc                                                                                                                                                                                                                            <  2:103262/239‑1 (MQ=255)                             ctgacgTACCATCGGGCAGCAAACACTCGTGGTGCCGCTTATCGCGATTACTGTGACCTAGCTTTATACGGTGTAGGTAGAGGATAAAACCGTGGATCGGACACGAAGGTTTATCCATTTTTGGGTACCTGAGTCATTGATTGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGAGTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATGTCGATGTCGCGCGTCtga                                                                                                                                                                                                                      <  2:577287/243‑3 (MQ=255)                             ccaaacttacatcaggCCGCAAAAACTCGTTGGGCCGCGTATCGCTATTCCTGTTCCCGAGCTTTATACGGTGCAGTAATAGGATAAATCCGTTGAGCGGTAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCTATCAGAACGTGGAAGAACGCAAGTTATGGATTTATTCCGGTCAGGTCTGGCACAGTGATTTAACGCGAACAGCCGTTTTCGATGTTTCGCGTCTTAGCg                                                                                                                                                                                                                   <  2:352120/237‑1 (MQ=255)                                 aGTACCAGCGGGCAGCCAAAAATAGGTGGGCCGCTTATCGATATTACTGTGCCCGAGATTTATACGGTGCAGGAAGAGGATAAAACCGTGGCGCGGAAACGAAGATTAATCCATTATTGGTTACCTGAGTCATTGATGGGTGATGGCAATCAGAACGTGGAAGACAGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTGAACGCTAAAAGCAGAATTCGATGTTTCGCGTCTTAGCGagct                                                                                                                                                                                                               >  2:386817/1‑248 (MQ=255)                                  gTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTc                                                                                                                                                                                                              >  1:485521/1‑251 (MQ=255)                                  gTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTc                                                                                                                                                                                                              <  1:487096/251‑1 (MQ=255)                                    aCCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCaa                                                                                                                                                                                                            <  1:59453/251‑1 (MQ=255)                                          gggCCGCAAAAACTCGTTAGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGCGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTAGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGCTTTAACGCTAAAAGCCAATTTCGATGTTTAGCGTCTTAGCGAGCTCACCGCGcc                                                                                                                                                                                                      >  2:380662/1‑251 (MQ=255)                                           ggCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGt                                                                                                                                                                                                                                 <  2:85220/223‑1 (MQ=255)                                           ggCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGt                                                                                                                                                                                                                                 >  1:85220/1‑223 (MQ=255)                                           ggCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCa                                                                                                                                                                                                     <  1:29485/251‑1 (MQ=255)                                                     acgcGTTGGGCCGCTTATCGCTATTCCTGTGTCCGAGCTTTAGTCGGTTCAGGAAGAGGATAAAACTGTGGAGCGGAAACGAAGTTTTATACTTTTTTGGTTACCTGAGTCATTGATTGTTGATGGCAATCAGAACGTGTAAGAACGCAATATAGGTATTTATAACTGTCAGGTCTGGACCAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAtatagcagct                                                                                                                                                                                           <  2:29651/248‑7 (MQ=255)                                                         tttGGGCCGCTTATCGCTATTCCTGTGACAGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGTGTTTTATACCGGTCAGGTCTGGCACAGGGATTTTACGCTAAAAGCCGCTTTCGATGTTTCGCGTCTTAGCGAGCACAACGCTACAAATATAATCTTAgg                                                                                                                                                                                       >  2:132353/2‑251 (MQ=255)                                                                 gcttaTCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGccct                                                                                                                                                                               <  1:262900/247‑3 (MQ=255)                                                                   ttaTCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAAcg                                                                                                                                                                                                          >  2:477850/3‑222 (MQ=255)                                                                   ttaTCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAAcg                                                                                                                                                                                                          <  1:477850/220‑1 (MQ=255)                                                                    ttttgctcttcttgtgACCGAGATTTGTACGTTGCAGGAAGAGGATAAATCCGTGGATCGGAAACGAAGGTTTATCAATTTTTGGTTACCTGAGTCATTGATGGTTGATTGCAATCAGAACGTGGTAGAACTCAAGATAGGGATTTATACCGGTCAGGTCTGGCACTGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGTGCTCAACGCGCCAAAGATCATCTTTTGCAAGctcttt                                                                                                                                                                             <  2:146983/238‑6 (MQ=255)                                                                                  tgACAGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGAAGCGGAAACGCTGTTTTATCCAGTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGAGAGGGATTTAAACCGGTCAGGTCTGGCACAGTGCTTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAATGCACCACATATCATCTTAGGAACGCACTTTATTGTGATTAGCGTc                                                                                                                                                              >  2:357736/1‑251 (MQ=255)                                                                                     ccGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCggg                                                                                                                                                           >  2:577543/1‑251 (MQ=255)                                                                                     ccGAGCTTTATACGGTGAAGGTAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTGACCTGAGGCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGTGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCATTTATTGTGATTAGCGTCggg                                                                                                                                                           <  2:294568/251‑1 (MQ=255)                                                                                       gAGCTTTATTCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAAATAAGGTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATTCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGa                                                                                                                                                         <  2:485521/251‑1 (MQ=255)                                                                                                  ggTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCggcg                                                                                                                                                             <  2:82922/236‑4 (MQ=255)                                                                                                  ggTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCg                                                                                                                                                             >  1:82922/1‑236 (MQ=255)                                                                                                      cAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATtgg                                                                                                                                          >  2:282900/1‑251 (MQ=255)                                                                                                                  aaaaCCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAgc                                                                                                                              >  1:18843/1‑251 (MQ=255)                                                                                                                                 aaaCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGagc                                                                                                                                                                                                                <  2:508824/154‑3 (MQ=255)                                                                                                                                 aaaCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGagc                                                                                                                                                                                                                >  1:508824/1‑152 (MQ=255)                                                                                                                                   aCGAAGGTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaacggaa                                                                                                             <  2:201852/251‑1 (MQ=255)                                                                                                                                     gAAGTTTTATCCATTTTTGGTTAACTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGCTAGGGATTTATACCGGTCAGGTCTGGCAAAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGACAAATATTATCTTAGGCAAGCCCTATATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGAGCCTGAAGTTaacggcacg                                                                                                           >  1:162306/1‑251 (MQ=255)                                                                                                                                      atgTTTTATCTATTTGTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCTCTTTATTGTGATGTGCTTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaatggaacgg                                                                                                          <  2:102153/249‑1 (MQ=255)                                                                                                                                                       ggTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAccc                                                                                         >  1:153704/1‑251 (MQ=255)                                                                                                                                                       ggTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAccc                                                                                         <  1:563015/251‑1 (MQ=255)                                                                                                                                                       ggTTACCTGAGGCATTGATGGTTGATGGCAAGCAGAACGTGGAAGAACGCAATAGAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGTTGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGGTAACGGAACGGCGCTGATCATTGAAcc                                                                                          <  2:164429/250‑1 (MQ=255)                                                                                                                                                           aCCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGc                                                                                      <  1:631852/250‑1 (MQ=255)                                                                                                                                                                 gTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGt                                                                               >  1:404862/1‑251 (MQ=255)                                                                                                                                                                                 ggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGGACCCGGCACCGGGTTAGAGCAAGGCGGGc                                                                >  2:112640/1‑250 (MQ=255)                                                                                                                                                                                     aTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGc                                                           <  1:282900/251‑1 (MQ=255)                                                                                                                                                                                        aGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTg                                                        <  1:600655/251‑1 (MQ=255)                                                                                                                                                                                                                  aTTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCTTGGTGTGGTGAAAGCGCCTGAAGTTACAGGAACGGCGCTGACCATTGAACCCGGCCACGGGGTAGAGCAAGGCGGGCAAGGCGTGCATAACCCTTTACCCGCAGGGGGCTg                              >  2:379258/1‑251 (MQ=255)                                                                                                                                                                                                                     tataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGc                            <  1:180003/250‑1 (MQ=255)                                                                                                                                                                                                                        aCCGGTCAGGTCTGGCACAGTTATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGAGAGCTCAACGCGCAAAATATCATCTTAGGCAAGCCCTTTATTGAGATTAGCGTCGGGGAGGCGCGGGGTATTGGTGGGGTGAAAGCGCCTGAAGTCAACGTAAAGGCGCTGACCATTGAACACGGCCCCGGGTTAGAGCAAGGCGGGGCGGGGGGGCATATCCCTTCACCCGAAaggggcatgagggc                        >  2:4996/1‑241 (MQ=255)                                                                                                                                                                                                                           ggTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGa                                                                                                                        <  2:291215/152‑1 (MQ=255)                                                                                                                                                                                                                           ggTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGa                                                                                                                        >  1:291215/1‑152 (MQ=255)                                                                                                                                                                                                                           ggTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGGAGc                      >  1:14395/1‑250 (MQ=255)                                                                                                                                                                                                                               aGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAc                 <  1:380662/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGATAACGGAACGGCGCTGACCATTGAACCCGGCCCCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGGCTGGCGGAAGCAAAACCTg              >  2:219591/1‑248 (MQ=255)                                                                                                                                                                                                                                          aGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCgg                                                                                       >  1:429680/1‑170 (MQ=255)                                                                                                                                                                                                                                          aGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCgg                                                                                       <  2:429680/170‑1 (MQ=255)                                                                                                                                                                                                                                              aTTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCGGGGCGTGCATATCCCTTTACCCGCAGGGGACTGGCgggagaggaaacttaaccgtaatatg  >  2:477872/1‑227 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   ACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATG  >  NC_000913/4638287‑4638774

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