Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,073,698	A→C	intergenic (‑127/‑41)	yihU ← / → yihV	3‑sulfolactaldehyde reductase/6‑deoxy‑6‑sulfofructose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,073,698	0	A	C	100.0%	177.6 / NA	49	intergenic (‑127/‑41)	yihU/yihV	3‑sulfolactaldehyde reductase/6‑deoxy‑6‑sulfofructose kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/28);  total (21/28)

CATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGAAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGT  >  NC_000913/4073569‑4073824                                                                                                                                  |                                                                                                                               cATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAAt                                                                                                                             <  1:107478/133‑1 (MQ=255)  aTAACCCTCTCCTGATTACAGTGATGTCACATTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAAtt                                                                                                                            <  1:161752/133‑1 (MQ=255)          ctcCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGaa                                                                                                                    >  1:106338/1‑133 (MQ=255)           tcCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAAc                                                                                                                   >  1:506145/1‑133 (MQ=255)                    cAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCaa                                                                                                          <  1:365387/133‑1 (MQ=255)                     aGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCaaa                                                                                                         >  2:376792/1‑133 (MQ=255)                     aGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCaaa                                                                                                         >  2:364526/1‑133 (MQ=255)                       ttATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAAt                                                                                                       >  2:91340/1‑133 (MQ=255)                       ttATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAAt                                                                                                       >  1:23581/1‑133 (MQ=255)                          tGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGtt                                                                                                    >  1:405915/1‑133 (MQ=255)                                  tttGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAAt                                                                                            <  2:106338/133‑1 (MQ=255)                                     gTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGAc                                                                                         >  2:60993/1‑133 (MQ=255)                                     gTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGAc                                                                                         >  2:452345/1‑133 (MQ=255)                                     gTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGGGCAATTAACATGTACGCATCCCAAACTGTTTAAGGAATGAc                                                                                         >  2:215267/1‑133 (MQ=255)                                      tCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGAcc                                                                                        <  2:23581/133‑1 (MQ=255)                                               aCATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCgt                                                                               <  1:452345/133‑1 (MQ=255)                                                cATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCgtg                                                                              <  1:60993/133‑1 (MQ=255)                                                   gctttgcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCGAAAGGTTGAAGGAATGGCCATGATTCGTGTTg                                                                           >  2:184244/1‑133 (MQ=255)                                                   gctttgcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTg                                                                           >  1:380787/1‑133 (MQ=255)                                                    ctttgcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGc                                                                          <  1:12379/133‑1 (MQ=255)                                                     tttgcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCt                                                                         <  1:115469/133‑1 (MQ=255)                                                      ttgcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCtt                                                                        <  2:380787/133‑1 (MQ=255)                                                       tgcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTtg                                                                       >  1:566571/1‑133 (MQ=255)                                                       tgcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTtg                                                                       >  1:75654/1‑133 (MQ=255)                                                        gcttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTtgt                                                                      <  1:364526/133‑1 (MQ=255)                                                          ttGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTgtgt                                                                    <  2:405915/133‑1 (MQ=255)                                                             tCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCAC‑TTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGt                                                                >  2:36907/1‑133 (MQ=255)                                                                 tttttGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTtgt                                                                      >  1:15873/1‑124 (MQ=255)                                                                 tttttGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTtgt                                                                      <  2:15873/124‑1 (MQ=255)                                                                   tttGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAAc                                                           <  1:549528/133‑1 (MQ=255)                                                                         cGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGAt                                                     <  2:432337/133‑1 (MQ=255)                                                                            aTTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGa                                                  <  2:506145/133‑1 (MQ=255)                                                                             tttGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGAt                                                 <  1:376792/133‑1 (MQ=255)                                                                             tttGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGAt                                                 <  1:91340/133‑1 (MQ=255)                                                                             tttGTAATCTATCGTCAAAAAAGTGACAGCAGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGAt                                                 <  1:184244/133‑1 (MQ=255)                                                                                      tATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACGGTAATGGATGGCATCTAt                                        >  1:221966/1‑133 (MQ=255)                                                                                         cGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTCGGTATAACCGTAATGGATCGCATCTATTAc                                     >  1:268278/1‑133 (MQ=255)                                                                                         cGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTAc                                     >  1:110286/1‑133 (MQ=255)                                                                                                    tgccAGCCGTCACTTTTTAAACATTTGGTGCAATTAACATGAACGCATCCCAAAATG‑TTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTa                         <  1:387265/130‑1 (MQ=255)                                                                                                       cAGCCGTCACTTTTTCAAGAATAGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTAcc                       <  1:215267/133‑1 (MQ=255)                                                                                                             tCACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGga                 <  2:75654/133‑1 (MQ=255)                                                                                                              cACTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGgag                <  2:378008/133‑1 (MQ=255)                                                                                                                cTTTTTAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGagag              <  2:110286/133‑1 (MQ=255)                                                                                                                acttttAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGagag              <  1:36907/131‑1 (MQ=255)                                                                                                                   tttAAACAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCgg           <  2:141595/133‑1 (MQ=255)                                                                                                                     tAAACAATTGGTGCAATTGAAAGGAACGCACCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTa         <  1:214903/133‑1 (MQ=255)                                                                                                                      aaaCAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTGCGTGGAAAGGTGACGGACGGAGAGCGGTaa        >  1:190144/1‑133 (MQ=255)                                                                                                                       aaCAATTGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTaaa       <  1:85894/133‑1 (MQ=255)                                                                                                                            ttGGTGCAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATATATTACGTGGAATGGTAACCGACGGGGAGCGGTAAATACGt  >  2:214903/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                               CATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGAAATTAAAATGAACGCATCCCAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGT  >  NC_000913/4073569‑4073824

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