Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,500,060	C→T	pseudogene (475/942 nt)	yjgX ←	pseudogene fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,500,060	0	C	T	96.5%	83.0 / ‑3.9	30	pseudogene (475/942 nt)	yjgX	pseudogene fragment
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (15/13);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCC  >  NC_000913/4499924‑4500186                                                                                                                                         |                                                                                                                               aaaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtac                                                                                                                              <  1:190762/139‑1 (MQ=255) aaaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtac                                                                                                                              <  2:747838/139‑1 (MQ=255)  aaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                             <  1:454953/139‑1 (MQ=255)        gTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAtt                                                                                                                       <  1:413648/139‑1 (MQ=255)                   tgAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTCGGATGAGAAc                                                                                                            <  2:382224/138‑1 (MQ=255)                        gATACCATCAAGCTGGTCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCcatc                                                                                                       <  2:396641/139‑1 (MQ=255)                                cAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAg                                                                                               >  1:370856/1‑139 (MQ=255)                                  aGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                               >  1:631894/1‑121 (MQ=255)                                  aGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                               <  2:631894/121‑1 (MQ=255)                                                    tCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                           <  2:370856/139‑1 (MQ=255)                                                    tCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                           <  1:816557/139‑1 (MQ=255)                                                                   cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCtt                                                                                   >  1:104655/1‑116 (MQ=255)                                                                   cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCtt                                                                                   <  2:104655/116‑1 (MQ=255)                                                                   cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTACTCGGATAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcg                                                            >  1:317486/1‑139 (MQ=255)                                                                    tGGCGAGGTAACAGGCAATTTCTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTCATACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                           <  2:647767/139‑1 (MQ=255)                                                                         aGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGtt                                                      <  2:619772/139‑1 (MQ=255)                                                                           gTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCgg                                                                                <  1:758351/111‑1 (MQ=255)                                                                           gTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCgg                                                                                >  2:758351/1‑111 (MQ=255)                                                                           gTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCt                                                    >  2:23977/1‑139 (MQ=255)                                                                                     aTTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTa                                          >  2:380368/1‑139 (MQ=255)                                                                                      tttttCGGGGATACTGCTCCAGGTAATTACACTGTGAGGGGTTAAGGCTGTGACAGGGATGTAGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGGACGCTCTGTTTCCTTTaa                                         <  2:317486/139‑1 (MQ=255)                                                                                                 tACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACATGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGt                              <  2:910552/139‑1 (MQ=255)                                                                                                      cTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCt                         >  2:44968/1‑139 (MQ=255)                                                                                                      cTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATT‑TCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTg                        >  2:368535/1‑139 (MQ=255)                                                                                                         cAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGtt                      >  1:858298/1‑139 (MQ=255)                                                                                                          aGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGttt                     >  2:456228/1‑139 (MQ=255)                                                                                                                   ttCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaag            <  1:380368/139‑1 (MQ=255)                                                                                                                      ggttAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaagcga         >  2:157165/4‑139 (MQ=255)                                                                                                                             aGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaagcgaag       <  1:983516/134‑1 (MQ=255)                                                                                                                             aGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaagcgaag       >  2:983516/1‑134 (MQ=255)                                                                                                                             aGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGcc  >  2:380210/1‑139 (MQ=255)                                                                                                                                     cTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCAt                           <  1:1014477/106‑1 (MQ=255)                                                                                                                                     cTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCAt                           >  2:1014477/1‑106 (MQ=255)                                                                                                                                         |                                                                                                                               AAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCC  >  NC_000913/4499924‑4500186

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