Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	504,535	T→C	E54E (GAA→GAG)	fsr ←	fosmidomycin efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	504,535	0	T	C	100.0%	195.4 / NA	52	E54E (GAA→GAG)	fsr	fosmidomycin efflux pump
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/22);  total (30/22)

AAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGC  >  NC_000913/504400‑504666                                                                                                                                        |                                                                                                                                    aaaGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTctg                                                                                                                                 >  1:198294/1‑140 (MQ=255)   aGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTctgac                                                                                                                               >  2:351121/1‑140 (MQ=255)   aGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTctgac                                                                                                                               >  2:869506/1‑140 (MQ=255)                   cAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATaaa                                                                                                               <  2:198294/140‑1 (MQ=255)                    aaCGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAAt                                                                                                              >  1:84463/1‑140 (MQ=255)                    aaCGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAAt                                                                                                              >  1:705367/1‑140 (MQ=255)                            tGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGaa                                                                                                      >  2:645837/1‑140 (MQ=255)                            tGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGaa                                                                                                      >  2:684666/1‑140 (MQ=255)                              gCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATc                                                                                                    <  2:705367/140‑1 (MQ=255)                               cATCGGGTATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCa                                                                                                   <  1:79710/140‑1 (MQ=255)                               cATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCa                                                                                                   <  2:36863/140‑1 (MQ=255)                                    gATATTTATCGGTCCAGTAGCCGACTACTGGATGCAGTAGCGACGAGTCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGAt                                                                                              >  2:901898/1‑140 (MQ=255)                                        tttATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGa                                                                                          >  2:625261/1‑140 (MQ=255)                                               gTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGa                                                                                          >  1:783102/1‑133 (MQ=255)                                               gTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGa                                                                                          <  2:783102/133‑1 (MQ=255)                                               gTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTc                                                                                   >  2:937153/1‑140 (MQ=255)                                               gTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTc                                                                                   >  1:486832/1‑140 (MQ=255)                                                tCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCg                                                                                  >  1:730204/1‑140 (MQ=255)                                                tCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCg                                                                                  <  1:367249/140‑1 (MQ=255)                                                  cAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGtt                                                                                <  1:869506/140‑1 (MQ=255)                                                   aGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTc                                                                               >  1:54215/1‑140 (MQ=255)                                                     tAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTcag                                                                             <  1:625261/140‑1 (MQ=255)                                                     tAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTcag                                                                             <  1:936261/140‑1 (MQ=255)                                                        ccGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTcagcag                                                                          <  1:645837/140‑1 (MQ=255)                                                               cTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTgag                                                                   >  1:1047277/1‑140 (MQ=255)                                                               cTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTgag                                                                   >  1:846427/1‑140 (MQ=255)                                                                  gTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAgg                                                                >  2:189695/1‑140 (MQ=255)                                                                     gCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTg                                                             >  2:524617/1‑140 (MQ=255)                                                                        gTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCGGAATCGGCGATTGGATCGTGTCGTGCAGCAGATGTGAGAGGCTGAtc                                                          <  2:1135143/140‑2 (MQ=255)                                                                        gTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATa                                                          <  2:846427/140‑1 (MQ=255)                                                                        gTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATa                                                          <  2:1047277/140‑1 (MQ=255)                                                                        gTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATa                                                          <  2:54215/140‑1 (MQ=255)                                                                         tAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAg                                                         >  1:874207/1‑140 (MQ=255)                                                                                   gCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGc                                                        <  2:830220/131‑1 (MQ=255)                                                                                   gCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGc                                                        >  1:830220/1‑131 (MQ=255)                                                                                           gAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACg                                       >  1:205540/1‑140 (MQ=255)                                                                                             aGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGAt                                     >  1:1089603/1‑140 (MQ=255)                                                                                                tGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGtt                                  <  2:486832/140‑1 (MQ=255)                                                                                                tGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGtt                                  <  2:730204/140‑1 (MQ=255)                                                                                                 gAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTc                                 >  2:218936/1‑140 (MQ=255)                                                                                                     gTTATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGc                             >  1:408117/1‑140 (MQ=255)                                                                                                       tATCATGCCAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCt                           <  1:189695/140‑1 (MQ=255)                                                                                                              ccAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGt                    >  1:777463/1‑140 (MQ=255)                                                                                                              ccAATCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGt                    <  1:351121/140‑1 (MQ=255)                                                                                                                  tCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGAc                <  2:801388/140‑1 (MQ=255)                                                                                                                  tCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGAc                >  2:1037638/1‑140 (MQ=255)                                                                                                                  tCTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGAc                <  2:127145/140‑1 (MQ=255)                                                                                                                   cTGCATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTAGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACg               >  1:229056/1‑140 (MQ=255)                                                                                                                      cATAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAg            <  1:937153/140‑1 (MQ=255)                                                                                                                        tAAATGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTAGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGcc          <  2:229056/140‑1 (MQ=255)                                                                                                                           aTGTCAGAGAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCcgcg       >  2:413530/1‑140 (MQ=255)                                                                                                                                agagAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGcgcccgc  >  1:865864/1‑140 (MQ=255)                                                                                                                                agagAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGcgcccgc  >  1:657156/1‑140 (MQ=255)                                                                                                                                agagAAAACTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGcgcccgc  >  1:400894/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                    AAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGC  >  NC_000913/504400‑504666

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