Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,965,021	G→A	Q113* (CAA→TAA)	fleQ ←	transcriptional regulator FleQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,965,021	0	G	A	100.0%	125.3 / NA	37	Q113* (CAA→TAA)	fleQ	transcriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/17);  total (20/17)

TGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATT  >  NC_002947/4964876‑4965165                                                                                                                                                  |                                                                                                                                                 tgcATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGc                                                                                                                                                >  4:71915/1‑148 (MQ=255)        tcTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCt                                                                                                                                         >  8:103033/1‑148 (MQ=255)        tcTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCt                                                                                                                                         >  8:178255/1‑148 (MQ=255)                    gCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTc                                                                                                                             <  8:92078/148‑1 (MQ=255)                     cTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCAttt                                                                                                                              >  5:81046/1‑146 (MQ=255)                      tGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCa                                                                                                                           >  4:77716/1‑148 (MQ=255)                           ggCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTg                                                                                                                      <  2:105025/148‑1 (MQ=255)                               cGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGaaa                                                                                                                  <  7:178255/148‑1 (MQ=255)                               cGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGaaa                                                                                                                  <  3:142823/148‑1 (MQ=255)                                     gTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCAccc                                                                                                            >  4:24152/1‑148 (MQ=255)                                              aGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGaa                                                                                                   >  8:4499/1‑148 (MQ=255)                                                     ggAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCtt                                                                                             <  2:19696/147‑1 (MQ=255)                                                     ggAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTc                                                                                            <  3:89449/148‑1 (MQ=255)                                                       aaCAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAACGCACCCGGCGGCGAAGGTCTTCCg                                                                                          >  5:61890/1‑148 (MQ=255)                                                                  ttCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTcggcg                                                                               <  7:117238/148‑1 (MQ=255)                                                                             gACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTc                                                                    <  2:144466/148‑1 (MQ=255)                                                                               ccccGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCAc                                                                  >  8:129345/1‑148 (MQ=255)                                                                                     tCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTaaaa                                                            <  4:89291/148‑1 (MQ=255)                                                                                     tCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTaaaa                                                            <  2:109479/148‑1 (MQ=255)                                                                                       gcgcgcCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAg                                                           >  3:73186/1‑147 (MQ=255)                                                                                         gcgcCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTCGCTTGGCGGCATTTCCAGGTTGGAACGCACCCGGCGGCGATGGTCTTCCGGCAGCTCGGGGGAAGAAATTTCACCTAAAAGCa                                                         >  3:5035/1‑147 (MQ=255)                                                                                             cTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAAcc                                                    >  3:159752/1‑148 (MQ=255)                                                                                                   aTACATTTCGCGGTAGACCTGGGCACGGTGCAGCGCGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAACGCAAAACCggaag                                               >  3:151377/1‑147 (MQ=255)                                                                                                   aTACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCggaag                                               >  5:177219/1‑147 (MQ=255)                                                                                                     aCATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCggaaggaa                                            >  3:69619/1‑148 (MQ=255)                                                                                                       aTTTCGCGGTAGACCTGGGCACTGTTCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAATCACCCGGCGGCGAAGTTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAAc                                           <  4:153332/147‑1 (MQ=255)                                                                                                             cgGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATcc                                     <  7:103033/147‑1 (MQ=255)                                                                                                             cgGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATc                                      >  5:136575/1‑146 (MQ=255)                                                                                                               gTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAg                                  <  1:152132/148‑1 (MQ=255)                                                                                                                   aCCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAg                              >  1:149487/1‑148 (MQ=255)                                                                                                                   aCCTGGGCACGGTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAg                              >  2:199940/1‑148 (MQ=255)                                                                                                                             ggTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCtt                    <  4:58079/148‑1 (MQ=255)                                                                                                                             ggTGCAGCGAGTCCAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCtt                    <  3:71915/148‑1 (MQ=255)                                                                                                                               tGCAGCGAGTCCAGCAATTAGCTGTAGCATGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTaa                  >  5:128436/1‑148 (MQ=255)                                                                                                                                          cAGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAACTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCAACCCAGGTAGCCACTGTCTTaa                  >  7:131108/1‑137 (MQ=255)                                                                                                                                           aGCAATTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAACTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCAACCCAGGTAGCCACTGTCTTAag                 <  8:131108/137‑1 (MQ=255)                                                                                                                                               aTTAGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGAtt  <  1:110605/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATT  >  NC_002947/4964876‑4965165

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