Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,161,091	T→C	Y19C (TAC→TGC)	ygiQ ←	radical SAM superfamily protein YgiQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,161,091	0	T	C	100.0%	180.2 / NA	51	Y19C (TAC→TGC)	ygiQ	radical SAM superfamily protein YgiQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/28);  total (23/28)

GCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGTACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGT  >  NC_000913/3160956‑3161221                                                                                                                                        |                                                                                                                                   gCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTg                                                                                                                                <  2:890941/140‑1 (MQ=255)          aaGCTTGGGTGATCGACATACGGGCCGCCAGTAACCAAAATGATGCCGCAGCTAT‑CCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGgag                                                                                                                      <  2:1172689/140‑1 (MQ=255)          aaGCTTGGGTGATCGACATACGCGTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGga                                                                                                                       >  1:455119/1‑140 (MQ=255)             cTTGGGTGATCGACATACGCGTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGaa                                                                                                                    <  2:411340/140‑1 (MQ=255)                 ggTGATCGACATACGCGTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAgg                                                                                                                >  1:1031552/1‑140 (MQ=255)                            tACGCGTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCgg                                                                                                     <  1:863092/140‑1 (MQ=255)                               gcgTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTcgcg                                                                                                           >  2:710918/1‑131 (MQ=255)                               gcgTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTcgcg                                                                                                           <  1:710918/131‑1 (MQ=255)                               gcgTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTg                                                                                                  <  2:1209170/140‑1 (MQ=255)                                cgTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTgg                                                                                                 <  1:679271/140‑1 (MQ=255)                                  tCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGAt                                                                                               <  2:536807/140‑1 (MQ=255)                                  tCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGAt                                                                                               >  1:180254/1‑140 (MQ=255)                                    gCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTgg                                                                                                 >  1:1106219/1‑136 (MQ=255)                                    gCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTgg                                                                                                 <  2:1106219/136‑1 (MQ=255)                                     ccAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAg                                                                                            <  1:1071690/140‑1 (MQ=255)                                       aGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTgg                                                                                                 <  1:1198857/133‑1 (MQ=255)                                       aGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTgg                                                                                                 >  2:1198857/1‑133 (MQ=255)                                       aGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAggg                                                                                          <  2:455119/140‑1 (MQ=255)                                       aGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCAGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGAATCAggg                                                                                          <  1:765516/140‑1 (MQ=255)                                                     tgTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCAtt                                                                            >  1:1125066/1‑140 (MQ=255)                                                       tCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGc                                                                                 >  1:1109475/1‑133 (MQ=255)                                                       tCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGc                                                                                 <  2:1109475/133‑1 (MQ=255)                                                       tCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCAtttt                                                                          <  2:1031552/140‑1 (MQ=255)                                                          cAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCt                                                                       >  1:935479/1‑140 (MQ=255)                                                                 cccAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGaa                                                                >  2:85996/1‑140 (MQ=255)                                                                     gCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGa                                                            >  2:973190/1‑140 (MQ=255)                                                                              aTCCATCTCTTCACGAGACATCGGCAAAAACGGTGTCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAaat                                                   >  1:21595/1‑140 (MQ=255)                                                                                   tctcTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCa                                              <  2:180254/140‑1 (MQ=255)                                                                                    ctctTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCaa                                             >  2:427998/1‑140 (MQ=255)                                                                                    ctctTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCa                                              >  1:167089/1‑139 (MQ=255)                                                                                    ctctTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCa                                              <  2:167089/139‑1 (MQ=255)                                                                                        tCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGc                                         >  2:455772/1‑140 (MQ=255)                                                                                         cACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCGGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGcc                                        <  2:617796/140‑1 (MQ=255)                                                                                         cACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGcc                                        <  1:973190/140‑1 (MQ=255)                                                                                            gagaCATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCggg                                     <  1:427998/140‑1 (MQ=255)                                                                                            gagaCATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCggg                                     <  2:24019/140‑1 (MQ=255)                                                                                              gaCATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGa                                   >  2:1302172/1‑140 (MQ=255)                                                                                              gaCATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGa                                   >  2:1096708/1‑140 (MQ=255)                                                                                                 aTCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCggg                                     >  1:1147845/1‑135 (MQ=255)                                                                                                 aTCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCggg                                     <  2:1147845/135‑1 (MQ=255)                                                                                                             ggTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACg                    <  1:1302172/140‑1 (MQ=255)                                                                                                             ggTGCCGGTCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACg                    <  1:1096708/140‑1 (MQ=255)                                                                                                                     tCCAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGaa            >  2:926320/1‑140 (MQ=255)                                                                                                                       cAAAACAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAAt          >  2:477509/1‑140 (MQ=255)                                                                                                                        aaaaCAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATc         <  1:926320/140‑1 (MQ=255)                                                                                                                         aaaCAGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCg        >  1:845904/1‑140 (MQ=255)                                                                                                                             aGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAg    <  2:21595/140‑1 (MQ=255)                                                                                                                             aGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAg    <  2:845904/140‑1 (MQ=255)                                                                                                                             aGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAg    <  1:314758/140‑1 (MQ=255)                                                                                                                             aGGCGGCCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGCTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATACTCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAg    <  2:935479/140‑1 (MQ=255)                                                                                                                              ggcggcCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg   >  2:1258486/1‑140 (MQ=255)                                                                                                                               gcggcCCAGCACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGt  >  1:687870/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                   GCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAGTAACCAAAATGATGTCGCAGCTATCCCAGCCAAGTTGATCCATCTCTTCACGAGACATCGGCAAAAACGGTGCCGGTCCAAAACAGGCGGCCCAGTACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGT  >  NC_000913/3160956‑3161221

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: ‑