Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	437,107	G→A	*173* (TAG→TAA)	pgpA →	phosphatidylglycerophosphatase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	437,107	0	G	A	100.0%	174.3 / NA	55	*173* (TAG→TAA)	pgpA	phosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/30);  total (25/30)

GGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCA  >  NC_000913/436974‑437240                                                                                                                                      |                                                                                                                                      ggCCGAACCGCTGGTCAGATCGCAATGTGCATGGCGGCATGGGGATGATGATCGTCGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgtg                                                                                                                                  <  2:467951/139‑1 (MQ=255)        ccGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGAt                                                                                                                           >  2:414427/1‑139 (MQ=255)                 tGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtgt                                                                                                                  <  1:292818/139‑1 (MQ=255)                 tGAGCGAAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATTCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtgt                                                                                                                  <  1:677975/139‑1 (MQ=255)                   aTCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGa                                                                                                                >  2:997570/1‑139 (MQ=255)                    tCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGAc                                                                                                               <  1:670237/139‑1 (MQ=255)                          tgtgCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         <  1:997570/139‑1 (MQ=255)                           gtgCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTc                                                                                                        <  1:540033/139‑1 (MQ=255)                                    tgcATGGGGATCATGTTCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGAGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  2:1008202/138‑1 (MQ=255)                                    ggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  1:970340/1‑139 (MQ=255)                                    ggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:858865/139‑1 (MQ=255)                                      cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcatca                                                                                             >  2:232013/1‑139 (MQ=255)                                      cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:1019700/137‑1 (MQ=255)                                      cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:1019700/1‑137 (MQ=255)                                      cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGACtcatca                                                                                             <  1:414427/139‑1 (MQ=255)                                                    tCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:455621/1‑123 (MQ=255)                                                    tCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:455621/123‑1 (MQ=255)                                                       acgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtgt                                                                                                                  <  1:254058/101‑1 (MQ=255)                                                       acgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtgt                                                                                                                  >  2:254058/1‑101 (MQ=255)                                                          atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         <  2:400471/107‑1 (MQ=255)                                                          atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         >  1:400471/1‑107 (MQ=255)                                                          atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTaaa                                                                         >  1:1306307/1‑139 (MQ=255)                                                           tatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  1:1303813/1‑116 (MQ=255)                                                           tatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  2:1303813/116‑1 (MQ=255)                                                                  gCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:1202742/1‑109 (MQ=255)                                                                  gCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:1202742/109‑1 (MQ=255)                                                                    cGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACa                                                               <  2:970340/139‑1 (MQ=255)                                                                     ggggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                       >  1:415616/1‑114 (MQ=255)                                                                     ggggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                       <  2:415616/114‑1 (MQ=255)                                                                          gATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGc                                                         <  1:316612/139‑1 (MQ=255)                                                                                         ccTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  2:851522/86‑1 (MQ=255)                                                                                         ccTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  1:851522/1‑86 (MQ=255)                                                                                         ccTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGtt                                          <  2:1306307/139‑1 (MQ=255)                                                                                            gTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTAttt                                                                            <  2:1032206/102‑1 (MQ=255)                                                                                            gTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTAttt                                                                            >  1:1032206/1‑102 (MQ=255)                                                                                               ttttATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGt                                    <  1:905736/139‑1 (MQ=255)                                                                                                   aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCgg                                <  1:543052/139‑1 (MQ=255)                                                                                                   aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCgg                                >  2:543052/1‑139 (MQ=255)                                                                                                   aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCgg                                >  1:1277926/1‑139 (MQ=255)                                                                                                    tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGc                               >  2:109572/1‑139 (MQ=255)                                                                                                    tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGc                               >  2:215679/1‑139 (MQ=255)                                                                                                       gTCATCACTGGCCGCTGGGTGTTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGtt                                          >  2:723123/1‑125 (MQ=255)                                                                                                       gTCATCACTGGCCGCTGGGTGTTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGtt                                          <  1:723123/125‑1 (MQ=255)                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAct                                                                                   >  2:859119/1‑84 (MQ=255)                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAct                                                                                   <  1:859119/84‑1 (MQ=255)                                                                                                                 gCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGc                  >  2:212851/1‑139 (MQ=255)                                                                                                                      tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACg                                                                                                     >  2:765589/1‑51 (MQ=255)                                                                                                                      tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACg                                                                                                     <  1:765589/51‑1 (MQ=255)                                                                                                                      tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCaa             >  2:390548/1‑139 (MQ=255)                                                                                                                           aTTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         >  2:426753/1‑42 (MQ=255)                                                                                                                           aTTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         <  1:426753/42‑1 (MQ=255)                                                                                                                                gTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                                  <  2:1008140/60‑1 (MQ=255)                                                                                                                                gTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                                  >  1:1008140/1‑60 (MQ=255)                                                                                                                                 tCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCa  <  2:1277926/139‑1 (MQ=255)                                                                                                                                 tCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCa  <  1:232013/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      GGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCA  >  NC_000913/436974‑437240

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