Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	961,476	T→C	R92R (CGT→CGC)	cmk →	cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	961,476	0	T	C	97.9%	161.1 / NA	48	R92R (CGT→CGC)	cmk	cytidylate kinase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (22/25);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTA  >  NC_000913/961344‑961611                                                                                                                                     |                                                                                                                                        aCATCACCATGTTGAAGTTGCGTCGTAAGATGCGCTGTTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAgg                                                                                                                                  <  1:215668/140‑1 (MQ=255)   aTCACCATGTTGATGTTACGTCCGAAGATGCGCTGGTACTGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGaa                                                                                                                                <  2:286348/140‑1 (MQ=255)         atgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCg                                                                                                                          >  2:17085/1‑140 (MQ=255)          tgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAGGATGTCACCGGCGAGATTCGCGCTCAGGAAGTGGCGa                                                                                                                         >  1:300232/1‑140 (MQ=255)          tgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGTGAAATTCGCACTCAGGAAGTGGCGa                                                                                                                         >  1:295799/1‑140 (MQ=255)          tgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCAc                                                                                                                                       >  1:616002/1‑126 (MQ=255)          tgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCAc                                                                                                                                       <  2:616002/126‑1 (MQ=255)          tgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGa                                                                                                                         >  1:439978/1‑140 (MQ=255)          tgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTAAGCGGCGAAATTCGCACTCAGGAAGTGGCGa                                                                                                                         >  1:583955/1‑140 (MQ=255)            ttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAAt                                                                                                                       >  1:37801/1‑140 (MQ=255)             tgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATg                                                                                                                      <  2:37801/140‑1 (MQ=255)             tgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATg                                                                                                                      <  2:556707/140‑1 (MQ=255)                tgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGGGAAATTCGTACTCAGGAAGTGGCGAATGCAg                                                                                                                   >  2:261642/1‑140 (MQ=255)                    gCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTc                                                                                                               >  2:355661/1‑140 (MQ=255)                              aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAGGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTcgcg                                                                                                     <  2:78158/140‑1 (MQ=255)                              aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTcgcg                                                                                                     <  1:390514/140‑1 (MQ=255)                              aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTcgcg                                                                                                     <  1:51221/140‑1 (MQ=255)                                gcgcTGGTACCGCTGGCATCCCATCTGGAGGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGc                                                                                                   <  1:218000/140‑1 (MQ=255)                                         ccGCTGGCATCCCATCTGGATGTACTTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCAcg                                                                                          <  1:232926/140‑1 (MQ=255)                                          cGCTGGCATCCCATCTGTATGTTCGTTTTGTATCGACCACTGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCAcgc                                                                                         <  2:300232/140‑1 (MQ=255)                                            cTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGt                                                                                       >  2:87985/1‑140 (MQ=255)                                               gCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCg                                                                                    <  1:87985/140‑1 (MQ=255)                                                cATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGt                                                                                   >  2:547819/1‑140 (MQ=255)                                                 aTCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTg                                                                                  >  2:6626/1‑140 (MQ=255)                                                 aTCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTg                                                                                  <  1:368200/140‑1 (MQ=255)                                                    ccATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACCCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAg                                                                               <  2:534478/140‑1 (MQ=255)                                                             tGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATCCAGCTTCACAAGTCGCCGCATTCCCCCGCGTTCGTGAAGCATTATTGc                                                                      >  2:260999/1‑140 (MQ=255)                                                                aCGTTTTGTGTCGACTAATGGCCATTTGGTAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTc                                                                   <  1:180509/140‑1 (MQ=255)                                                                     ttGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCaa                                                              >  2:621762/1‑140 (MQ=255)                                                                      tgtgTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAAc                                                             >  2:378701/1‑140 (MQ=255)                                                                              ccAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGttt                                                     <  1:607887/140‑1 (MQ=255)                                                                                  tGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTGCGTGAAGGATTTTTTCGTCGCAAACGCGCGTTTcgcg                                                 >  1:201371/1‑140 (MQ=255)                                                                                   ggCAATCTGGAAGTGATCCTCGCAGGGGAAGATGTTAGCCGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGa                                                <  2:201371/140‑1 (MQ=255)                                                                                                  aTCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGGAGCATTATGGCGTCGCGAACGCGCGTTTCGTGAATTACCAGGTCTGAt                                 >  2:190048/1‑140 (MQ=255)                                                                                                  aTCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGAt                                 >  1:149151/1‑140 (MQ=255)                                                                                                            gggAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGc                       >  1:151098/1‑140 (MQ=255)                                                                                                              gAAGATGTCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCcg                     <  1:643450/140‑1 (MQ=255)                                                                                                                     tCAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATg              <  1:1571/140‑1 (MQ=255)                                                                                                                      cAGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGAGAATTACGAGGGCTGATTGCCGATGGCCGCGACATgg             >  1:211417/1‑140 (MQ=255)                                                                                                                       aGCGGCGAAATTCGCACTCAGGCAGCGGCCAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATggg            <  1:179871/140‑1 (MQ=255)                                                                                                                       aGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATggg            <  2:540687/140‑1 (MQ=255)                                                                                                                       aGCGGCGAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCCCGAATTACCAGGTCTGCTTGCTGATGGCCGCGACATggg            >  2:217872/1‑140 (MQ=255)                                                                                                                           gcgAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACg        <  1:381661/140‑1 (MQ=255)                                                                                                                           gcgAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACg        <  2:18291/140‑1 (MQ=255)                                                                                                                            cgAAATTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGGTTCGCGAAATACCAGGCCTGATTGCCGATGGCCGCGACATGGGAACgg       >  2:285548/1‑140 (MQ=255)                                                                                                                              aaaTTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTCCGCGAATTGCCAGGTCCGATTGCCGGTGGCTGCGACATGGGAACggtg     >  2:259980/1‑140 (MQ=255)                                                                                                                                aTTCGCACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACggtggt   <  1:621762/140‑1 (MQ=255)                                                                                                                                 ttCGCACTCAGGAAGTGCCGAATGCAGCTTCTCAATTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTa  <  2:211417/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        ACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTA  >  NC_000913/961344‑961611

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