Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,163,010	G→A	P158S (CCG→TCG)	plsC ←	1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,163,010	0	G	A	100.0%	187.2 / NA	48	P158S (CCG→TCG)	plsC	1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/23);  total (25/23)

GACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGGTAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTG  >  NC_000913/3162864‑3163146                                                                                                                                                    |                                                                                                                                         gacGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATa                                                                                                                                        <  1:304242/149‑1 (MQ=255)       aaTTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCcg                                                                                                                                >  2:30210/1‑151 (MQ=255)          tGGGGGCAGCATTCCGTCAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGGGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCcgcg                                                                                                                              <  1:26007/150‑1 (MQ=255)           ggcggcAGCATTTCGACAATCCCCAGACCGTTTTGCATGCGATTAAGATTAATCTTATCCGAAGTTGTAGAGACGCACACGGGAATAAGCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCcgcg                                                                                                                              <  2:472741/149‑1 (MQ=255)                      ttCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCg                                                                                                                  <  2:58285/150‑1 (MQ=255)                       tCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCggg                                                                                                                >  2:263311/1‑151 (MQ=255)                       tCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGATGTAGAGACGCACACGGGAATAATCGGGACGC‑ACGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCggg                                                                                                                >  2:423071/1‑151 (MQ=255)                                          tgtgCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATcc                                                                                              >  1:673483/1‑150 (MQ=255)                                           gtgCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCa                                                                                             <  2:658273/150‑1 (MQ=255)                                                 tCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGa                                                                                       >  2:563251/1‑150 (MQ=255)                                                   gataaagattagTTTAGTACGAGGGTGTAGCGACCCACTGGGGAATAA‑AGGCATGCTCCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCGCGGCTGCGGGTTCCTTCCGGGAACATCCAg                                                                                            <  2:347391/143‑1 (MQ=255)                                                    attaagattatTCTTATTCGAATTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATTCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATa                                                                                   <  2:512799/151‑1 (MQ=255)                                                    attaagattaaTCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAAt                                                                                    <  1:585234/150‑1 (MQ=255)                                                     ttaagattaaTCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATAc                                                                                  >  2:560380/1‑151 (MQ=255)                                                      taagattaaTCTTATTCGAAGTTGTAGAGACGCACACGGAAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATAc                                                                                  <  1:213256/150‑1 (MQ=255)                                                        agattaaTCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGc                                                                               <  1:372516/151‑1 (MQ=255)                                                             aaTCTTATTCGAAGTTGTAGAGCCGCACCGGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGtttt                                                                          <  2:534863/151‑1 (MQ=255)                                                                  tATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGa                                                                      >  2:400976/1‑150 (MQ=255)                                                                     tCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTg                                                                  >  2:257175/1‑151 (MQ=255)                                                                           ttGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCac                                                            >  2:280094/1‑151 (MQ=255)                                                                              tAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCact                                                           <  1:193346/149‑1 (MQ=255)                                                                                  gaCGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTcc                                                     <  1:156867/151‑1 (MQ=255)                                                                                     gCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGc                                                   >  1:343737/1‑150 (MQ=255)                                                                                      cgcacGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGc                                                   >  2:479439/3‑149 (MQ=255)                                                                                      cacaCGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCa                                                  >  1:212125/1‑150 (MQ=255)                                                                                       acacGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCTGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCa                                                  <  2:530073/149‑1 (MQ=255)                                                                                       acacGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAAt                                                >  1:328476/1‑151 (MQ=255)                                                                                         acGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATgg                                              >  2:225709/1‑151 (MQ=255)                                                                                           gggAATAATCGGCACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGt                                             <  1:632159/150‑1 (MQ=255)                                                                                            ggAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGc                                           <  1:599321/151‑1 (MQ=255)                                                                                             gAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGcc                                          >  2:186209/1‑151 (MQ=255)                                                                                                taatCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTg                                       >  1:546371/1‑151 (MQ=255)                                                                                                 aatCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTg                                       >  1:201698/1‑150 (MQ=255)                                                                                                     gggACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCtt                                  <  1:506420/151‑1 (MQ=255)                                                                                                      ggACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCttt                                 >  2:294879/1‑151 (MQ=255)                                                                                                       gACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTa                                <  1:227003/151‑1 (MQ=255)                                                                                                        aCGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCttt                                 >  1:196507/1‑149 (MQ=255)                                                                                                        aCGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCttt                                 >  2:142151/1‑149 (MQ=255)                                                                                                        aCGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGAGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCttt                                 >  1:61979/1‑149 (MQ=255)                                                                                                             ccgccgCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCg                           <  1:536235/151‑1 (MQ=255)                                                                                                                  gcAATTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGt                      <  1:492041/151‑1 (MQ=255)                                                                                                                    aaTTGCCGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGttt                    <  2:253021/151‑1 (MQ=255)                                                                                                                         ccGCGTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTc               >  1:498252/1‑151 (MQ=255)                                                                                                                           gcgTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCg              >  2:466048/1‑150 (MQ=255)                                                                                                                           gcgTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCg              >  1:442551/1‑150 (MQ=255)                                                                                                                           gcgTGAAATGCTCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGa             <  2:297359/151‑1 (MQ=255)                                                                                                                                      tCCAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTg  >  1:12647/1‑151 (MQ=255)                                                                                                                                       ccAGTCTTGAACGATAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTg  <  2:486093/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                         GACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGACGC‑CCGCCGCAATTGCCGCGTGAAATGCTCCAGTCTTGAACGGTAGCAGGCCGCGACCACGGCTGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTG  >  NC_000913/3162864‑3163146

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