Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,428,720	T→C	intergenic (‑3/‑215)	ytfB ← / → fklB	OapA family protein/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,428,720	0	T	C	100.0%	154.5 / NA	48	intergenic (‑3/‑215)	ytfB/fklB	OapA family protein/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/31);  total (17/31)

GCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTCCCGCGA  >  NC_000913/4428576‑4428865                                                                                                                                                 |                                                                                                                                                  gCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCCACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCCCGTGCCAGACTTTCTCCAGGGCTGGTTTTAGTTCAAAGCGCCAGGGCATGGCTtcctc                                                                                                                                             >  2:113061/1‑151 (MQ=255)   aGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTcctcct                                                                                                                                           <  2:219023/151‑1 (MQ=255)      aaTCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGc                                                                                                                                         >  1:847083/1‑150 (MQ=255)      aaTCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCt                                                                                                                                        <  1:987012/151‑1 (MQ=255)                  aaTACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGtt                                                                                                                            >  2:802823/1‑151 (MQ=255)                       aaTGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGGGCCAGACTTTCTCCAGGGGTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTc                                                                                                                        >  2:147020/1‑150 (MQ=255)                        aTGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCaa                                                                                                                      <  2:780928/151‑1 (MQ=255)                        aTGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCaa                                                                                                                      <  2:426586/151‑1 (MQ=255)                            cggcAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCAt                                                                                                                  >  1:12817/1‑151 (MQ=255)                              gcAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCAt                                                                                                                  >  1:172488/1‑149 (MQ=255)                                     tGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATg                                                                                                          >  2:399505/1‑150 (MQ=255)                                      gATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGa                                                                                                         <  2:788924/150‑1 (MQ=255)                                         cccccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATa                                                                                                      <  2:99076/150‑1 (MQ=255)                                         cccccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGGGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATa                                                                                                      >  2:980826/1‑150 (MQ=255)                                            ccGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCt                                                                                                  <  2:533291/151‑1 (MQ=255)                                             cgacgaTGCATCGGCGGCAGTGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTa                                                                                                 <  2:483061/151‑1 (MQ=255)                                                        cggcggCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTatgatg                                                                                      >  2:751075/1‑151 (MQ=255)                                                         ggcggcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGt                                                                                     <  1:72913/151‑1 (MQ=255)                                                          gcggcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATaa                                                                                                     <  2:851236/134‑1 (MQ=255)                                                          gcggcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATaa                                                                                                     >  1:851284/1‑134 (MQ=255)                                                            ggcAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCg                                                                                  <  1:518684/151‑1 (MQ=255)                                                                gCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTgg                                                                               <  1:812672/150‑1 (MQ=255)                                                                  gggTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTaa                                                                            <  2:867916/151‑1 (MQ=255)                                                                               gAAAATTGTCCGGCGCGTGCCAGTCTTTCCCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAGCGCTATGATGCCCGTGGTAACCCCGCCTttatt                                                               <  2:969963/151‑1 (MQ=255)                                                                                aaaaTTGTCCGGCGCGTGCCAACCTTTCTCCAGGGTTGGTTTTAGTTCAAACCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATGGCCTAGGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCttt                                                                  <  1:12055/147‑1 (MQ=255)                                                                                  aaTTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATtata                                                            <  1:820111/151‑1 (MQ=255)                                                                                      gTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTTTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATAtt                                                          <  2:774144/149‑1 (MQ=255)                                                                                      gTCCGGCGCGGGCCAGTCTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTAGGATGTCCGTGGTAAACCCGCCTTTATTATATTCa                                                        <  2:172487/151‑1 (MQ=255)                                                                                          ggCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGAATGGCTTCCTCCTGCTCCGCGTCTTTTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGGCCGTGGTAAGCCCGCCTTTATTATATTCATACg                                                    <  1:755078/151‑1 (MQ=255)                                                                                                tGCCAGACTTTCCCCGGGGTGGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTACCGCTATGATGTCCGTGGTAACCCCGCCTTTATTATATTCATACg                                                    <  2:872775/145‑1 (MQ=255)                                                                                                     gACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCt                                          >  1:838681/1‑150 (MQ=255)                                                                                                      aCTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTAt                                        <  1:335146/151‑1 (MQ=255)                                                                                                         ttCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTg                                      <  1:126469/150‑1 (MQ=255)                                                                                                             ccAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCTTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGttt                                   <  1:987542/149‑1 (MQ=255)                                                                                                             ccAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGtt                                    <  2:265339/148‑1 (MQ=255)                                                                                                             ccAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCt                                 >  2:287941/1‑151 (MQ=255)                                                                                                               aGGGTTGGTTTTAGTTCAAAGCGCCGGGGAAGGGCTTCCTCCTGCTCCGCTTCTTGTTCCTCATTCATAGCCTATGAATAATCTAACTTTAGTTTTTCCTTGGTAACTCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTACtt                                <  2:868310/150‑1 (MQ=255)                                                                                                               aGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCttt                               <  2:167511/151‑1 (MQ=255)                                                                                                                gggTTGGTTTTAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCttt                               <  1:1047582/150‑1 (MQ=255)                                                                                                                         ttAGTTCAAATCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGGCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTgg                       >  2:141891/1‑149 (MQ=255)                                                                                                                         ttAGTTCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTggg                      <  1:359522/150‑1 (MQ=255)                                                                                                                             ttCAAAGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTg                 >  2:633243/1‑151 (MQ=255)                                                                                                                                aaaGCGCCCGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGtc               >  2:992729/1‑150 (MQ=255)                                                                                                                                   gcgcCCGGGCATGGCCTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGtcatc            >  2:1038352/1‑150 (MQ=255)                                                                                                                                      cccGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGtcat             >  2:499447/1‑146 (MQ=255)                                                                                                                                      cccGGGCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGtcat             <  1:499473/146‑1 (MQ=255)                                                                                                                                         gggCATGGCTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTccc      <  1:447432/150‑1 (MQ=255)                                                                                                                                            cATGGCTTCCTCCTGCTCCGCGTCTTGCTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTCCCGCGa  <  1:620167/151‑1 (MQ=255)                                                                                                                                               gttttCCCCCTGTTCCGCGTTTTTTTCCTCAATCATAGCTTATAACGAAGCTAACGTTATTATTTCTTTTTTTTTTCCGCCTTTAGTACACTCATACGGGATATTGCTATTGTCCCCTTTTCCCTCGCACTTGTCCTCATTccctcc   <  2:524114/145‑4 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GCAGAAATCCGACCACCAATACAATGGCGGCAATAATGATCCCCCGACGATGCATCGGCGGCAGCGGGTCCATGAAGCGAAAATTGTCCGGCGCGTGCCAGACTTTCTCCAGGGTTGGTTTTAGTTCAAAGCGCCCGGGCATGGTTTCCTCCTGCTCCGCGTCTTGTTCCTCAATCATAGCCTATGAATAAGCTAACGCTATGATGTCCGTGGTAAACCCGCCTTTATTATATTCATACGGGATATTGCTATTGTTTCTTTTTCCCTGGGATTTGTCATCATTCCCGCGA  >  NC_000913/4428576‑4428865

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