Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	560,293	(T)7→6	coding (597/2604 nt)	sfmD →	putative fimbrial usher protein SfmD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	560,287	0	T	.	100.0%	88.9 / NA	21	coding (591/2604 nt)	sfmD	putative fimbrial usher protein SfmD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (6/15);  total (6/15)

AGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTC  >  NC_000913/560209‑560360                                                                               |                                                                          agagTGGGATGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTAttt                                                                         <  2:122062/81‑1 (MQ=255)  gagTGGGATGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTAtttt                                                                        <  1:122063/81‑1 (MQ=255)     tGGGATGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAGCGATAGCTAtttttta                                                                     <  2:122064/81‑2 (MQ=255)         aTGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGa                                                                <  1:122065/81‑1 (MQ=255)            aaGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCt                                                             <  1:122066/81‑1 (MQ=255)             aGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTc                                                            <  1:122067/81‑1 (MQ=255)                       gCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGAtt                                                   >  1:122112/1‑80 (MQ=255)                          ttACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAata                                               >  2:122111/1‑81 (MQ=255)                              tCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTgg                                           >  1:122100/1‑81 (MQ=255)                                cTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGcc                                         <  1:122068/81‑1 (MQ=255)                                   aaTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGt                                      <  1:122069/81‑1 (MQ=255)                                       ataATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGc                                   <  1:122071/80‑1 (MQ=255)                                       ataATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGc                                   <  2:122070/80‑1 (MQ=255)                                        taATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGt                                 <  2:122072/81‑1 (MQ=255)                                         aaTCTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGtt                                >  1:122102/1‑81 (MQ=255)                                          aTTTCACCGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGtt                                <  1:122073/80‑1 (MQ=255)                                                 cGGTAACAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGcaac                        <  1:122074/81‑1 (MQ=255)                                                      aCAGAGGTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATg                    <  2:122075/80‑1 (MQ=255)                                                            gTAATGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTg             <  1:122076/81‑1 (MQ=255)                                                              aaTGGTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTgg            >  2:122104/1‑80 (MQ=255)                                                                  gTAACGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTa       <  2:122077/81‑1 (MQ=255)                                                                    aaCGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt     >  2:122108/1‑81 (MQ=255)                                                                      cGATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtttt   <  1:122078/81‑1 (MQ=255)                                                                       gATAGCTA‑TTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc  <  1:122079/81‑1 (MQ=255)                                                                               |                                                                          AGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTC  >  NC_000913/560209‑560360

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