Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,143,377	C→T	S283L (TCA→TTA)	frwC →	putative PTS enzyme IIC component FrwC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,143,377	0	C	T	100.0%	160.2 / NA	46	S283L (TCA→TTA)	frwC	putative PTS enzyme IIC component FrwC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/24);  total (22/24)

GCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAAT  >  NC_000913/4143141‑4143606                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                      gCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGc                                                                                                                                                                                                                           >  1:378301/1‑249 (MQ=255)    gATTTGCGTTGCTCAGGGTGTTTATCCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTa                                                                                                                                                                                                                       <  2:264612/250‑1 (MQ=255)                   gggTGTTTATACCGGTGTGGCTAGCGCGTCGGTGGGTGTTTGTATCCCGCCGCTGGGTATGGGGCGGGCGACGCTGATTTGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGGGATGGGCTGTGTTGGCGTTACTGAATGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCGTCCATCATGGTCTGTTTAGTTTGTGGTGCAGTAACTGCTGCGCTGGTCg                                                                                                                                                                                                       <  2:262492/251‑1 (MQ=255)                     gtgtTTATAACGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGt                                                                                                                                                                                                     >  1:132532/1‑251 (MQ=255)                         ttATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTgcgc                                                                                                                                                                                                 >  1:70160/1‑251 (MQ=255)                          tataCCGGGGGGGCCATGGCTGCCGGGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTgcgc                                                                                                                                                                                                 <  2:325796/250‑1 (MQ=255)                                  ggtggCCATCGCTGCGGTGTTTATTGGTATCCCGCCGATGGGTATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGTGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCTCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCTGTGTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTAt                                                                                                                                                                                         <  1:399457/250‑1 (MQ=255)                                            gCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTgggg                                                                                                                                                                              <  1:38003/251‑1 (MQ=255)                                                  gTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCt                                                                                                                                                                        >  1:232038/1‑251 (MQ=255)                                                         tttGTATCCCGCCGCTGGGGATGGGGCTGGCGCCGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCGGAGtgtg                                                                                                                                                                 >  2:335315/1‑251 (MQ=255)                                                             tATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGtgctgc                                                                                                                                                             >  1:341236/1‑251 (MQ=255)                                                                     cgcTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGtt                                                                                                                                                     >  2:299118/1‑251 (MQ=255)                                                                        tGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTgc                                                                                                                                                   >  2:185505/1‑249 (MQ=255)                                                                             aTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCATGCTGGGGTGGTCTGATTGTGCTGGCGGTGGTTGAAGGCaa                                                                                                                                             <  1:148994/251‑1 (MQ=255)                                                                             aTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCaa                                                                                                                                             <  1:99034/251‑1 (MQ=255)                                                                               ggggCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAg                                                                                                                                            <  1:299118/250‑1 (MQ=255)                                                                                 ggCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTg                                                                                                                                         >  1:188070/1‑251 (MQ=255)                                                                                          cGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATc                                                                                                                                >  1:211662/1‑251 (MQ=255)                                                                                             tGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCgc                                                                                                                              >  2:203322/1‑250 (MQ=255)                                                                                                ttGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCgcagc                                                                                                                           >  1:255732/1‑250 (MQ=255)                                                                                                     cGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTgg                                                                                                                      <  1:203322/250‑1 (MQ=255)                                                                                                                     cAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTAGGAGccgggggg                                                                                                     >  2:104598/1‑244 (MQ=255)                                                                                                                       gcagagCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGACGGGGCGATTCCTTTCGCCGCTGCCGCTCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGTCGCTGGTCGGTGCGCAGTGCTATGCCGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAgtgggtac                                                                                                   >  2:329319/3‑247 (MQ=255)                                                                                                                       gaagaGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATGCCTTCCATCAGGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGAc                                                                                                   <  1:188648/251‑1 (MQ=255)                                                                                                                                 aaaCCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGGGGTGACGGCTGTTtgg                                                                                         >  1:346643/1‑250 (MQ=255)                                                                                                                                      ggCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGaa                                                                                    <  1:404092/251‑1 (MQ=255)                                                                                                                                       gcaaggctgCACTGGTGATGGGCTGTGTTGGCGTTACGGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGGGTCATTCCGTCCATCATGGTCGGTTTAGTGTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGTCGGCTGTTTGTGTGAAc                                                                                   <  2:145568/243‑1 (MQ=255)                                                                                                                                             cGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGAGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTg                                                                             <  1:101012/251‑1 (MQ=255)                                                                                                                                              ggCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTg                                                                             >  2:287094/1‑250 (MQ=255)                                                                                                                                                 aCTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGaaa                                                                          >  1:257477/1‑250 (MQ=255)                                                                                                                                                                ggtTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTaaaa                                                          <  2:90940/250‑1 (MQ=255)                                                                                                                                                                   tGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCTGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGGGAACGTGCTGAAATGTCTGGCGCGTAAAAAt                                                        <  2:174677/250‑1 (MQ=255)                                                                                                                                                                               aGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACtg                                           >  1:161104/1‑251 (MQ=255)                                                                                                                                                                                            ttCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAaagaag                               <  2:10055/250‑1 (MQ=255)                                                                                                                                                                                            ttCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAaagaa                                <  2:255456/249‑1 (MQ=255)                                                                                                                                                                                            ttCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGGGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGGGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGAGGGGGCGGGTGTTTGTGGGCACGCGCTGAAAAGGAAGGCGCGTAAAAATGGGTCGTTGACTGATGAAAaaga                                 >  2:400876/1‑248 (MQ=255)                                                                                                                                                                                                    tGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGTCTGATGAAAAAGAAGACGACCtgg                      >  2:125359/1‑251 (MQ=255)                                                                                                                                                                                                     gCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGGAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtg                       <  2:289896/249‑1 (MQ=255)                                                                                                                                                                                                       cGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGATGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggatt                   <  2:124052/251‑1 (MQ=255)                                                                                                                                                                                                       cGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGGTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggatt                   <  2:15301/251‑1 (MQ=255)                                                                                                                                                                                                               ggCGTGGCATTCATTCCATCATGGTCGGTGTAGTTGGTGGTGCAGTAACTGAGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATGGTGCTGCCGGTGGTTGAAGGCAAGATGGGTTATATCGCAGCAGTGGATGTAGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACTACCGGGATTTGGAtttt           <  1:329319/250‑1 (MQ=255)                                                                                                                                                                                                                         tccttccATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAaattaat  <  2:119872/250‑1 (MQ=255)                                                                                                                                                                                                                                  catGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAgt                                                                                                         >  2:11668/1‑138 (MQ=255)                                                                                                                                                                                                                                  catGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAgt                                                                                                         <  1:11668/138‑1 (MQ=255)                                                                                                                                                                                                                                    tGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTgggg                                                                                                                                                                              >  1:66029/1‑67 (MQ=255)                                                                                                                                                                                                                                    tGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTgggg                                                                                                                                                                              <  2:66029/67‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                      GCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAAT  >  NC_000913/4143141‑4143606

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