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	98.6%	215.9 / NA	71	*173* (TAG→TAA)	pgpA	phosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (34/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCA  >  NC_000913/436973‑437240                                                                                                                                       |                                                                                                                                      tGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgt                                                                                                                                   >  2:1010162/1‑139 (MQ=255) tGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgt                                                                                                                                   >  1:199614/1‑139 (MQ=255)  ggCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgtg                                                                                                                                  >  1:76220/1‑139 (MQ=255)  ggCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgtg                                                                                                                                  >  1:1388545/1‑139 (MQ=255)  ggCCGATCCGCTGGTTTGATCACAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCATGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgtg                                                                                                                                  >  2:995960/1‑139 (MQ=255)    ccGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                                <  1:995960/139‑1 (MQ=255)    ccGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                                <  1:1554650/139‑1 (MQ=255)     cGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAc                                                                                                                               <  1:342774/139‑1 (MQ=255)      gATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                              <  2:770880/139‑1 (MQ=255)                         aaTGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCtt                                                                                                           >  1:1551044/1‑139 (MQ=255)                         aaTGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCtt                                                                                                           >  1:785090/1‑139 (MQ=255)                          aTGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTg                                                                                                          >  1:1133057/1‑139 (MQ=255)                              gCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCAc                                                                                                      <  1:299062/139‑1 (MQ=255)                               cATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACg                                                                                                     >  1:234764/1‑139 (MQ=255)                                aTGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                    <  1:464968/139‑1 (MQ=255)                                  ggcggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGAAGGCATACTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCAccgct                                                                                                  >  2:1407528/1‑135 (MQ=255)                                   gcggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGtctc                                                                                                 >  2:875678/1‑139 (MQ=255)                                     ggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:1010162/139‑1 (MQ=255)                                     ggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  2:785090/139‑1 (MQ=255)                                       cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                    <  2:1123542/132‑1 (MQ=255)                                       cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                    >  1:1123542/1‑132 (MQ=255)                                       cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcatca                                                                                             <  2:101368/139‑1 (MQ=255)                                                  tGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                                  <  2:234764/139‑1 (MQ=255)                                                        acgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGa                                                                                     <  1:1027606/130‑1 (MQ=255)                                                        acgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGa                                                                                     >  2:1027606/1‑130 (MQ=255)                                                           atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTaaa                                                                         <  1:1590683/139‑1 (MQ=255)                                                             atTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATc                                                                       <  2:1388545/139‑1 (MQ=255)                                                                   gCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:1482797/1‑109 (MQ=255)                                                                   gCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:1482797/109‑1 (MQ=255)                                                                    ccGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGAc                                                                <  2:76220/139‑1 (MQ=255)                                                                                   cAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:957733/93‑1 (MQ=255)                                                                                   cAGGCATCCTGTATTTTATCGGTCATCACAGGCCGCTGGGTATTCCGTCGTATTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:957733/1‑93 (MQ=255)                                                                                        aTCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTaaa                                                                         <  2:1029149/110‑1 (MQ=255)                                                                                        aTCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGAGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTaaa                                                                         >  1:1029149/1‑110 (MQ=255)                                                                                                  ttATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:125385/1‑78 (MQ=255)                                                                                                  ttATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:125385/78‑1 (MQ=255)                                                                                                   tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTa                                                                           <  2:1459435/97‑1 (MQ=255)                                                                                                   tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTa                                                                           >  1:1459435/1‑97 (MQ=255)                                                                                                   tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCg                                 <  1:86513/139‑1 (MQ=255)                                                                                                     tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                                <  2:645198/42‑1 (MQ=255)                                                                                                     tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                                >  1:645198/1‑42 (MQ=255)                                                                                                     tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                                           >  1:768285/1‑111 (MQ=255)                                                                                                     tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                                           <  2:768285/111‑1 (MQ=255)                                                                                                     tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACAGCTGCCTGATGTTACGCTTTTCAGGTCTCATCAGGCCTGGACTCTTATTTAAATCCTGCGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGc                               >  2:598052/1‑139 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                              >  1:1342055/1‑41 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                              <  2:1342055/41‑1 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                              >  2:23963/1‑41 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                              <  1:23963/41‑1 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:611091/1‑72 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:611091/72‑1 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                       <  2:437452/80‑1 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTg                                                                                       >  1:437452/1‑80 (MQ=255)                                                                                                        gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAAt                            >  1:292040/1‑139 (MQ=255)                                                                                                        gTCATCACTGGCCGCT‑GGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         >  2:594555/1‑61 (MQ=38)                                                                                                        gTCATCACTGGCCGCT‑GGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         <  1:594555/61‑1 (MQ=37)                                                                                                                  gCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCg                                                        >  1:506145/1‑101 (MQ=255)                                                                                                                  gCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCg                                                        <  2:506145/101‑1 (MQ=255)                                                                                                                  gCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGc                  >  2:547420/1‑139 (MQ=255)                                                                                                                   ccGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                                  <  2:17441/74‑1 (MQ=255)                                                                                                                   ccGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGActc                                                                                  >  1:17441/1‑74 (MQ=255)                                                                                                                   ccGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCt                 >  1:826024/1‑139 (MQ=255)                                                                                                                      cTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCa              <  1:875678/139‑1 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGTCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATcc                                                                      <  1:978322/82‑1 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  2:102692/1‑57 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  1:102692/57‑1 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATca                                                                      >  2:978322/1‑81 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCa                                      >  2:1307954/1‑114 (MQ=255)                                                                                                                       tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCa                                      <  1:1307954/114‑1 (MQ=255)                                                                                                                                  tCGTAATTGTGCACCGATGCCTTATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCa  <  1:547420/139‑1 (MQ=255)                                                                                                                                   cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACg                                               <  2:537276/93‑1 (MQ=255)                                                                                                                                   cGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACg                                               >  1:537276/1‑93 (MQ=255)                                                                                                                                       |                                                                                                                                      TGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCA  >  NC_000913/436973‑437240

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