Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,310,841	C→T	E110K (GAA→AAA)	pnp ←	polynucleotide phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,310,841	0	C	T	97.5%	108.0 / NA	41	E110K (GAA→AAA)	pnp	polynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (22/17);  total (24/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 = 1.00e+00

ACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCC  >  NC_000913/3310708‑3310977                                                                                                                                      |                                                                                                                                         aCGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGaa                                                                                                                                 >  1:431082/1‑143 (MQ=255)   gcgGGCAGCCCCAATCGCGCCATTGAACGGAATACCCGACAGAGCCAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGc                                                                                                                               <  1:873637/143‑1 (MQ=255)   gcgGACAGCGTCAATCGGTCCATTGAATGTAATACCACATAGAGACAGCGCTGCGGACGCACCAATCGGCGCGACGGTAGCCGGGTTAACGTGCGGGTTAACAGAAACCACGTTGGCGATAACCTGAACTTTGTTGACGAAGc                                                                                                                               <  2:267950/143‑1 (MQ=255)       gaagcaCCAAGCGGGCCATTGAACGGAATACGAGACAGAGACAGCGCTGCGGATGCTCCAATCAAAGCGACGGTATCAGGGGTAACATGCGGGTTAACAGAAATCACGGTGGCGATAACCTGCGCGACGCTGACGAAGCCtac                                                                                                                           >  2:244385/3‑141 (MQ=255)              cAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAAc                                                                                                                    >  1:92650/1‑143 (MQ=255)                    ggCGATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCggg                                                                                                              <  1:679050/143‑1 (MQ=255)                    ggCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCggg                                                                                                              >  1:473299/1‑143 (MQ=255)                      ccATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGACCAGCGGGCg                                                                                                            >  2:378217/1‑143 (MQ=255)                        aTTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGATAACTGGCGGGATAACAGAAACCACGATGGCGAAAACCTGAACTTTGTTGACGAGGCCTTCCAGGAACAGCGGGCgga                                                                                                          >  2:645581/1‑141 (MQ=255)                             gcGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGAAGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCggg                                                                                                     <  1:803949/142‑1 (MQ=255)                                 aatgcgggacagagggagGGCTGCGGCAGGACTCAGGAGCGCGACGATATGCGGGTTAACTTGCGGGTTTACAGAAACTACGGTGGCGGTAACCTGAACTGTGTTGAGGAGGGCTTCCGGGAACAGCGGGCGAATCGGGCGGt                                                                                                 <  2:620264/136‑1 (MQ=255)                                 aaTACCAGGCAGAGGCAGCGCTGCGGACGCCCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTACCAGAAACCGCGGTGGCGATAACCTGAACTTTGTGGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGt                                                                                                 <  2:600102/143‑1 (MQ=255)                                      cagacagaGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGGTAACCTGGACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGCCaagc                                                                                            >  2:369149/1‑141 (MQ=255)                                        gacagagacagCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGACGCCTTCCGGGAACAGCGGGCGAATCGGGCGGtcagtcag                                                                                          >  1:510404/1‑143 (MQ=255)                                         acagagacagCGCTGCGGAAGCACCACTCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTCGCGATAACCTGAACTAAGCTGACGAAGCCTTCCGGCAACAGCGGGCGCACCGGGCGATCAATCAGa                                                                                         >  2:840483/1‑143 (MQ=255)                                             agacagCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACAAGAGGGTAAAGAGAACGCCCGGTGCCGATAACCTGAACTTTGTTGACGAGGCGTTCAGGGAACGGCGGGCGGAACGGGCGGGCAATCAGAcgcg                                                                                     >  1:211552/1‑143 (MQ=255)                                              gtccgcgcTGCGAAAGCGCCAACCGTCGCCACGATATCCGGGTTAACTTGCGGGTCAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCGTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGAcgcgc                                                                                    <  2:586516/139‑1 (MQ=255)                                                    gcTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAAATTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTATGTTGACCAAGCCTTCCGGGAACAGCGGGCGAATCGGCAGGGAGATCAGACGCGCGagcgc                                                                              >  2:227263/1‑139 (MQ=255)                                                     cttcGGAAGCGCCGATCATCACCACTGTACCCGGCTGAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCCAAAATATGTTTACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAgg                                                                             <  2:366106/140‑1 (MQ=255)                                                       gCGTAAGCACCCATCATCCCCAGACTATCCGGTTTAACTTGCGGGTTAACAGAAACCACCGTGGCGATAATCTGCACTTTGTTGACGAAGCCGTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGt                                                                           <  2:571806/143‑1 (MQ=255)                                                             gCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGGATCTGGCGGCCAATCAGACGCGCGATCAGGGGTTGGcc                                                                     >  2:289972/1‑143 (MQ=255)                                                             gCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGATAACAGAAACCACGGTGACGATAACCTGAACGTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGGATCGGGCGGGCAATCAGACGCGCGATCAgggggtcgac                                                                     >  2:491317/1‑136 (MQ=255)                                                               aCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCtt                                                                   >  1:624311/1‑143 (MQ=255)                                                                  aaTCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGc                                                                <  1:852566/143‑1 (MQ=255)                                                                         gcgACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGggcg                                                         <  2:92650/143‑1 (MQ=255)                                                                         gcgACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGAGCGGGGTTTCGCCTTCGCTTGggcg                                                         >  1:8603/1‑143 (MQ=255)                                                                             cgaTATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCt                                                     >  1:145595/1‑143 (MQ=255)                                                                               atatCCGGGTTAACTGGCGGGTTATGGTAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTc                                                   <  2:434507/143‑1 (MQ=255)                                                                                  tCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCGTGGGCGGCCtccccg                                                >  2:539749/1‑138 (MQ=255)                                                                                           acgtGCTGGTTAACAGGAAGCACGGTGGCGATAACCTGAACTTTGCTGACGAAGCCTGCCGGGAACAGCGCGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGgaagaa                                       <  1:837205/140‑1 (MQ=255)                                                                                             ttGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAGTCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGGAGc                                     >  1:264346/1‑143 (MQ=255)                                                                                                  ggTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCAGGAAGGAAGAAGCTAccc                                >  2:417984/1‑143 (MQ=255)                                                                                                   gTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGAAGGAAGAAGCTACCCg                               >  1:270925/1‑143 (MQ=255)                                                                                                           aaaCCACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACGGCGGGCGAATCGGGCGGACAATCAGACGCGCGATCAGGGTTTCGCCGTCGCTTGGGCGGCCTTCACGCCGGAAGAAGCTCCCCGGGATACGa                       >  1:476498/1‑143 (MQ=255)                                                                                                            aaCCACGCGGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGAc                      <  2:510404/143‑1 (MQ=255)                                                                                                              ccACGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACca                    <  1:144360/143‑1 (MQ=255)                                                                                                                aCGGTGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACcagc                  >  1:438314/1‑143 (MQ=255)                                                                                                                    tGGCGATAACCTGAACTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGCTACGACCAGCAGCg              >  2:381624/1‑143 (MQ=255)                                                                                                                               tGAACTTTGTTGACGGAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACcag                   >  2:366697/1‑127 (MQ=255)                                                                                                                               tGAACTTTGTTGACGGAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACcag                   <  1:366697/127‑1 (MQ=255)                                                                                                                                ggaCTTTGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTcc  <  1:727705/141‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         ACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCC  >  NC_000913/3310708‑3310977

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