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%	166.9 / NA	51	*173* (TAG→TAA)	pgpA	phosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/24);  total (27/24)

TGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCAC  >  NC_000913/436973‑437241                                                                                                                                       |                                                                                                                                       tGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgt                                                                                                                                    >  1:1110534/1‑139 (MQ=255)            cTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGcc                                                                                                                         <  1:389252/139‑1 (MQ=255)                  tGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtgt                                                                                                                   >  2:343352/1‑139 (MQ=255)                     tCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGAc                                                                                                                >  2:782834/1‑139 (MQ=255)                         aaTGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGcc                                                                                                                         >  2:325990/1‑126 (MQ=255)                         aaTGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGcc                                                                                                                         <  1:325990/126‑1 (MQ=255)                                aTGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACg                                                                                                               <  2:261776/129‑1 (MQ=255)                                aTGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACg                                                                                                               >  1:261776/1‑129 (MQ=255)                                     ggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                                <  1:343352/139‑1 (MQ=255)                                              aTCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGa                                                                                                                 <  2:424088/113‑1 (MQ=255)                                              aTCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGa                                                                                                                 >  1:424088/1‑113 (MQ=255)                                                      cgacgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGAt                                                                                                                            <  2:1095938/94‑1 (MQ=255)                                                      cgacgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGAt                                                                                                                            >  1:1095938/1‑94 (MQ=255)                                                                      ggggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                     <  1:297427/101‑1 (MQ=255)                                                                      ggggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                     >  2:297427/1‑101 (MQ=255)                                                                      ggggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGAc                                                                 >  2:641402/1‑137 (MQ=255)                                                                      ggggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGAc                                                                 <  1:641402/137‑1 (MQ=255)                                                                        ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                                <  1:1505587/104‑1 (MQ=255)                                                                        ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                                >  2:1505587/1‑104 (MQ=255)                                                                        ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACa                                                                <  2:654592/136‑1 (MQ=255)                                                                        ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACa                                                                >  1:654592/1‑136 (MQ=255)                                                                              ttCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGt                                                       >  2:1179358/1‑139 (MQ=255)                                                                                  gCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTtat                                                   >  1:109930/1‑139 (MQ=255)                                                                                  gCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTtat                                                   <  2:1110534/139‑1 (MQ=255)                                                                                    aGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATAc                                                 <  2:109930/139‑1 (MQ=255)                                                                                         tCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcatca                                                                                              >  2:509977/1‑89 (MQ=255)                                                                                         tCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcatca                                                                                              <  1:509977/89‑1 (MQ=255)                                                                                                   tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCg                                  >  1:431775/1‑139 (MQ=255)                                                                                                   tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCg                                  >  2:923702/1‑139 (MQ=255)                                                                                                    aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                        <  2:563586/84‑1 (MQ=255)                                                                                                    aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                        >  1:563586/1‑84 (MQ=255)                                                                                                    aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCgg                                 <  1:923702/139‑1 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCt                                                                                                                        >  2:494814/1‑48 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCt                                                                                                                        <  1:494814/48‑1 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                        <  2:720679/80‑1 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                        >  1:720679/1‑80 (MQ=255)                                                                                                           atcaCTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtgt                                                                                                                   <  1:437827/50‑1 (MQ=255)                                                                                                           atcaCTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtgt                                                                                                                   >  2:437827/1‑50 (MQ=255)                                                                                                                  gCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACg                                                                                                               >  2:1035418/1‑47 (MQ=255)                                                                                                                  gCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACg                                                                                                               <  1:1035418/47‑1 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                        >  2:853882/1‑65 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                        <  1:853882/65‑1 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                                   <  2:926450/70‑1 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                                   >  1:926450/1‑70 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCaa              >  1:1152990/1‑139 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGCGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                                            <  2:571670/93‑1 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGCGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                                            >  1:571670/1‑93 (MQ=255)                                                                                                                        gggTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCaaa             <  1:782834/139‑1 (MQ=255)                                                                                                                                   cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                                            >  1:263119/1‑81 (MQ=255)                                                                                                                                   cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                                            <  2:263119/81‑1 (MQ=255)                                                                                                                                   cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCAc  >  1:860104/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       TGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCAC  >  NC_000913/436973‑437241

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