Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,085,675	Δ1 bp	intergenic (+38/+13)	yiiG → / ← frvR	hypothetical protein/putative frv operon regulatory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,085,675	0	G	.	100.0%	362.0 / NA	79	intergenic (+38/+13)	yiiG/frvR	hypothetical protein/putative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (43/36);  total (43/36)

CGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATC  >  NC_000913/4085604‑4085745                                                                        |                                                                       cGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCg                                                                   <  2:1887424/76‑1 (MQ=255) cGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCg                                                                   <  2:1855340/76‑1 (MQ=255)  gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGa                                                                  >  1:2006241/1‑76 (MQ=255)  gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGa                                                                  >  1:782972/1‑76 (MQ=255)      tggtggATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGAcaca                                                              <  2:1272225/76‑1 (MQ=255)          tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTc                                                          <  1:4478818/76‑1 (MQ=255)          tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTc                                                          <  2:1945554/76‑1 (MQ=255)          tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTc                                                          >  1:2592379/1‑76 (MQ=255)            atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAc                                                        >  2:3214671/1‑76 (MQ=255)             tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAcc                                                       >  1:753981/1‑76 (MQ=255)               aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAcccc                                                     >  1:3506441/1‑76 (MQ=255)               aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAcccc                                                     <  1:2130554/76‑1 (MQ=255)                 tatgATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCtt                                                   <  2:3278554/76‑1 (MQ=255)                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCtt                                                   >  1:3523751/1‑76 (MQ=255)                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCtt                                                   >  1:3964401/1‑76 (MQ=255)                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCtt                                                   >  1:1904183/1‑76 (MQ=255)                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGGTCACCCCtt                                                   >  1:1967314/1‑76 (MQ=255)                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTtg                                                  <  1:4073421/76‑1 (MQ=255)                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGa                                               <  1:2496044/76‑1 (MQ=255)                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGt                                            <  2:773837/76‑1 (MQ=255)                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGt                                            >  1:3696684/1‑76 (MQ=255)                         ccTGCATTAATCGTATTTGATGCTCAATGCAGTGTCTGATGCGGCGC‑ATCCGACACATTTCACCCCTTGTGaattc                                           <  1:3379615/76‑4 (MQ=255)                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAAc                                         >  1:4111152/1‑76 (MQ=255)                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACt                                        >  2:2697499/1‑76 (MQ=255)                             cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                       <  1:831766/76‑1 (MQ=255)                             cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                       <  2:4137181/76‑1 (MQ=255)                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTcc                                      >  2:3841668/1‑76 (MQ=255)                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGaa                                              <  2:3167010/67‑1 (MQ=255)                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGaa                                              >  1:3167010/1‑67 (MQ=255)                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCa                                     >  1:224326/1‑76 (MQ=255)                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCa                                     >  1:650528/1‑76 (MQ=255)                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCa                                     <  2:1255408/76‑1 (MQ=255)                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCAt                                    <  2:2852026/76‑1 (MQ=255)                                tAAACGTTTTGGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCAt                                    <  2:2854647/76‑1 (MQ=255)                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGCTCACCCCTTGTGAAGTAACTCCATaa                                  >  1:890419/1‑76 (MQ=255)                                   aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                 <  2:3863381/76‑1 (MQ=255)                                    cttttttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAAt                                <  2:2216649/74‑1 (MQ=255)                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                             <  2:4243079/76‑1 (MQ=255)                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                             >  2:2759204/1‑76 (MQ=255)                                        tttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCt                            <  2:2593153/76‑1 (MQ=255)                                         ttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                           >  1:3031819/1‑76 (MQ=255)                                         ttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                           <  2:3045192/76‑1 (MQ=255)                                            atgatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCata                        >  2:1953678/1‑76 (MQ=255)                                              gatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                      <  2:942711/76‑1 (MQ=255)                                              gatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                      >  1:780688/1‑76 (MQ=255)                                              gatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                      >  2:4222342/1‑76 (MQ=255)                                              gatAAATGCAGTGGCGGATGCGGCGC‑ATCCGACACAGTTCACCACTTGTGAAGTAACTCCATAAATGCCTCATATc                      <  2:974385/76‑1 (MQ=255)                                               atAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                     >  1:4249895/1‑76 (MQ=255)                                               atAAATGCAGTGTCGGAAGCGGCGC‑AACCGACACAGTTCACCCCTTGTGAAGGAACTCCATAAACGCCTCACATcc                     >  2:2691687/1‑76 (MQ=255)                                                tAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                    >  1:4312587/1‑76 (MQ=255)                                                 aaaaGCAGGGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                   <  1:330489/76‑1 (MQ=255)                                                 aaaTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                   <  1:1835389/76‑1 (MQ=255)                                                 aaaTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                   <  1:1540989/76‑1 (MQ=255)                                                 aaaTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                   <  2:1255792/76‑1 (MQ=255)                                                  aaTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                  >  2:3823677/1‑76 (MQ=255)                                                   aTGCAGTGTCGGATGCGGGGC‑ATCCGACGCAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                 <  1:2326294/76‑1 (MQ=255)                                                   aTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                 >  1:4299029/1‑76 (MQ=255)                                                    tGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                <  2:3743205/76‑1 (MQ=255)                                                      cAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc              >  1:963223/1‑76 (MQ=255)                                                      cAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc              <  2:4111152/76‑1 (MQ=255)                                                       aGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc             >  1:2094298/1‑76 (MQ=255)                                                        gtgtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc            <  1:2561354/76‑1 (MQ=255)                                                          gtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg          >  2:3871845/1‑76 (MQ=255)                                                          gtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg          >  1:1614862/1‑76 (MQ=255)                                                          gtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg          >  1:788464/1‑76 (MQ=255)                                                          gtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg          >  1:647517/1‑76 (MQ=255)                                                           tCGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc         <  1:2759204/76‑1 (MQ=255)                                                            cGGATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa        >  2:1200392/1‑76 (MQ=255)                                                             ggATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa       >  2:2921496/1‑76 (MQ=255)                                                             ggATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa       >  2:2159770/1‑76 (MQ=255)                                                             ggATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACGCCATAAATGCCTCATATCCTGTTAACCCGGCaa       >  1:2956927/1‑76 (MQ=255)                                                             ggATGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGCAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa       >  1:1895750/1‑76 (MQ=255)                                                               aTGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACCCCAAACATCCCTCATACCCTTTTAACACCGCacta     >  2:4673582/1‑76 (MQ=255)                                                                tGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa    >  2:1683195/1‑76 (MQ=255)                                                                tGCGGCGC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa    >  2:1430293/1‑76 (MQ=255)                                                                 gcggcgC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat   >  1:790862/1‑76 (MQ=255)                                                                 gcggcgC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat   <  1:3947588/76‑1 (MQ=255)                                                                 gcggcgC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat   <  1:2624705/76‑1 (MQ=255)                                                                  cggcgC‑ATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc  <  1:4237339/76‑1 (MQ=255)                                                                        |                                                                       CGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATC  >  NC_000913/4085604‑4085745

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