Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,622,019	C→T	P84L (CCG→CTG)	yjjJ →	putative protein kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,019	0	C	T	100.0%	176.9 / NA	48	P84L (CCG→CTG)	yjjJ	putative protein kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/18);  total (30/18)

AGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTG  >  NC_000913/4621873‑4622162                                                                                                                                                   |                                                                                                                                                agaagaTCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAg                                                                                                                                             <  1:558490/151‑1 (MQ=255)   aagaTCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAggg                                                                                                                                           >  1:445352/1‑151 (MQ=255)    agaTCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAggg                                                                                                                                           >  1:226528/1‑150 (MQ=255)        cGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGtt                                                                                                                                      >  1:765859/1‑151 (MQ=255)                tCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGt                                                                                                                               >  2:276023/1‑150 (MQ=255)                tCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGt                                                                                                                               <  1:29860/150‑1 (MQ=255)                    tttGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAAc                                                                                                                            >  2:324380/1‑149 (MQ=255)                     ttGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAg                                                                                                                          <  2:206390/150‑1 (MQ=255)                        gTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGgcgc                                                                                                                      >  2:350555/1‑151 (MQ=255)                          aaaGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAg                                                                                                                    <  1:301315/151‑1 (MQ=255)                            aGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGgcgc                                                                                                                      >  1:302473/1‑147 (MQ=255)                                 gggCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcg                                                                                                              >  1:218718/1‑150 (MQ=255)                                  ggCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGAATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcgac                                                                                                            >  1:1551/1‑151 (MQ=255)                                          gATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGGGGACGATACCGGAAAGGGGCATAAATTTGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGCCTGGTAACAGGCGCAGATGGCGACGCACAGTg                                                                                                    >  1:183179/1‑151 (MQ=255)                                           atatGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTgg                                                                                                   >  1:370576/1‑151 (MQ=255)                                           atatGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTgg                                                                                                   >  2:419253/1‑151 (MQ=255)                                             atGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGt                                                                                                  <  2:445386/150‑1 (MQ=255)                                                  ctgctgCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGCTACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATg                                                                                            >  1:411081/1‑151 (MQ=255)                                                       gcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGttt                                                                                        <  1:443694/150‑1 (MQ=255)                                                       gcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGtt                                                                                         <  2:272163/149‑1 (MQ=255)                                                       gcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTg                                                                                       <  1:372849/151‑1 (MQ=255)                                                              tATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTgg                                                                                 <  2:150967/150‑1 (MQ=255)                                                              tATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTgg                                                                                 >  1:91463/1‑150 (MQ=255)                                                                   cgGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTg                                                                           <  2:827760/151‑1 (MQ=255)                                                                     gAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGAc                                                                         >  2:692384/1‑151 (MQ=255)                                                                            gCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACCAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGAtctc                                                                  >  1:656519/1‑151 (MQ=255)                                                                             cGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTTGTTTGATGGTTTGCCCTGGTATTTGACCGAtctc                                                                  <  2:216698/150‑1 (MQ=255)                                                                                   cccGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGc                                                           >  1:587114/1‑151 (MQ=255)                                                                                            cGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCtttt                                                  <  2:472495/151‑1 (MQ=255)                                                                                             gggTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCttttt                                                 <  1:42678/151‑1 (MQ=255)                                                                                              ggTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCttttt                                                 >  1:222201/1‑150 (MQ=255)                                                                                                 ggACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAgg                                              <  1:742274/150‑1 (MQ=255)                                                                                                   aCGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGgcg                                           >  2:40234/1‑151 (MQ=255)                                                                                                         ccGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTg                                     >  1:174809/1‑151 (MQ=255)                                                                                                              aaGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGc                                 >  1:679021/1‑150 (MQ=255)                                                                                                                ggCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAgg                              >  1:508205/1‑151 (MQ=255)                                                                                                                ggCGCATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAgg                              >  2:12202/1‑151 (MQ=255)                                                                                                                 gcgcATAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGa                             <  2:181208/151‑1 (MQ=255)                                                                                                                     aTAAATTCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGtt                         >  2:233442/1‑151 (MQ=255)                                                                                                                          ttCGCCGACATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCcg                    >  2:818077/1‑151 (MQ=255)                                                                                                                                 aCATCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACt             >  1:494807/1‑151 (MQ=255)                                                                                                                                   aTCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGa           <  2:494843/151‑1 (MQ=255)                                                                                                                                   aTCCGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGa           <  1:656294/151‑1 (MQ=255)                                                                                                                                     ccGGTTGTGCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAAt         >  1:600270/1‑151 (MQ=255)                                                                                                                                          tgtgCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctga     <  1:305942/150‑1 (MQ=255)                                                                                                                                           gtgCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctgact   >  2:130374/1‑151 (MQ=255)                                                                                                                                            tgCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctgactg  >  1:35657/1‑151 (MQ=255)                                                                                                                                            tgCTGGCTGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAAAAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGCTAGCAGCGCAACTGAATctgactg  >  1:266064/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                AGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTG  >  NC_000913/4621873‑4622162

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