Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,143,280	G→A	A251T (GCG→ACG)	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,280	0	G	A	100.0%	131.4 / NA	39	A251T (GCG→ACG)	frwC	putative PTS enzyme IIC component FrwC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/21);  total (18/21)

GGGATGCAGCAGGGCAGCATTGTTATGCTGGCGGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGG  >  NC_000913/4143040‑4143511                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        gggATGCAGCAGGGCAGCATTGTTATGCTGGCGGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGt                                                                                                                                                                                                                               <  1:227229/251‑1 (MQ=255)               cagcaTTGTTATGCTGGCGGTGATACTGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCCGGGTGTTTCTCCCGGGGTGGCCATCGCTGCCGTGGGTCTTTGTATCCCGCCGCTGGGGATGGGGCTGGCGCCGCTGATTGGTCGTAAAAATTTCTCCGCTGGAGCGCGCGAACCCGGCAAAACGGCACTGGTGagtgggtgggttg                                                                                                                                                                                                                 >  2:244722/1‑239 (MQ=255)                     tGTTATGCTGGCGGTGATCATGGTTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGtt                                                                                                                                                                                                            <  1:59691/250‑1 (MQ=255)                              ggCGGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGTCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAggg                                                                                                                                                                                                   <  1:143775/250‑1 (MQ=255)                                gggTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGATGGCGACGCTGATTGGGCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCg                                                                                                                                                                                                <  1:81808/250‑1 (MQ=255)                                        aTGGGTATGATGCTAGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGAGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGAGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAGCCGGCACAACGGCCCTGGTGGTGGGCTGTGGTGGCGTTACTGAAGGGGAGATTACttt                                                                                                                                                                                        >  2:68723/1‑251 (MQ=255)                                                                atGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCgtgg                                                                                                                                                                >  2:140819/1‑250 (MQ=255)                                                                    gCGGTCCGGTTAACAAAGTGTCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAtt                                                                                                                                                            <  1:304822/251‑1 (MQ=255)                                                                                  aatgTGGCCTATGCCTTCATGCTGATGTGCGTTGCTCAGGGTGTTTATCCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGTATGGGGCTGGCGACGCTGATTGGTCTTAAAAATTTCTCCTCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCTATTCCTTTCGCCGTTGCCTATCCGCTGCGTGTCATTCCTTCcatcattgt                                                                                                                                              <  1:361033/248‑4 (MQ=255)                                                                                       ggCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCTGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGt                                                                                                                                          <  1:140819/250‑1 (MQ=255)                                                                                          cTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGAtt                                                                                                                                                                                             <  1:222508/196‑1 (MQ=255)                                                                                          cTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGAACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGAtt                                                                                                                                                                                             >  2:222508/1‑196 (MQ=255)                                                                                             tGCCTTAATGCTGATTTGCGTTGCTCCGGATGTTTATACCGTGGTGGCCATAGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGttt                                                                                                                                   >  2:276631/1‑251 (MQ=255)                                                                                                cTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATAACGTGGTGGCCATCGCTGCCGTGGGTCTTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAAAGGCAATGGTGATGGGCTGTGtt                                                                                                                                                                                                                  >  2:3417/1‑169 (MQ=255)                                                                                                cTTCATGCTGATGTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGGGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGATCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGtt                                                                                                                                                                                                                  <  1:3417/169‑1 (MQ=255)                                                                                                   caggtggctgtGCGTTGCTCAGGTTGTTTATTCCTTGGGGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTgtggtg                                                                                                                             <  2:4357/241‑1 (MQ=255)                                                                                                     tGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCCTGGTCGGTTCAGTTTGTGGTGc                                                                                                                            >  2:108414/1‑250 (MQ=255)                                                                                                            ttGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTg                                                                                                                    >  1:372337/1‑251 (MQ=255)                                                                                                             ggCGTTGCGCAGGGTGTTTATCCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTg                                                                                                                    <  2:306796/249‑1 (MQ=255)                                                                                                               cGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTgcgg                                                                                                                 <  1:113765/251‑1 (MQ=255)                                                                                                                    cTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTg                                                                                                            >  1:4743/1‑251 (MQ=255)                                                                                                                             tttATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTgcg                                                                                                   >  1:363263/1‑251 (MQ=255)                                                                                                                             tttATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTATCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTgcg                                                                                                   >  1:257894/1‑251 (MQ=255)                                                                                                                                       ggtggCCGGCGCTGCCGTGGGTATTTTTTTCCCGCCGCTGGGGATGGGGCGGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCGTCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTTCTATg                                                                                         <  2:398950/251‑1 (MQ=255)                                                                                                                                          ggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGc                                                                                        >  1:159711/1‑249 (MQ=255)                                                                                                                                             cATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCt                                                                                   >  1:93361/1‑251 (MQ=255)                                                                                                                                              aTCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTg                                                                                  >  1:121003/1‑251 (MQ=255)                                                                                                                                                             ttttggtgcccGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATtg                                                                    <  2:324316/242‑1 (MQ=255)                                                                                                                                                                        gccgcTGGGGATGGGGCTGGCGACGCTGATTGGTCGAAAAAATTTCTACGCAGAAGAGAGCGAAACCGGCAACACGGCACTGGTGATGGGATGAGTTGGCGTTAATGAAGGGGCGATTCCTTTCCCCGCTGCCGATCAGCGGCGTGTCATTACTTCCACCCGGGTAGGTTCAGTTTTAGGCGCAGAAACTGCGGCGCTGGTAGGTGCGCAGTGCTCTGCAGGCTGGGGTGGACAGATTGTGATGCAggtgg                                                        >  2:339701/1‑251 (MQ=255)                                                                                                                                                                                     gggatGGCGGCGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACGGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGGTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGcg                                           <  2:366085/247‑2 (MQ=255)                                                                                                                                                                                      ggCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTg                                          >  1:365502/1‑251 (MQ=255)                                                                                                                                                                                           gCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTta                                     <  2:4743/251‑1 (MQ=255)                                                                                                                                                                                                     ttGGTCGTAAAAATTGCTCCTCAGAAGAGCGCGAAACCGGCAAAACGGCACTTGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGAGCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTTTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCTGGTGGTTGAAGGCAAGCTGGGTTATATCgcagc                            <  1:7320/250‑1 (MQ=255)                                                                                                                                                                                                     ttGGTCGCAACAATTTCTCCGCAGAAGAGAGAGAAACCGGCCAAACGTCACTGGTGCTGGGATGTGGTGGCGTAACTGCAGGGACGATTCCTTTCGCCGCTGCAGATCCGATGAGTGTCATTCCTTCCATCATGGACGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGAAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGGTATATCgaagc                            >  2:271695/1‑250 (MQ=255)                                                                                                                                                                                                          cGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCTGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTgg                       <  1:129187/250‑1 (MQ=255)                                                                                                                                                                                                             aaaaaTTTCTCCGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTg                    <  1:35214/250‑1 (MQ=255)                                                                                                                                                                                                                       ccGCAGAAGAGCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAg           <  1:276631/249‑1 (MQ=255)                                                                                                                                                                                                                              agagCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACgg  <  1:12802/251‑1 (MQ=255)                                                                                                                                                                                                                              agagCGCGAAACCGGCAAAACGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGAGCTGCAGGTGGTTGAAGGCAAGCGGGGTTATATAGCAGCAGTGGCTGGTGGAGCAGGGGGGACgg  >  2:223295/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        GGGATGCAGCAGGGCAGCATTGTTATGCTGGCGGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGG  >  NC_000913/4143040‑4143511

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