Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,933,617	Δ1 bp	coding (72/648 nt)	fucA ←	L‑fuculose‑1‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,933,614	0	C	.	100.0%	163.0 / NA	35	coding (75/648 nt)	fucA	L‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (20/15);  total (20/15)

TCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTGAATTAGCCGTTTAA  >  NC_000913/2933471‑2933759                                                                                                                                                |                                                                                                                                                  tCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGtt                                                                                                                                           >  1:417905/1‑151 (MQ=255) tCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGtt                                                                                                                                           >  2:372278/1‑151 (MQ=255)  cGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGtt                                                                                                                                           <  2:412391/150‑1 (MQ=255)   gCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTc                                                                                                                                          >  2:242639/1‑150 (MQ=255)     ttGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGt                                                                                                                                       >  2:15445/1‑151 (MQ=255)                      tcctCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTc                                                                                                                      >  2:564017/1‑151 (MQ=255)                        ctCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGc                                                                                                                                >  1:434051/1‑139 (MQ=255)                        ctCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGc                                                                                                                                <  2:434056/139‑1 (MQ=255)                        ctCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTcc                                                                                                                     >  1:274300/1‑150 (MQ=255)                                        ccATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCaataat                                                                                                    >  1:337606/1‑151 (MQ=255)                                          aTCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCt                                                                                                  >  2:84623/1‑151 (MQ=255)                                             aaTAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGAc                                                                                               >  2:424589/1‑151 (MQ=255)                                                 aaGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGc                                                                                           <  1:493523/151‑1 (MQ=255)                                                 aaGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGc                                                                                           >  2:393140/1‑151 (MQ=255)                                                         atGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTAt                                                                                   <  2:366414/151‑1 (MQ=255)                                                           gCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTAttt                                                                                 <  1:182555/151‑1 (MQ=255)                                                                     cAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCAtt                                                                        <  2:295402/150‑1 (MQ=255)                                                                      aGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAg                                                                      <  1:66061/151‑1 (MQ=255)                                                                           ttCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTAcc                                                                 >  2:77913/1‑151 (MQ=255)                                                                                  ggAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACctctctct                                                          <  2:166008/151‑1 (MQ=255)                                                                                            tAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCaaaa                                                 >  2:357089/1‑150 (MQ=255)                                                                                                 gTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGc                                           <  1:628898/151‑1 (MQ=255)                                                                                                        gCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAAtg                                    >  2:51251/1‑151 (MQ=255)                                                                                                          atcccatccTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAAtgtt                                  >  1:91645/1‑151 (MQ=255)                                                                                                              catccTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTc                              >  2:88490/1‑151 (MQ=255)                                                                                                               atccTGATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTcc                             <  2:274300/151‑1 (MQ=255)                                                                                                                    gATAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTc                         >  2:646126/1‑150 (MQ=255)                                                                                                                      tAACGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTcac                       >  2:303980/1‑150 (MQ=255)                                                                                                                        aCGTACACTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACt                    >  2:409307/1‑151 (MQ=255)                                                                                                                            acacTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTg                <  2:471030/151‑1 (MQ=255)                                                                                                                            acacTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTg                <  1:13528/151‑1 (MQ=255)                                                                                                                             cacTGACGTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTGa               <  2:18151/151‑1 (MQ=255)                                                                                                                                    gTTCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCTGGTCATTTCCAGGCAAGTGTCAAGAATCTGACGATCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTGAAt             <  1:661416/146‑1 (MQ=255)                                                                                                                                     ttCCCCGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTGAATTAGCCg       >  1:218738/1‑151 (MQ=255)                                                                                                                                          cGCTGT‑CCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTGAATTAGCCGTTTaa  <  2:305956/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  TCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAAAACAGGGCAATAATGTTGTTCCTTTCACACTATTGAATTAGCCGTTTAA  >  NC_000913/2933471‑2933759

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