Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,722,946	G→A	A74T (GCT→ACT)	pssA →	phosphatidylserine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,722,946	0	G	A	96.4%	73.5 / ‑4.2	29	A74T (GCT→ACT)	pssA	phosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/16);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

AACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGC  >  NC_000913/2722705‑2723188                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   aaCAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGc                                                                                                                                                                                                                                           <  2:134867/251‑1 (MQ=255)         agaaaagaaCTGTGATGGTGCAAAAAATAAAGCGAAATAAACATCAACAAAACCGTGCCCAAATGCCCAAGATGTCTCAATAAGTTGATGATGGGGATGTCTTTTACGAGACCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGAGAAGCAGCGCATTTGCAGTGTCGCCCTGTATCGAGAACAGGATGACGGTGGCAAAGGCAGTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCg                                                                                                                                                                                                                                   <  2:125262/243‑1 (MQ=255)               gCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACt                                                                                                                                                                                                                              <  1:315802/250‑1 (MQ=255)                                   tttAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGCCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGt                                                                                                                                                                                                            >  1:316374/1‑248 (MQ=255)                                   tttAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAACGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGttg                                                                                                                                                                                                          >  2:172634/1‑248 (MQ=255)                                              aaaaCATCAACAACCCCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTg                                                                                                                                                                                              <  1:245545/250‑1 (MQ=255)                                                                tGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTAGCTTTCTTTAACGCTCCGGACGAATTCCGGGAGACGCTGCTGGAAAAAATAACCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGAGTGACGGTGGCAAAGGCATTCTGAACGCGATGAATGAGACTAAAAGGCAGCGTCCGGAACTGGCTGTGCGGGGGATGGTCGACTGGCATGGTGCACaagggggacacgttg                                                                                                                                                                            >  2:305818/1‑238 (MQ=255)                                                                     aacgaCCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGgcgcg                                                                                                                                                                       <  1:151071/247‑2 (MQ=255)                                                                        tACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCg                                                                                                                                                                    >  1:143770/1‑251 (MQ=255)                                                                                                 gtatgcctattccTTTTCCGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGt                                                                                                                                           <  2:55216/239‑1 (MQ=255)                                                                                                           ttCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTCTCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCAt                                                                                                                                 >  1:230601/1‑251 (MQ=255)                                                                                                                ttACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGgcg                                                                                                                             <  1:174941/250‑1 (MQ=255)                                                                                                                                          cGCTGCTGGAAAAAATAGCCAGCGAGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGGGCTGGTAGACTGGCATCGTGCACAACGGGGACGCATTGGCGCTGCGGCATCGAACACTACAGCTGACTGGTACTGCCGCGTGGCGCAGGAAAATCCGGGCGGAGATGttagg                                                                                                  >  2:202309/1‑248 (MQ=255)                                                                                                                                                     aaaaTATCCAGCGCTAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCGAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGt                                                                                       <  1:16547/251‑1 (MQ=255)                                                                                                                                                                             aTTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGATGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGATGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                               >  2:308651/1‑251 (MQ=255)                                                                                                                                                                                     tGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCGTGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTc                                                       >  1:3343/1‑251 (MQ=255)                                                                                                                                                                                        cGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGCTGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCAGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGc                                                    >  1:133693/1‑251 (MQ=255)                                                                                                                                                                                           ccTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACt                                                 <  1:202309/251‑1 (MQ=255)                                                                                                                                                                                              gTATCTCTAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCGGCGTGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACttt                                               <  2:3343/250‑1 (MQ=255)                                                                                                                                                                                                    tgAACAGGATGACGGTGGCAAAGGCATTCTGAACGCTTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAgg                                          <  1:263386/248‑1 (MQ=255)                                                                                                                                                                                                           gATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTatcat                                 <  1:308651/251‑1 (MQ=255)                                                                                                                                                                                                                    ggCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGCTGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTatcattgatcgga                         >  2:314724/1‑242 (MQ=255)                                                                                                                                                                                                                    ggCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTatcat                                 <  1:314724/242‑1 (MQ=255)                                                                                                                                                                                                                         aGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTAc                   >  2:135556/1‑251 (MQ=255)                                                                                                                                                                                                                              ttCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTtat              >  1:16589/1‑251 (MQ=255)                                                                                                                                                                                                                              ttCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGAGCTGGACGACTGGCCTCGTGCACAACGTGGACGACTTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAACGGCTTTATCATCGACGATAGCGTACTTtat              >  1:367149/1‑251 (MQ=255)                                                                                                                                                                                                                                  taaCGGGTTGTATGGGACTAAAAGGCAGCGTCCGGAACTGGATGTTCGGGTGCTGGTCGACTTGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGc           <  2:79332/249‑1 (MQ=255)                                                                                                                                                                                                                                     cgcgTTGTATGAGACTAAAAGGCAGCGTCCGGAACTTGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTg       <  1:238452/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGc  <  1:227988/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   AACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGC  >  NC_000913/2722705‑2723188

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