Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,237,668	T→C	L12P (CTG→CCG)	yjbH →	YjbH family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,237,668	0	T	C	100.0%	189.4 / NA	54	L12P (CTG→CCG)	yjbH	YjbH family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/25);  total (29/25)

CAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGAACTAT  >  NC_000913/4237532‑4237804                                                                                                                                         |                                                                                                                                         cagcaTTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCggg                                                                                                                                       <  2:484424/140‑1 (MQ=255)       tatttatGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAg                                                                                                                                 >  2:199449/1‑140 (MQ=255)         tttatGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGca                                                                                                                               <  1:469280/140‑1 (MQ=255)              gTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCt                                                                                                                          >  1:424366/1‑140 (MQ=255)                 ggCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGc                                                                                                                       >  1:193160/1‑140 (MQ=255)                  gCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCt                                                                                                                      <  2:583675/140‑1 (MQ=255)                  gCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCt                                                                                                                      >  1:1006893/1‑140 (MQ=255)                   ccTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTa                                                                                                                     >  2:484447/1‑140 (MQ=255)                     tCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGACGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACg                                                                                                                   <  2:1006893/140‑1 (MQ=255)                       gcgGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGc                                                                                                                 <  2:621062/140‑1 (MQ=255)                       gcgGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGc                                                                                                                 <  1:484447/140‑1 (MQ=255)                       gcgGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGc                                                                                                                 <  2:1035840/140‑1 (MQ=255)                                 tCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATc                                                                                                       >  2:222064/1‑140 (MQ=255)                                     gAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTgc                                                                                                   >  1:465515/1‑140 (MQ=255)                                      aGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTgcg                                                                                                  <  2:942462/140‑1 (MQ=255)                                         gagaCGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGccc                                                                                               >  1:704517/1‑140 (MQ=255)                                          agaCGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCa                                                                                              >  1:987668/1‑140 (MQ=255)                                                      cccTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGc                                                                                  <  2:424366/140‑1 (MQ=255)                                                         ttAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGt                                                                               <  2:465515/140‑1 (MQ=255)                                                                cGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGAttt                                                                        >  1:236961/1‑140 (MQ=255)                                                                  aCATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCg                                                                      >  2:343087/1‑140 (MQ=255)                                                                  aCATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCg                                                                      >  2:161944/1‑140 (MQ=255)                                                                        ttcAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCg                                                                <  1:233972/140‑1 (MQ=255)                                                                        ttcAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCg                                                                >  2:688877/1‑140 (MQ=255)                                                                         tcAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGt                                                               <  2:926377/140‑1 (MQ=255)                                                                             aCTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGa                                                           >  1:807766/1‑140 (MQ=255)                                                                               tctGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGatt                                                         >  2:611673/1‑140 (MQ=255)                                                                                         cGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACg                                               >  1:1046810/1‑140 (MQ=255)                                                                                                   aataatGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACcgcgc                                     >  1:19355/1‑140 (MQ=255)                                                                                                    ataatGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACcgcgcg                                    <  2:884979/140‑1 (MQ=255)                                                                                                     taatGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGt                                   >  1:357205/1‑140 (MQ=255)                                                                                                      aatGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTa                                  <  1:660807/140‑1 (MQ=255)                                                                                                      aatGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTa                                  <  1:222064/140‑1 (MQ=255)                                                                                                      aatGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCCCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTa                                  <  2:743000/140‑1 (MQ=255)                                                                                                        tGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATg                                >  2:922744/1‑140 (MQ=255)                                                                                                        tGAAAAAAAGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATg                                >  1:749603/1‑140 (MQ=255)                                                                                                          aaaaaaaGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTa                                  >  1:954625/1‑136 (MQ=255)                                                                                                          aaaaaaaGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTa                                  <  2:954625/136‑1 (MQ=255)                                                                                                          aaaaaaaGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGc                              >  1:752489/1‑140 (MQ=255)                                                                                                          aaaaaaaGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGc                              >  1:581015/1‑140 (MQ=255)                                                                                                               aaGACATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACggg                         <  1:477092/140‑1 (MQ=255)                                                                                                                   cATCTGCTTAGCTTACTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAgg                     <  1:107177/140‑1 (MQ=255)                                                                                                                        gcttagcttaCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGt                <  2:236961/140‑1 (MQ=255)                                                                                                                        gcttagcttaCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGt                <  1:967619/140‑1 (MQ=255)                                                                                                                            agcttaCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGc                              >  2:1045908/1‑122 (MQ=255)                                                                                                                            agcttaCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGc                              <  1:1045908/122‑1 (MQ=255)                                                                                                                              cttaCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTc          >  2:157947/1‑140 (MQ=255)                                                                                                                                 aCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGt                >  2:983579/1‑131 (MQ=255)                                                                                                                                 aCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGt                <  1:983579/131‑1 (MQ=255)                                                                                                                                 aCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGa       >  1:372579/1‑140 (MQ=255)                                                                                                                                 aCTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGa       >  1:190475/1‑140 (MQ=255)                                                                                                                                  cTGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGaa      <  1:611673/140‑1 (MQ=255)                                                                                                                                   tGGCGCCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGAAc     >  2:717366/1‑140 (MQ=255)                                                                                                                                      cgcCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGAACTAt  <  1:922744/140‑1 (MQ=255)                                                                                                                                      cgcCGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGAACTAt  <  2:94643/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         CAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACGCCCACCGCGCGTATGGCACGGGAAGGGGAGTTGAGTCTGAACTAT  >  NC_000913/4237532‑4237804

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