Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_KO_pgi	4,133,781	C→T	intergenic (+275/‑54)	metF → / → katG	5,10‑methylenetetrahydrofolate reductase/catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_KO_pgi	4,133,781	0	C	T	100.0%	167.6 / NA	46	intergenic (+275/‑54)	metF/katG	5,10‑methylenetetrahydrofolate reductase/catalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/23);  total (23/23)

AGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACA  >  NC_000913_3_KO_pgi/4133638‑4133927                                                                                                                                                |                                                                                                                                                   aGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTtctct                                                                                                                                              <  3:75403/150‑1 (MQ=255) aGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTtctct                                                                                                                                              <  1:628260/150‑1 (MQ=255) aGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTtct                                                                                                                                                >  1:119463/1‑148 (MQ=255) aGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTtct                                                                                                                                                <  2:119463/148‑1 (MQ=255)  gCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTtctctc                                                                                                                                             >  4:96636/1‑150 (MQ=255)  gCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTtctctc                                                                                                                                             >  2:505813/1‑150 (MQ=255)             tttttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTgtgt                                                                                                                                  >  1:94920/1‑150 (MQ=255)                      aGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCtgt                                                                                                                                    >  3:416338/1‑139 (MQ=255)                      aGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCtgt                                                                                                                                    <  4:416338/139‑1 (MQ=255)                       gattatGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACgg                                                                                                                        >  2:603360/6‑150 (MQ=255)                        cATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGt                                                                                                                       <  1:505813/150‑1 (MQ=255)                         ttatGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTaa                                                                                                                           <  4:400926/142‑1 (MQ=255)                         ttatGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTaa                                                                                                                           >  3:400926/4‑145 (MQ=255)                                 gAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTACACATGAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTa                                                                                                              <  1:362646/150‑1 (MQ=255)                                      cGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAggg                                                                                                         <  2:460941/150‑1 (MQ=255)                                             aaTCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCg                                                                                                                              >  3:118246/1‑122 (MQ=255)                                             aaTCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCg                                                                                                                              <  4:118246/122‑1 (MQ=255)                                              aTCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACAt                                                                                                 >  2:320538/1‑150 (MQ=255)                                                 aaaaaaGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                                                                                                       <  1:426194/141‑1 (MQ=255)                                                 aaaaaaGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                                                                                                       >  2:426194/1‑141 (MQ=255)                                                         ttaattaaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACg                                                                                      >  4:201014/1‑150 (MQ=255)                                                           aattaaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgatga                                                                                           >  4:24597/1‑143 (MQ=255)                                                           aattaaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgatga                                                                                           <  3:24597/143‑1 (MQ=255)                                                           aattaaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTc                                                                                    <  2:94920/150‑1 (MQ=255)                                                                       tttGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCa                                                                        <  3:201014/150‑1 (MQ=255)                                                                        ttGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCAt                                                                       <  3:96636/150‑1 (MQ=255)                                                                           aTCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAac                                                                    >  2:47221/1‑150 (MQ=255)                                                                                  tctcTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAg                                                                                                      <  4:339790/109‑1 (MQ=255)                                                                                  tctcTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAg                                                                                                      >  3:339790/1‑109 (MQ=255)                                                                                     ctTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCa                                                          >  2:96051/1‑150 (MQ=255)                                                                                        tAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTg                                                       <  1:96051/150‑1 (MQ=255)                                                                                          aCCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGc                                                     <  1:607090/150‑1 (MQ=255)                                                                                            caacaaTATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat                                                                                             >  1:532837/1‑108 (MQ=255)                                                                                            caacaaTATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat                                                                                             <  2:532837/108‑1 (MQ=255)                                                                                                   atGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat                                                                                             >  1:580107/1‑101 (MQ=255)                                                                                                   atGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat                                                                                             <  2:580107/101‑1 (MQ=255)                                                                                                      tAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAaca                                                                   <  1:523584/124‑1 (MQ=255)                                                                                                      tAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAaca                                                                   >  2:523584/1‑124 (MQ=255)                                                                                                      tAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTc                                         >  4:373793/1‑150 (MQ=255)                                                                                                             tcAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATc                                                                                                                               >  3:275867/1‑57 (MQ=255)                                                                                                             tcAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATc                                                                                                                               <  4:275867/57‑1 (MQ=255)                                                                                                             tcAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATaa                                                                     >  1:692359/1‑115 (MQ=255)                                                                                                             tcAACTATCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATaa                                                                     <  2:692359/115‑1 (MQ=255)                                                                                                                   aTCGCATCCGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTc                            >  2:669817/1‑150 (MQ=255)                                                                                                                          ccGTGGATTAATTCAATTATAATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCa                     <  3:373793/150‑1 (MQ=255)                                                                                                                                             taATTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGcaca  >  1:299179/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                   AGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACA  >  NC_000913_3_KO_pgi/4133638‑4133927

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