Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,009,257	T→C	G327G (GGT→GGC)	yedE →	UPF0394 family sulphur transport domain‑containing inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,009,257	0	T	C	100.0%	148.2 / NA	46	G327G (GGT→GGC)	yedE	UPF0394 family sulphur transport domain‑containing inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/26);  total (20/26)

GGCGTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGTCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCTGTTGTTTGCT  >  NC_000913/2009111‑2009404                                                                                                                                                   |                                                                                                                                                    ggCGTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTgg                                                                                                                                                 <  2:352109/151‑1 (MQ=255) ggCGTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTgg                                                                                                                                                 <  1:499870/151‑1 (MQ=255)    gTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCa                                                                                                                                              <  2:581766/151‑1 (MQ=255)       gAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAAtg                                                                                                                                           >  1:610027/1‑151 (MQ=255)       gAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAAt                                                                                                                                            <  1:162230/150‑1 (MQ=255)       gAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAAt                                                                                                                                            <  2:694178/150‑1 (MQ=255)        aaCCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAAtgt                                                                                                                                          >  2:480811/1‑151 (MQ=255)                aTCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCt                                                                                                                                  >  2:584825/1‑151 (MQ=255)                  cATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCa                                                                                                                                <  2:612078/151‑1 (MQ=255)                                     gcgGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTAc                                                                                                              >  1:632375/1‑150 (MQ=255)                                       ggTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTgg                                                                                                           >  2:564549/1‑151 (MQ=255)                                         tAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTggg                                                                                                          >  2:596420/1‑150 (MQ=255)                                         tAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTggg                                                                                                          >  2:518134/1‑150 (MQ=255)                                           aTTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGgatg                                                                                                       <  1:176464/151‑1 (MQ=255)                                                  gtTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATCTCGc                                                                                                >  1:3147/1‑151 (MQ=255)                                                      aCTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCg                                                                                            <  1:518096/151‑1 (MQ=255)                                                      aCTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGCGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCg                                                                                            >  1:574435/1‑151 (MQ=255)                                                       cTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTcc                                                                                             <  1:779926/149‑1 (MQ=255)                                                         gTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGgc                                                                                          <  2:469754/150‑1 (MQ=255)                                                            tGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCt                                                                                       <  2:303140/150‑1 (MQ=255)                                                                      aTCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACt                                                                            <  1:503387/151‑1 (MQ=255)                                                                            ctggctggCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGAc                                                                       >  2:432390/1‑150 (MQ=255)                                                                            ctggctggCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGGGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGa                                                                        >  1:500769/1‑149 (MQ=255)                                                                               gctggcGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATGTCGCTCCGGCGCTGGCCACCGACTGGGACaa                                                                     <  1:374253/149‑1 (MQ=255)                                                                                  ggcggcTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATc                                                                 >  1:496318/1‑150 (MQ=255)                                                                                  ggcggcTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCa                                                                <  1:48702/151‑1 (MQ=255)                                                                                         gCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACctgct                                                         >  1:108412/1‑151 (MQ=255)                                                                                            aaaCCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGa                                                       >  2:795468/1‑150 (MQ=255)                                                                                             aaCCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGaa                                                      >  1:700148/1‑150 (MQ=255)                                                                                                cGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAAcc                                                  <  1:365645/151‑1 (MQ=255)                                                                                                 ggCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCt                                                 >  2:539279/1‑151 (MQ=255)                                                                                                  gCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCtt                                                >  2:332113/1‑151 (MQ=255)                                                                                                     ggATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCttt                                               <  2:125830/149‑1 (MQ=255)                                                                                                          tACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCg                                         >  2:303060/1‑150 (MQ=255)                                                                                                            ccGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATg                                      <  2:241488/151‑1 (MQ=255)                                                                                                                gcgGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTg                                  <  1:640873/151‑1 (MQ=255)                                                                                                                   gTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGcc                               <  1:552806/151‑1 (MQ=255)                                                                                                                     aGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCct                              <  2:108435/150‑1 (MQ=255)                                                                                                                      gAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCctgc                            <  2:542847/151‑1 (MQ=255)                                                                                                                             aGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGAc                     <  1:810835/151‑1 (MQ=255)                                                                                                                              ggTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACa                    <  2:118836/151‑1 (MQ=255)                                                                                                                                tGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACat                   >  1:33973/1‑150 (MQ=255)                                                                                                                                 gCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACatat                 <  2:252807/151‑1 (MQ=255)                                                                                                                                        tggtggGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCtgt          >  1:109010/1‑151 (MQ=255)                                                                                                                                         ggtggGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCtgtt         <  1:236226/151‑1 (MQ=255)                                                                                                                                         ggtggGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCtgt          <  2:830225/150‑1 (MQ=255)                                                                                                                                             ggtcggccTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCTGTTGttt     <  1:609724/151‑1 (MQ=255)                                                                                                                                              gtcggccTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCTGTTGTTtg    <  1:332090/151‑1 (MQ=255)                                                                                                                                                cggccTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCTGTTGTTtgct  <  1:527082/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GGCGTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGTCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCTGTTGTTTGCT  >  NC_000913/2009111‑2009404

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