Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,095,367	C→T	A206T (GCA→ACA)	rhaA ←	L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,367	0	C	T	100.0%	135.4 / NA	39	A206T (GCA→ACA)	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/19);  total (20/19)

TGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAGGCCGAAACGCGACGGCTGGCTTTGCAG  >  NC_000913/4095247‑4095500                                                                                                                         |                                                                                                                                      tGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCt                                                                                                                            <  2:61043/132‑1 (MQ=255) tGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCt                                                                                                                            >  1:61043/1‑132 (MQ=255)           aTTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACgg                                                                                                           >  1:631648/1‑139 (MQ=255)           aTTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACgg                                                                                                           >  2:169129/1‑139 (MQ=255)                      ccGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGGTATGGTGCGCAGGTTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGat                                                                                                <  1:1072501/139‑1 (MQ=255)                          gtAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCt                                                                                            >  2:664685/1‑139 (MQ=255)                                ctctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATa                                                                                      <  1:387007/139‑1 (MQ=255)                                 tctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTc                                                                                                         <  2:854906/119‑1 (MQ=255)                                 tctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTc                                                                                                         >  1:854906/1‑119 (MQ=255)                                   tGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGat                                                                                                <  1:927350/126‑1 (MQ=255)                                   tGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGat                                                                                                >  2:927350/1‑126 (MQ=255)                                   tGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCa                                                                                   >  1:236503/1‑139 (MQ=255)                                       ccAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCg                                                                               <  2:1158211/139‑1 (MQ=255)                                        cAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCgg                                                                              <  1:1142546/139‑1 (MQ=255)                                          aTGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGAcg                                                                                                                        >  2:944142/1‑95 (MQ=255)                                          aTGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGAcg                                                                                                                        <  1:944142/95‑1 (MQ=255)                                          aTGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGat                                                                                                <  2:1131887/119‑1 (MQ=255)                                          aTGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGat                                                                                                >  1:1131887/1‑119 (MQ=255)                                           tGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGAt                                                                           <  1:50949/139‑1 (MQ=255)                                           tGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGAt                                                                           >  2:50949/1‑139 (MQ=255)                                               aaaCAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAg                                                                       <  2:631648/139‑1 (MQ=255)                                                aaCAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGa                                                                      <  1:169129/139‑1 (MQ=255)                                                               aCGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATg                                                       <  2:858864/139‑1 (MQ=255)                                                                 ggCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGgt                                                     >  1:577832/1‑139 (MQ=255)                                                                        atatGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGa                                              >  2:69682/1‑139 (MQ=255)                                                                            ggTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTg                                          <  1:664685/139‑1 (MQ=255)                                                                                  cAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGAcg                                                                                                                        >  1:1030146/1‑55 (MQ=255)                                                                                  cAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGAcg                                                                                                                        <  2:1030146/55‑1 (MQ=255)                                                                                  cAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGcc                                    >  1:1115739/1‑139 (MQ=255)                                                                                   aGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCAc                                                            >  2:403793/1‑114 (MQ=255)                                                                                   aGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCAc                                                            <  1:403793/114‑1 (MQ=255)                                                                                       tttAGCTTCTCGCTGATCACCTCGTCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATa                                                                                      <  2:667501/84‑1 (MQ=255)                                                                                       tttAGCTTCTCGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATa                                                                                      >  1:667501/1‑84 (MQ=255)                                                                                               ctcGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTc                                                                                                         >  1:1050740/1‑57 (MQ=255)                                                                                               ctcGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTc                                                                                                         <  2:1050740/57‑1 (MQ=255)                                                                                               ctcGCTGATCACCTCATCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCAGCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAGGCCGaa                       >  1:741961/1‑139 (MQ=255)                                                                                                              aTCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAGGCCGAAACGCGACggctggct        >  2:430802/1‑139 (MQ=255)                                                                                                              aTCCAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAGGCCGAAACGCGACggctggct        <  1:430802/139‑1 (MQ=255)                                                                                                                ccAGTGCTGTCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAGGCCGAAACGCGACGGCTGGCttt      >  2:948783/1‑139 (MQ=255)                                                                                                                    tgctgtCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAGGCCGAAACGCGACGGCTGGCTTTGCAg  >  1:900941/1‑139 (MQ=255)                                                                                                                     gctgtCAGCAGACGCTGACGCGGGGCGAGACGGTCAACg                                                                                                     >  2:572988/1‑39 (MQ=255)                                                                                                                     gctgtCAGCAGACGCTGACGCGGGGCGAGACGGTCAACg                                                                                                     <  1:572988/39‑1 (MQ=255)                                                                                                                         |                                                                                                                                      TGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAGGCCGAAACGCGACGGCTGGCTTTGCAG  >  NC_000913/4095247‑4095500

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