Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,697,915	C→T	*212* (TAG→TAA)	pgpC ←	phosphatidylglycerophosphatase C, membrane bound

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,697,915	0	C	T	96.2%	149.2 / ‑1.7	52	*212* (TAG→TAA)	pgpC	phosphatidylglycerophosphatase C, membrane bound
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (30/20);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

CTCATCCCAGGCACGTTTCGCCAGCGTCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTG  >  NC_000913/2697783‑2698041                                                                                                                                     |                                                                                                                               cTCATCCCAGGCACGTTTCGCCAGCGTCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTcc                                                                                                                          <  2:1572640/139‑1 (MQ=255)     tCCCAGGCACGTTTCGCCAGCGTCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGtt                                                                                                                      >  2:829358/1‑139 (MQ=255)                tttCGCCAGCGTCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCa                                                                                                           >  2:155824/1‑139 (MQ=255)                tttCGCCAGCGTCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCa                                                                                                           >  1:291261/1‑139 (MQ=255)                tttCGCCAGCGTCAGCGAGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCa                                                                                                           >  2:1125157/1‑139 (MQ=255)                         cGTCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGTCTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCG‑‑GCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGt                                                                                                <  2:1298397/139‑1 (MQ=255)                           tCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGt                                                                                                >  1:1660751/1‑139 (MQ=255)                             aGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTa                                                                                              <  2:793664/139‑1 (MQ=255)                             aGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTa                                                                                              >  1:1434036/1‑139 (MQ=255)                                                    gTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCt                                                                               >  2:1380082/1‑131 (MQ=255)                                                    gTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCt                                                                               <  1:1380082/131‑1 (MQ=255)                                                    gTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTAATTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGaaa                                                                       >  1:257283/1‑139 (MQ=255)                                                      ggCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATa                                                                     <  2:1434036/139‑1 (MQ=255)                                                       gCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTCTGCTTTACTTTATT‑CAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATg                                                                                 >  2:787995/1‑125 (MQ=255)                                                       gCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATaa                                                                    >  1:295859/1‑139 (MQ=255)                                                       gCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATT‑CAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATg                                                                                 <  1:787995/125‑1 (MQ=255)                                                         tAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAag                                                                  >  2:957872/1‑139 (MQ=255)                                                                gACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGa                                                           <  1:829358/139‑1 (MQ=255)                                                                     cAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTc                                                      <  2:1660751/139‑1 (MQ=255)                                                                     cAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTc                                                      <  2:291261/139‑1 (MQ=255)                                                                                  aaCTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCg                                         >  2:437414/1‑139 (MQ=255)                                                                                       cGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGt                                                                                                                       <  1:5225/56‑1 (MQ=255)                                                                                       cGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGt                                                                                                                       >  2:5225/1‑56 (MQ=255)                                                                                        ggTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAg                                   >  1:1143812/1‑139 (MQ=255)                                                                                         gTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGc                                  >  1:88060/1‑139 (MQ=255)                                                                                         gTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGc                                  >  2:1159357/1‑139 (MQ=255)                                                                                            ataAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCAc                               <  1:155824/139‑1 (MQ=255)                                                                                                   cgcgGCGCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAg                        >  2:1209734/1‑139 (MQ=255)                                                                                                    gcggcgCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCgggg                                                                                                 <  2:1735322/65‑1 (MQ=255)                                                                                                    gcggcgCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCgggg                                                                                                 >  1:1735322/1‑65 (MQ=255)                                                                                                    gcggcgCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGc                       >  1:545654/1‑139 (MQ=255)                                                                                                     cggcgCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACcgcg                                                                                                      <  2:370301/59‑1 (MQ=255)                                                                                                     cggcgCATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACcgcg                                                                                                      >  1:370301/1‑59 (MQ=255)                                                                                                        cgcATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTc                                          >  1:560641/1‑116 (MQ=255)                                                                                                        cgcATTATACACGGACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTc                                          <  2:560641/116‑1 (MQ=255)                                                                                                                    ggACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTc                                                      <  2:129910/92‑1 (MQ=255)                                                                                                                    ggACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTc                                                      >  1:129910/1‑92 (MQ=255)                                                                                                                     gACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTATTg                                            >  2:753750/1‑101 (MQ=255)                                                                                                                     gACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTg                                            <  1:753750/101‑1 (MQ=255)                                                                                                                     gACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGt                                           >  2:878765/1‑102 (MQ=255)                                                                                                                     gACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGt                                           <  1:878765/102‑1 (MQ=255)                                                                                                                     gACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGAt      >  2:1166545/1‑139 (MQ=255)                                                                                                                     gACGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGAt      >  2:389304/1‑139 (MQ=255)                                                                                                                     gACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCa                                <  1:879048/113‑1 (MQ=255)                                                                                                                     gACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCa                                >  2:879048/1‑113 (MQ=255)                                                                                                                      aCGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGc                                                                                         <  1:280194/55‑1 (MQ=255)                                                                                                                      aCGCTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGc                                                                                         >  2:280194/1‑55 (MQ=255)                                                                                                                         cTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTg  >  2:467323/1‑139 (MQ=255)                                                                                                                         cTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTg  >  1:462730/1‑139 (MQ=255)                                                                                                                         cTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTg  >  1:1396039/1‑139 (MQ=255)                                                                                                                         cTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTg  <  2:88060/139‑1 (MQ=255)                                                                                                                         cTATGCTTTACTTTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTg  <  2:1396039/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                               CTCATCCCAGGCACGTTTCGCCAGCGTCAGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTG  >  NC_000913/2697783‑2698041

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