Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,460,788	C→T	F324F (TTC→TTT)	gspF →	Type II secretion system protein GspF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,460,788	0	C	T	100.0%	131.2 / NA	36	F324F (TTC→TTT)	gspF	Type II secretion system protein GspF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/18);  total (18/18)

CCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCA  >  NC_000913/3460710‑3460863                                                                               |                                                                            ccAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTtcc                                                                           <  1:594285/81‑1 (MQ=255)   aGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTccgc                                                                         <  2:594285/81‑1 (MQ=255)    gCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTccgcc                                                                        <  2:594286/81‑1 (MQ=255)     cGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTccgccg                                                                       <  1:594288/81‑1 (MQ=255)      gTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCga                                                                      >  2:594338/1‑81 (MQ=255)      gTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCga                                                                      >  2:594328/1‑81 (MQ=255)        cTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCgat                                                                     >  2:594326/1‑80 (MQ=255)           gCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGatgatg                                                                 <  2:594288/81‑1 (MQ=255)            cAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTAATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGc                                                                <  1:594290/81‑1 (MQ=255)             aaaTGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGct                                                               <  2:594290/81‑1 (MQ=255)              aaTGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGctc                                                              >  1:594352/1‑81 (MQ=255)               aTGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGctct                                                             <  2:594291/81‑1 (MQ=255)                  cGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTAc                                                           <  1:594293/80‑1 (MQ=255)                    gCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACAtg                                                        <  2:594293/81‑1 (MQ=255)                       gagaACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACAtgttg                                                     >  2:594359/1‑81 (MQ=255)                       gagaACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACAtggtg                                                     >  1:594342/1‑81 (MQ=255)                         gaACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGc                                                   >  2:594352/1‑81 (MQ=255)                         gaACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGc                                                   >  2:594331/1‑81 (MQ=255)                          aaCGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGcc                                                  >  2:594357/1‑81 (MQ=255)                             gTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCctct                                               >  1:594330/1‑81 (MQ=255)                             gTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCctct                                               >  2:594343/1‑81 (MQ=255)                               tCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTgg                                             >  1:594337/1‑81 (MQ=255)                                cGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGc                                            >  2:594346/1‑81 (MQ=255)                                 gCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGc                                            >  1:594348/1‑80 (MQ=255)                                    aGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGaaa                                        <  1:594295/81‑1 (MQ=255)                                     gggTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGaaaa                                       <  1:594296/81‑1 (MQ=255)                                          aCAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGc                                   <  2:594296/80‑1 (MQ=255)                                           cAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCgg                                 >  1:594341/1‑81 (MQ=255)                                                 tCATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCt                           >  2:594356/1‑81 (MQ=255)                                                  cATCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCt                           <  1:594298/80‑1 (MQ=255)                                                   aTCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTc                          <  2:594298/80‑1 (MQ=255)                                                   aTCTTTCGCTGGAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCg                         >  1:594376/1‑81 (MQ=255)                                                              gAACAAACCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTaa                <  2:594299/79‑1 (MQ=255)                                                                   aaCCGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGt            >  1:594367/1‑78 (MQ=255)                                                                     ccGCAATTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCagag       <  2:594300/81‑1 (MQ=255)                                                                         aaTTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGc   <  1:594302/81‑1 (MQ=255)                                                                          aTTTTTCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCa  <  2:594302/81‑1 (MQ=255)                                                                            ttttCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCg    >  2:594350/1‑77 (MQ=255)                                                                            ttttCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCg    >  2:594354/1‑77 (MQ=255)                                                                               |                                                                            CCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCA  >  NC_000913/3460710‑3460863

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