Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,139,583	G→A	R286H (CGT→CAT)	yghU →	putative S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,139,583	0	G	A	100.0%	150.6 / NA	44	R286H (CGT→CAT)	yghU	putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/23);  total (21/23)

GCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCA  >  NC_000913/3139448‑3139718                                                                                                                                        |                                                                                                                                        gCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCa                                                                                                                                      <  2:349999/139‑1 (MQ=255)                tAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGgt                                                                                                                       >  1:868156/1‑139 (MQ=255)                      aaCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTtcgt                                                                                                                 <  1:98879/139‑1 (MQ=255)                       aCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGGATCAGGGGTa                                                                                                                                >  2:408324/1‑123 (MQ=255)                       aCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                                       <  2:682477/116‑1 (MQ=255)                       aCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                                       >  1:682477/1‑116 (MQ=255)                       aCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTa                                                                                                                                <  1:408324/123‑1 (MQ=255)                       aCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTa                                                                                                                                <  2:600764/123‑1 (MQ=255)                       aCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTa                                                                                                                                >  1:600764/1‑123 (MQ=255)                       aCGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTtcgtc                                                                                                                >  2:897802/1‑139 (MQ=255)                          tCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCgca                                                                                                             <  2:715723/139‑1 (MQ=255)                               cggTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAg                                                                                                        <  2:901736/139‑1 (MQ=255)                                                aTTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAg                                                                                                                              >  1:48878/1‑100 (MQ=255)                                                aTTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAg                                                                                                                              <  2:48878/100‑1 (MQ=255)                                                  tGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAgg                                                                                     <  2:868156/139‑1 (MQ=255)                                                      aaCCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                                       >  1:852390/1‑85 (MQ=255)                                                      aaCCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGt                                                                                                                                 >  1:647056/1‑91 (MQ=255)                                                      aaCCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGt                                                                                                                                 <  2:647056/91‑1 (MQ=255)                                                       aCCGCAACAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTc                                                                                <  1:605224/139‑1 (MQ=255)                                                           cACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGt                                                                                                                                 <  1:615401/86‑1 (MQ=255)                                                           cACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGt                                                                                                                                 >  2:615401/1‑86 (MQ=255)                                                            aCCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgcc                                                                           >  2:899155/1‑139 (MQ=255)                                                            aCCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgcc                                                                           >  1:277901/1‑139 (MQ=255)                                                                    aCCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                   <  2:277901/139‑1 (MQ=255)                                                                        cTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                               >  2:4684/1‑139 (MQ=255)                                                                              gAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGca                                                         >  1:890505/1‑139 (MQ=255)                                                                                   gTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCa                                                                                                                                      <  1:521845/57‑1 (MQ=255)                                                                                   gTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCa                                                                                                                                      >  2:521845/1‑57 (MQ=255)                                                                                     tGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAgg                                                                                               >  2:849669/1‑94 (MQ=255)                                                                                     tGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAgg                                                                                               <  1:849669/94‑1 (MQ=255)                                                                                         tGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGAt                                              >  1:880624/1‑139 (MQ=255)                                                                                             cgcCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGa                                          <  1:897802/139‑1 (MQ=255)                                                                                              gcCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                         >  1:256997/1‑139 (MQ=255)                                                                                                      gCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCcgc                                 <  2:880624/139‑1 (MQ=255)                                                                                                             aTTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                          <  2:890505/139‑1 (MQ=255)                                                                                                             aTTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                          <  1:899155/139‑1 (MQ=255)                                                                                                              tttCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAg                                                                                                        <  2:654301/60‑1 (MQ=255)                                                                                                              tttCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATc                         <  2:353891/139‑1 (MQ=255)                                                                                                              tttCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGCCGGTGTTCGTCGCAGCAAg                                                                                                        >  1:654301/1‑60 (MQ=255)                                                                                                                   agaCGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGcc                    >  2:88253/1‑139 (MQ=255)                                                                                                                        aaTACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa               >  1:268033/1‑139 (MQ=255)                                                                                                                                   tAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGa                                               >  1:4079/1‑96 (MQ=255)                                                                                                                                   tAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGa                                               <  2:4079/96‑1 (MQ=255)                                                                                                                                     aGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAAtcatca  <  1:4684/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCA  >  NC_000913/3139448‑3139718

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