Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 961,268 A→T 9.1% K23M (AAG→ATG)  cmk → cytidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913961,2680AT 9.1% 122.2 / 4.4 44K23M (AAG→ATG) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base A (25/15);  new base T (2/2);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATC  >  NC_000913/961143‑961399
                                                                                                                             |                                                                                                                                   
tCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAg                                                                                                                       >  1:945690/1‑140 (MQ=255)
  ggCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCg                                                                                                                     <  1:535759/140‑1 (MQ=255)
       gTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCa                                                                                                                <  1:856381/140‑1 (MQ=255)
       gTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCa                                                                                                                <  2:808723/140‑1 (MQ=255)
               cgcgcgGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATc                                                                                                        >  2:348407/1‑140 (MQ=255)
                 cgcgGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATctg                                                                                                      <  2:65248/140‑1 (MQ=255)
                        aTGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACt                                                                                               <  2:895146/140‑1 (MQ=255)
                              aCGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTTCCATTGATGGCCCAATCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTg                                                                                         >  2:322501/1‑140 (MQ=255)
                                                    aTAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCCGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCGGGACTCGGGTGCAATTTATCGCGTACTGGCAtt                                                                   >  2:1197883/1‑140 (MQ=255)
                                                    aTAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAATGGCACGTTGTGTAAGGCTATGGCGGTAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCAtt                                                                   >  1:251302/1‑140 (MQ=255)
                                                    aTAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCAtt                                                                   >  2:1031243/1‑140 (MQ=255)
                                                    aTAAAGATGACGGCAATGGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCAtt                                                                   <  2:1192275/140‑1 (MQ=255)
                                                            gACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCAt                                                           <  1:176031/140‑1 (MQ=255)
                                                            gACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCAt                                                           <  1:175862/140‑1 (MQ=255)
                                                            gACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCAt                                                           >  1:830522/1‑140 (MQ=255)
                                                            gACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCAt                                                           >  1:831676/1‑140 (MQ=255)
                                                              cGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTa                                                         >  1:669397/1‑140 (MQ=255)
                                                              cGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTa                                                         >  1:143606/1‑140 (MQ=255)
                                                                    ttGGCCCGGCTATTACCATTGATGGCGCCAGCGGGGCAGGGAAAGGCACCTTGTGTTAGGCTTTGGCGGAAGCGTTGCAAGGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCAGTGGCGGCATTACATCAc                                                   <  2:854166/140‑1 (MQ=255)
                                                                    ttGCACCGGTTATGACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCAc                                                   <  2:859609/140‑1 (MQ=255)
                                                                         ccGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgt                                              <  1:1173073/140‑1 (MQ=255)
                                                                                      ttGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCg                                 >  1:589715/1‑140 (MQ=255)
                                                                                      ttGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCg                                 >  1:1017037/1‑140 (MQ=255)
                                                                                          tGGCCCAAGCGGTGCAGGGAAAGGCACCGTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAg                             >  2:918788/1‑140 (MQ=255)
                                                                                          tGGCCAGAGCGGTGCAGGGATAGGCTCCTTTGGTTAGGCTATGGCGGAAGCGTTGCACTGGCATCTTCTGGGCTCGGGTGCAATTTATCGCGTTCTGTCATTGGCGGCTTTCCAGCCCCATCTGGAGGTTGCATCGTAAg                             >  2:1225998/1‑140 (MQ=255)
                                                                                            gCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGAt                           <  2:856034/140‑1 (MQ=255)
                                                                                            gCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGAt                           <  2:826838/140‑1 (MQ=255)
                                                                                               cAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATgcg                        <  2:812076/140‑1 (MQ=255)
                                                                                                aaGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATgcgc                       >  1:805656/1‑140 (MQ=255)
                                                                                                       gCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGAGGTTGCGTCGGAAGATGCGCTGGTAcc                >  1:289520/1‑140 (MQ=255)
                                                                                                        cAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCg               >  2:304374/1‑140 (MQ=255)
                                                                                                        cAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCg               >  1:519414/1‑140 (MQ=255)
                                                                                                          gggAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCt             <  1:386784/140‑1 (MQ=255)
                                                                                                                ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc       >  2:583688/1‑140 (MQ=255)
                                                                                                                ggCACCTTGGGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc       >  2:1107746/1‑140 (MQ=255)
                                                                                                                  cACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc     >  2:385585/1‑140 (MQ=255)
                                                                                                                  cACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc     >  1:1167020/1‑140 (MQ=255)
                                                                                                                  cACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc     >  1:1166811/1‑140 (MQ=255)
                                                                                                                  cACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc     >  1:107479/1‑140 (MQ=255)
                                                                                                                     cTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATc  >  2:1154711/1‑140 (MQ=255)
                                                                                                                       tgtgTATGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCAc                                                   >  1:145610/1‑89 (MQ=255)
                                                                                                                       tgtgTATGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCAc                                                   <  2:145610/89‑1 (MQ=255)
                                                                                                                       tgtgTATGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCAc                                                   <  2:145508/89‑1 (MQ=255)
                                                                                                                       tgtgTATGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCAc                                                   >  1:145508/1‑89 (MQ=255)
                                                                                                                             |                                                                                                                                   
TCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATC  >  NC_000913/961143‑961399

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