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

CCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTT  >  NC_000913/436976‑437226                                                                                                                                    |                                                                                                                        ccGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                  <  1:846683/139‑1 (MQ=255)                        tgtgCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                           >  1:565270/1‑139 (MQ=255)                                    cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                  >  2:104843/1‑104 (MQ=255)                                    cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                  <  1:104843/104‑1 (MQ=255)                                              aTGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAct                                                                     <  2:565270/139‑1 (MQ=255)                                                 aTCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCtt                                                                  >  1:795430/1‑139 (MQ=255)                                                  tCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTa                                                                 >  1:1412047/1‑139 (MQ=255)                                                   cgacgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTAt                                                                >  1:1024728/1‑139 (MQ=255)                                                        atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         >  2:241559/1‑125 (MQ=255)                                                        atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         <  1:241559/125‑1 (MQ=255)                                                        atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTaaa                                                           >  1:1352507/1‑139 (MQ=255)                                                                gCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                           >  1:1150007/1‑99 (MQ=255)                                                                gCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                           <  2:1150007/99‑1 (MQ=255)                                                                 ccGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         <  1:656054/116‑1 (MQ=255)                                                                 ccGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         >  2:656054/1‑116 (MQ=255)                                                                    gggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACg                                                                                                >  1:344547/1‑90 (MQ=255)                                                                    gggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACg                                                                                                <  2:344547/90‑1 (MQ=255)                                                                     ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGa                                                   >  2:699364/1‑134 (MQ=255)                                                                     ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGa                                                   <  1:699364/134‑1 (MQ=255)                                                                         aTTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCg                                          <  2:1412047/139‑1 (MQ=255)                                                                              cGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGc                                                                                                           >  1:1285834/1‑69 (MQ=255)                                                                              cGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGc                                                                                                           <  2:1285834/69‑1 (MQ=255)                                                                                      tCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGt                             >  2:1219579/1‑139 (MQ=255)                                                                                             ttttATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                    >  1:1037976/1‑93 (MQ=255)                                                                                             ttttATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                    <  2:1037976/93‑1 (MQ=255)                                                                                                tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcatca                                                                               <  2:6905/79‑1 (MQ=255)                                                                                                tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcatca                                                                               >  1:6905/1‑79 (MQ=255)                                                                                                tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCg                   >  2:1250758/1‑139 (MQ=255)                                                                                                 aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCgg                  <  1:1219579/139‑1 (MQ=255)                                                                                                  tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGc                 >  2:898850/1‑139 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                <  1:796761/41‑1 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                >  2:796761/1‑41 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         >  1:29706/1‑80 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         >  2:1339807/1‑80 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         <  1:1339807/80‑1 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         <  2:29706/80‑1 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActct                                                                   <  2:278941/86‑1 (MQ=255)                                                                                                     gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAAt              >  2:648451/1‑139 (MQ=255)                                                                                                     atcatcaCTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActct                                                                   >  1:278941/2‑86 (MQ=255)                                                                                                                 cGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTTGGCAATCTGTTCCGGCtt  <  2:795430/139‑1 (MQ=255)                                                                                                                    tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtg                                                                                                     >  2:1025243/1‑37 (MQ=38)                                                                                                                    tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtg                                                                                                     >  2:1372675/1‑37 (MQ=38)                                                                                                                    tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtg                                                                                                     <  1:1372675/37‑1 (MQ=255)                                                                                                                    tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGAtg                                                                                                     <  1:1025243/37‑1 (MQ=38)                                                                                                                    tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAct                                                                     <  1:763506/69‑1 (MQ=255)                                                                                                                    tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAct                                                                     >  2:763506/1‑69 (MQ=255)                                                                                                                           tCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                           >  2:585574/1‑40 (MQ=25)                                                                                                                           tCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                           <  1:585574/40‑1 (MQ=25)                                                                                                                              gTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAc                                                                      >  2:630920/1‑58 (MQ=255)                                                                                                                              gTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAc                                                                      <  1:630920/58‑1 (MQ=255)                                                                                                                                cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         <  2:563963/53‑1 (MQ=255)                                                                                                                                cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                         >  1:563963/1‑53 (MQ=255)                                                                                                                                cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                             <  2:11377/81‑1 (MQ=255)                                                                                                                                cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                             >  1:11377/1‑81 (MQ=255)                                                                                                                                    |                                                                                                                        CCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTT  >  NC_000913/436976‑437226

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