Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,869,637	C→A	L87L (CTC→CTA)	cbrA →	colicin M resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,869,637	0	C	A	100.0%	143.0 / NA	40	L87L (CTC→CTA)	cbrA	colicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/26);  total (14/26)

GTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTA  >  NC_000913/3869392‑3869873                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             gTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGAACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCCCCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGGCGTCGCCGCATCGCTAAcgcg                                                                                                                                                                                                                                         >  1:394890/1‑251 (MQ=255)      gATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCTCATCGCTAACGCGTAACt                                                                                                                                                                                                                                    <  1:111920/251‑1 (MQ=255)       attattGGCCTCGGCACGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGCAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTAGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCCGATTTTCCGCGTCAAAACCGACGACGTCGCCGCATCGCAAACGCGTAACt                                                                                                                                                                                                                                    >  2:261678/1‑250 (MQ=255)           ttGGCCTCGGCCCGGCAGGATCGGCGTTGACACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCCCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGACACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGCTGTGATCGCTCATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACtccc                                                                                                                                                                                                                                 >  2:444404/1‑246 (MQ=255)                   ggCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCt                                                                                                                                                                                                                        >  2:59846/1‑250 (MQ=255)                    gCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCtat                                                                                                                                                                                                                      <  1:59846/251‑1 (MQ=255)                    gCCCGGCAGGATCGGCGTTGGCACGAAAGGTAGCAGGCAAAAGGCAGGTGAGCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGGTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACTCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCATCGAAGCtat                                                                                                                                                                                                                      <  2:65433/251‑1 (MQ=255)                      ccGGCAGGATCGGCGTGGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACGGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAACCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCTCATCGCTAACGCGTAACTTCCAGCGAAGCTatat                                                                                                                                                                                                                    <  2:394890/251‑1 (MQ=255)                       cGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGCCGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATc                                                                                                                                                                                                                   >  1:53110/1‑251 (MQ=255)                             gATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATAtt                                                                                                                                                                                                             <  2:59172/251‑1 (MQ=255)                                    gTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGcc                                                                                                                                                                                                      >  2:311269/1‑251 (MQ=255)                                       ggCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTGTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTGCCAGCGAAGCTATATCAATATTAATCGCCATg                                                                                                                                                                                                   <  2:86720/251‑1 (MQ=255)                                          aCGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCtt                                                                                                                                                                                                <  1:356693/251‑1 (MQ=255)                                                gTTAGCAGGCAAAATGCAGGTGATCGCTGTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGAGCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACGACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACt                                                                                                                                                                                          <  2:53110/251‑1 (MQ=255)                                                           aaaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGaaa                                                                                                                                                                               <  2:289355/251‑1 (MQ=255)                                                               gCAGTTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCTTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCg                                                                                                                                                                            <  2:263000/250‑1 (MQ=255)                                                                cAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCt                                                                                                                                                                          <  1:108227/251‑1 (MQ=255)                                                                                 tAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTCACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGCTGAAATCGCTGATCCCCGCCAGCGTTg                                                                                                                                                         >  1:413640/1‑251 (MQ=255)                                                                                 aaaaaaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCGGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTGGTGGATGAAATCGCTGATCCCCGCCAGCGTTg                                                                                                                                                         <  2:379991/250‑1 (MQ=255)                                                                                           cAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCa                                                                                                                                               <  1:409744/251‑1 (MQ=255)                                                                                                    aCTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCt                                                                                                                                      >  2:277666/1‑251 (MQ=255)                                                                                                       gAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTTTTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCtgtg                                                                                                                                   <  2:209615/251‑1 (MQ=255)                                                                                                                   aaaCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTg                                                                                                                       >  1:285156/1‑251 (MQ=255)                                                                                                                           gggcgggttgcgggCACCGGACGCGCAGCGATCTTTTATCCGCGATGGGCTGACGCTTCCTGGCTATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAgg                                                                                                               <  2:150242/239‑1 (MQ=255)                                                                                                                           gggcggTCTGCTGGCACCGGACGCGCAGCTATCTTTTATTCGCGATTGGCGGACGATTCCTGTCGCGGTGATCGCTAATCCGCAGATTTTCAGCGGCAAAACCGTCGACGTCGCCGCATCCTTAACGCGTAATTCACAGCGAAGTTATAGCAATAGTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTGGAGTTTTATCACGACAGCCTGGGCCGAAAAATCTGGCGTGAgg                                                                                                               <  2:288156/250‑1 (MQ=255)                                                                                                                             gcggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgat                                                                                                             <  1:238923/251‑1 (MQ=255)                                                                                                                               ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatg                                                                                                            >  2:62194/1‑250 (MQ=255)                                                                                                                                gTCTGCTGGCACCGGACGCGCAGCGATCTTTTATCCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATAGTAATCGCCATGCTTCCGACTTGTGGAGGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatgat                                                                                                          <  2:202001/251‑1 (MQ=255)                                                                                                                                       ggCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGc                                                                                                   <  2:445194/251‑1 (MQ=255)                                                                                                                                            cGGACGCGCAGCGATCTTTTATGCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTc                                                                                              <  1:337143/251‑1 (MQ=255)                                                                                                                                                    cAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCg                                                                                      <  1:62194/251‑1 (MQ=255)                                                                                                                                                             tttATCCGCGAGGGGCTGACGTTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGAc                                                                               <  1:348449/249‑1 (MQ=255)                                                                                                                                                                aTTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTg                                                                          <  1:166629/251‑1 (MQ=255)                                                                                                                                                                     cgATGGCCTGACGCATCCTGTCGATGTGATCGATAATCCGCAGATTATCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGATATATAACTATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTAATCACGAACGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAGTGGCATGTCATTTTTCGTGCAGACGGATGGGAGc                                                                     >  2:427965/1‑251 (MQ=255)                                                                                                                                                                       aTGGGCTCACGCTTCCTGTCGCTGTGAGCGCTAATCCGCAGATTTCCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACTCGTAACTACCAGCGAAGCGAGACCAATATTAATCCCAATTTTTTCGACTTGTTGAGGAAATCGCTGATCCCCGCCTGCTTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCGGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGc                                                                     <  1:286938/249‑1 (MQ=255)                                                                                                                                                                                     ccTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgcccg                                                     >  1:398049/1‑251 (MQ=255)                                                                                                                                                                                                   cGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTAGTCGGTg                                       >  2:313822/1‑251 (MQ=255)                                                                                                                                                                                                        aTCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGAt                                  >  1:334754/1‑251 (MQ=255)                                                                                                                                                                                                                   ttCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACt                        <  1:403362/250‑1 (MQ=255)                                                                                                                                                                                                                                        gTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTa  <  1:92956/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             GTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTA  >  NC_000913/3869392‑3869873

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