Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,148,691	G→A	*357* (TAG→TAA)	plsX →	putative phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,148,691	0	G	A	100.0%	144.4 / NA	45	*357* (TAG→TAA)	plsX	putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/25);  total (20/25)

CGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTA  >  NC_000913/1148565‑1148790                                                                                                                               |                                                                                                    cgcgATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGcc                                                                                         >  2:621350/1‑139 (MQ=255)   cgATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAg                                                                                       >  1:1314014/1‑139 (MQ=255)           cAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCg                                                                               <  1:220417/139‑1 (MQ=255)                      ggCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg                                                                    <  2:445633/139‑1 (MQ=255)                        cGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGt                                                                           >  2:779343/1‑130 (MQ=255)                        cGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGt                                                                           <  1:779343/130‑1 (MQ=255)                         ggTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa                                                                 <  1:621350/139‑1 (MQ=255)                           tGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACg                                                               <  2:1590370/139‑1 (MQ=255)                                          cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGaaaa                                                <  1:1156463/139‑1 (MQ=255)                                               cGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCa                                                                             <  1:923493/105‑1 (MQ=255)                                               cGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCa                                                                             >  2:923493/1‑105 (MQ=255)                                                gAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACt                                          <  1:560623/139‑1 (MQ=255)                                                 aaTTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACt                                                                                 <  2:448192/99‑1 (MQ=255)                                                 aaTTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACt                                                                                 >  1:448192/1‑99 (MQ=255)                                                 aaTTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACt                                          >  2:1232849/1‑138 (MQ=255)                                                 aaTTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACt                                          <  1:1232849/138‑1 (MQ=255)                                                      ccGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGt                                    <  2:1314014/139‑1 (MQ=255)                                                       cGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCg                                                                                     <  1:1136570/89‑1 (MQ=255)                                                       cGCTCGCATGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCg                                                                                     >  2:1136570/1‑89 (MQ=255)                                                           cGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATg                               <  1:128474/139‑1 (MQ=255)                                                                            cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGa                                                                                                <  1:895906/57‑1 (MQ=255)                                                                            cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGa                                                                                                >  2:895906/1‑57 (MQ=255)                                                                            cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGta                                                            <  2:660791/93‑1 (MQ=255)                                                                            cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGta                                                            >  1:660791/1‑93 (MQ=255)                                                                                gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACt          >  1:844115/1‑139 (MQ=255)                                                                                  tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTgg        >  2:1355549/1‑139 (MQ=255)                                                                                        tGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTa  <  2:1699206/139‑1 (MQ=255)                                                                                           gctgctGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGc                                                                                      <  1:108409/52‑1 (MQ=255)                                                                                           gctgctGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGc                                                                                      >  2:108409/1‑52 (MQ=255)                                                                                              gctgGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGaa                                                  <  1:1002723/85‑1 (MQ=255)                                                                                              gctgGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGaa                                                  >  2:1002723/1‑85 (MQ=255)                                                                                               ctgGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTa                                   <  1:1377179/99‑1 (MQ=255)                                                                                               ctgGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTa                                   >  2:1377179/1‑99 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCa                                                                             >  1:70840/1‑47 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCa                                                                             <  2:70840/47‑1 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa                                                                 <  2:1270969/59‑1 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa                                                                 <  1:1194633/59‑1 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa                                                                 >  1:1270969/1‑59 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa                                                                 >  2:1194633/1‑59 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACg                                                               >  2:1264309/1‑61 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACg                                                               <  1:1264309/61‑1 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACt                                          >  2:110403/1‑82 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACt                                          <  1:110403/82‑1 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                     >  2:1659399/1‑87 (MQ=255)                                                                                                         gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                     <  1:1659399/87‑1 (MQ=255)                                                                                                                               |                                                                                                    CGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTA  >  NC_000913/1148565‑1148790

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