Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,320,333	T→C	E427E (GAA→GAG)	yhbX ←	putative EptAB family phosphoethanolamine transferase, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,320,333	0	T	C	100.0%	187.4 / NA	52	E427E (GAA→GAG)	yhbX	putative EptAB family phosphoethanolamine transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/28);  total (24/28)

TCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTGAAAC  >  NC_000913/3320208‑3320468                                                                                                                              |                                                                                                                                        tCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGt                                                                                                                            <  1:329107/139‑1 (MQ=255) tCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGt                                                                                                                            <  1:417681/139‑1 (MQ=255)   cGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCg                                                                                                                          <  1:937023/139‑1 (MQ=255)          cGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCAGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAAt                                                                                                                   >  1:983075/1‑139 (MQ=255)                            aGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                                 <  1:995520/139‑1 (MQ=255)                              aaCGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGa                                                                                               <  1:88708/139‑1 (MQ=255)                                   ggCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTg                                                                                                                              <  1:1361441/103‑1 (MQ=255)                                   ggCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTg                                                                                                                              >  2:1361441/1‑103 (MQ=255)                                                  aCATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAACTCAAAAAc                                                                           >  1:832156/1‑139 (MQ=255)                                                       gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaata                                                                                     <  2:411480/124‑1 (MQ=255)                                                       gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaata                                                                                     >  1:411480/1‑124 (MQ=255)                                                         gACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                    >  1:649761/1‑139 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGcc                                                                                                                                 >  1:747414/1‑68 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGcc                                                                                                                                 <  2:747414/68‑1 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                                 >  2:951942/1‑100 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                                 <  1:951942/100‑1 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCa                                                                                >  2:1173746/1‑117 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCa                                                                                <  1:1173746/117‑1 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                          >  2:45330/1‑139 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                          <  2:535699/139‑1 (MQ=255)                                                                   gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                          <  2:983075/139‑1 (MQ=255)                                                                        ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                                 <  1:770199/95‑1 (MQ=255)                                                                        ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                                 >  2:770199/1‑95 (MQ=255)                                                                        ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                                 <  1:238938/95‑1 (MQ=255)                                                                        ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                                 >  2:238938/1‑95 (MQ=255)                                                                         tgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                    <  2:1121749/123‑1 (MQ=255)                                                                         tgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                    >  1:1121749/1‑123 (MQ=255)                                                                          ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTAGGCAAAATACATGACCGAGGCGCGGCGATCTATTAATAATTCAAAAACCTGACCTAGTAAACTATCTGTGTa                                                   <  2:832156/139‑1 (MQ=255)                                                                                ccTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTaatgaat                                             >  1:479982/1‑139 (MQ=255)                                                                                 cTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTaatgaatg                                            <  2:1282779/139‑1 (MQ=255)                                                                                   cccTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                          <  2:649761/139‑1 (MQ=255)                                                                                         cTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCa                                                                                                                       >  1:378531/1‑56 (MQ=255)                                                                                         cTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCa                                                                                                                       <  2:378531/56‑1 (MQ=255)                                                                                           ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCAt                                  <  1:324007/139‑1 (MQ=255)                                                                                           ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAACTGTGTAATGAATGGAGTTGTCAt                                  <  1:928359/139‑1 (MQ=255)                                                                                              tGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGTTCTTTtaata                                                                                     <  2:688337/85‑1 (MQ=255)                                                                                              tGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGTTCTTTtaata                                                                                     >  1:688337/1‑85 (MQ=255)                                                                                               gAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCa                              <  1:553403/139‑1 (MQ=255)                                                                                                          ttcttcttcGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACa                                                                                                                >  2:343757/1‑46 (MQ=255)                                                                                                          ttcttcttcGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACa                                                                                                                <  1:343757/46‑1 (MQ=255)                                                                                                            cttcttcGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCg                 >  2:1351698/1‑139 (MQ=255)                                                                                                            cttcttcGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCg                 >  1:545635/1‑139 (MQ=255)                                                                                                                   gTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                          >  1:82638/1‑107 (MQ=255)                                                                                                                   gTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                          <  2:82638/107‑1 (MQ=255)                                                                                                                       ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAg                                                                                                        >  1:888633/1‑41 (MQ=255)                                                                                                                       ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAg                                                                                                        <  2:888633/41‑1 (MQ=255)                                                                                                                         tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCa                                                                                <  1:1052512/63‑1 (MQ=255)                                                                                                                         tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCa                                                                                >  2:1052512/1‑63 (MQ=255)                                                                                                                         tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                          <  1:21498/85‑1 (MQ=255)                                                                                                                         tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                          >  2:21498/1‑85 (MQ=255)                                                                                                                         tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTGaa    >  1:272571/1‑139 (MQ=255)                                                                                                                           aCGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTGAAac  >  1:295040/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                        TCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTGAAAC  >  NC_000913/3320208‑3320468

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