Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,070,891	G→T	12.2%	pseudogene (106/426 nt)	yeeP →	CP4‑44 prophage; putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,070,891	0	G	T	12.2%	107.0 / 4.5	33	pseudogene (106/426 nt)	yeeP	CP4‑44 prophage; putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base T (2/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CTCTCGCGTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCGCCAGTCCCCTGATGACCCGACTGCAGGACGAGCTGCGCACGGAATCTGTCCGGTCTCAGGCCCGTGAACAGTTTA  >  NC_000913/2070644‑2071077                                                                                                                                                                                                                                                        |                                                                                                                                                                                           actcggcgTAGGCGGGTCCCGTGCCATTTGGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGGGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACAt                                                                                                                                                                  <  2:128920/269‑1 (MQ=255)      gcgTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTCTTCAGCCTTCCCCTGCCCCGGCCCAGAAAATTCGCGAAAACCAGGGAGAGGTAATCCGTATGTTCCGGGCCGTACCACCggg                                                                                                                                                             >  2:320957/1‑273 (MQ=255)      gcgTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCCGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGGTATGGGCGGGATTCAGCCTTCTCCTGCACAGGCACAGAAAATTTGCGAAAAAACGGAGGGGGTATTCCGTCTGTTCCGGCCCGTa                                                                                                                                                                     >  1:128920/1‑266 (MQ=255)                                        tGACCTTGTACTGGGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGGGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTg                                                                                                                           <  1:320957/274‑1 (MQ=255)                                                      ggCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAg                                                                                                             <  1:29641/274‑1 (MQ=255)                                                                             gCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGTCAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGTTGTTCCGTCTGTTCCGGc                                                                                                                                                                          >  1:184893/1‑190 (MQ=255)                                                                             gCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGTTAtt                                                                                                                                                                                        >  2:226487/1‑176 (MQ=255)                                                                                ctgtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGTTGTTCCGTCTGTTCCGGCCCg                                                                                                                                                                       <  2:184893/190‑1 (MQ=255)                                                                                ctgtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGTTATTCCg                                                                                                                                                                                     <  1:226487/176‑1 (MQ=255)                                                                                              gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTcc                                                                                                                                                                             >  1:230433/1‑170 (MQ=255)                                                                                              gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTTCGCTCATGCCAGCGCTTCCCGACCATGCCGCCacggc                                                                      >  2:44015/1‑269 (MQ=255)                                                                                              gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCACCAGTccc                                                                     >  1:190074/1‑274 (MQ=255)                                                                                                 tttCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTa                                                                                                                                                                                          <  2:103113/154‑1 (MQ=255)                                                                                                 tttCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGc                                                                                                                                                                          <  2:230433/170‑1 (MQ=255)                                                                                                                          ggACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGt                                                                                                                >  2:252679/1‑203 (MQ=255)                                                                                                                             cATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAg                                                                                                             <  1:252679/203‑1 (MQ=255)                                                                                                                                   cagGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCGCCAGTCCCCTGATGACCCGACTGCAGGACGAGCTGCGCACGGaa                                 <  1:44015/273‑1 (MQ=255)                                                                                                                                         cTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTg                                                                                                                                  >  1:137279/1‑170 (MQ=255)                                                                                                                                         cTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGa                                                                                                                         >  1:224612/1‑179 (MQ=255)                                                                                                                                            tttGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGa                                                                                                                               <  2:137279/170‑1 (MQ=255)                                                                                                                                            tttGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATAc                                                                                                                      <  2:224612/179‑1 (MQ=255)                                                                                                                                                                cAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGt                                                                                                                >  1:119484/1‑165 (MQ=255)                                                                                                                                                                   aaCGGAGCCCGGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAg                                                                                                             <  2:119484/165‑1 (MQ=255)                                                                                                                                                                                cATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGccc                                                                                                                                            >  2:41569/1‑121 (MQ=255)                                                                                                                                                                                   gAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCa                                                                                                                                         <  1:41569/121‑1 (MQ=255)                                                                                                                                                                                                     aTTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGt                                                                                                                >  2:358954/1‑128 (MQ=255)                                                                                                                                                                                                        cAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAg                                                                                                             <  1:358954/128‑1 (MQ=255)                                                                                                                                                                                                                           cacagAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATg                                                                                                   >  1:258131/1‑119 (MQ=255)                                                                                                                                                                                                                              agAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACa                                                                                                <  2:258131/119‑1 (MQ=255)                                                                                                                                                                                                                                     tCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCGCCAGTCCCCTGATGACCCGACTGCAGGACGa                                              >  2:342657/1‑162 (MQ=255)                                                                                                                                                                                                                                      cgcgAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCGCCAGTCCCCTGATGACCCGACTGCAGGACGAGCTGCGCACGGAATCTGTCCGGTCTCAGGCCCGTGAACAGt     >  2:239569/1‑202 (MQ=255)                                                                                                                                                                                                                                        cgAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCGCCAGTCCCCTGATGACCCGACTGCAGGACGAGCt                                           <  1:342657/162‑1 (MQ=255)                                                                                                                                                                                                                                         gAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCGCCAGTCCCCTGATGACCCGACTGCAGGACGAGCTGCGCACGGAATCTGTCCGGTCTCAGGCCCGTGAACAGTTTa  <  1:239569/202‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                           CTCTCGCGTAAGCGGGTACCGTGACATTCTGCCTGAACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATATTCGCGAAAAAACGGAGGCGGTATTCCGTCTGTTCCGGCCCGTACATCCGGTTGTGGCCGTATCGGCCCGCACCGGCTGGGAACTGGATACGCTGGTCAGTGCGCTCATGACAGCGCTTCCCGACCATGCCGCCAGTCCCCTGATGACCCGACTGCAGGACGAGCTGCGCACGGAATCTGTCCGGTCTCAGGCCCGTGAACAGTTTA  >  NC_000913/2070644‑2071077

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