Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,235,173 C→G 8.3% A185A (GCG→GCC)  tdk ← thymidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,235,1730CG 8.3% 130.1 / 2.4 47A185A (GCG→GCC) tdkthymidine kinase
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base G (2/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTTATCTGAATCATTTACAACATCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCG  >  CP000730/2235038‑2235312
                                                                                                                                       |                                                                                                                                           
tttATCTGAATCATTTACAACATCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGcgcgcg                                                                                                                                        >  2:749602/1‑141 (MQ=255)
 ttATCTGAATCATTTACAACATCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGa                                                                                                                                       <  1:391273/141‑1 (MQ=255)
 ttATCTGAATCATTTACAACATCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAACTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATTTTGcgcgcgg                                                                                                                                       >  1:511919/1‑140 (MQ=255)
       gAATCATTTACAACATCTGGGTCACTTAACAGTTCCTTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGt                                                                                                                                 <  1:327643/141‑1 (MQ=255)
               tACAACATCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATACCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGACGCGGAGCGCTCGCAAACTTGGTTCATa                                                                                                                         >  1:465470/1‑141 (MQ=255)
                 caacaTCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGCTATGGCGCGCGCGACAACGCGGCTCacgac                                                                                                                       >  2:480649/1‑137 (MQ=255)
                  aacaTCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACt                                                                                                                      <  1:249294/141‑1 (MQ=255)
                   acaTCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTTCTACGATGTGGTGGCCGCTGCAACGTGGCTCCTAACtt                                                                                                                     >  2:390917/1‑141 (MQ=255)
                   acaTCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGGGCGCGACAACGTGGCTCATAACtt                                                                                                                     >  2:95240/1‑141 (MQ=255)
                   acaACTGGATCACTTAACAGTTCCTTTTACTGGTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACtt                                                                                                                     <  1:356221/141‑1 (MQ=255)
                        tGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCAtt                                                                                                                <  1:675096/141‑1 (MQ=255)
                         gggTCACTTAATAGTGCATTGAACTGTTCGTATCTGTCTTCAACAATATCTAAGTGATCAAACACTTATAATTCCGCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTg                                                                                                               <  1:685315/141‑1 (MQ=255)
                           ggcACTTAACAGTTCATTTAACTGTTCGTATCTTTCCTCTACAATATCTAATTGATCACACAATCATAATTCTTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGc                                                                                                             <  1:361294/139‑1 (MQ=255)
                           gTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGc                                                                                                             <  1:186203/141‑1 (MQ=255)
                           gTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGc                                                                                                             <  1:81507/141‑1 (MQ=255)
                             cACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCAc                                                                                                           <  1:5787/141‑1 (MQ=255)
                              aCTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCAcc                                                                                                          >  2:654060/1‑141 (MQ=255)
                                   acagttaaatgaaCTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATAGTGCGCGCGACAACGTGGCTCATAACTTTCATTCGCACatacta                                                                                                     >  2:387631/12‑141 (MQ=255)
                                         cATTTAACTGTTCGTATCTTTCTTCTACAATATCTAAGTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATaa                                                                                               <  2:82436/141‑1 (MQ=255)
                                           tttAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACATATAATTCCTCCTTATTATTATCACTAGGTGCTACGCTATGGTGAGCGCGACAACATGGCTCATAACTTTCTTTGGCACCTACTAAGATACTc                                                                                             >  2:714146/1‑141 (MQ=255)
                                                cTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCCCCTACTCAGATAATCGGatc                                                                                        >  2:295377/1‑141 (MQ=255)
                                                     cGTATCTTTCTTCTACAATATCTAATTGGTCCAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCTACAACATGGCTACTAACTTTCATAGGCAACAACTAAGATAATCGGatcatcta                                                                                   >  1:754601/1‑139 (MQ=255)
                                                               ttctATAATATCTAATTGATCAGACACTTATAATTTCTAATTCTTCTTCCCACTTGGTTCTTCGATATGGTGCGCGCGACAACGTCGCTCATATCTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTgg                                                                         <  1:432391/141‑1 (MQ=255)
                                                                       tatCTAATTCATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCAt                                                                 <  2:449596/141‑1 (MQ=255)
                                                                            ttttGATCAAACTCTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCAGGTTTACCATttatt                                                            <  1:380493/139‑1 (MQ=255)
                                                                             aTTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATttataa                                                           >  2:645146/1‑139 (MQ=255)
                                                                                            acaaTCCCTCCTTATTATTATCACTAGGGGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTctac                                            <  1:403213/139‑1 (MQ=255)
                                                                                              aaTTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTctacta                                          >  2:658182/1‑141 (MQ=255)
                                                                                              aaTTCCTCCTTATTATTATCAATAGGTGCTACGATATGGTGCGCGCGACAACGTGTCTCATAAATTTCTTTGGCACCTACTAAGATAATCTGATCAACGTTTTTAGCTGGTTTCACCTTTATTAACCGTTGCGTTtgacta                                          >  2:130487/1‑141 (MQ=255)
                                                                                                  cctcctTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTGATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATg                                      <  1:51332/141‑1 (MQ=255)
                                                                                                     cctTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATgaag                                   <  1:450319/141‑1 (MQ=255)
                                                                                                           tattatCACTAGGTGCTACGATATGGTGCGCGCGACAACGTAGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCAATTTTAGCTGGTTTACCATTTATTAATCGTTGCGGTCTACTTGATGAAGTACcac                             >  2:298239/1‑141 (MQ=255)
                                                                                                             ttatCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGGTCTACTAGAGGAAGAACCACaa                           >  2:745714/1‑141 (MQ=255)
                                                                                                                 cACTAGGTGCTACGATAGGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAg                       <  2:324302/141‑1 (MQ=255)
                                                                                                                        tGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCCTTAATTAATCGTTGCGTTCTACTAGATGAAAAACCACAAAAGGCCCAAAc                >  2:95541/1‑141 (MQ=255)
                                                                                                                            aCGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCt            <  1:654060/141‑1 (MQ=255)
                                                                                                                            aCGATATGGGGCGCGCGACAACGTGGCTCATAACTGTCATTAGCACCAACGAATGTAAGCGGGTCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCt            >  2:732728/1‑141 (MQ=255)
                                                                                                                              gATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTg          <  1:263534/141‑1 (MQ=255)
                                                                                                                                tatGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGaa        <  1:311245/141‑1 (MQ=255)
                                                                                                                                 atGGTGGGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAg       <  1:95240/141‑1 (MQ=255)
                                                                                                                                 atGGTGGGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAg       <  1:390917/141‑1 (MQ=255)
                                                                                                                                   ggTGCGGGCGAGAACGTGGCTCATAACTTTCATTGGCACCTACTAAGGTAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGtt     <  2:106816/141‑1 (MQ=255)
                                                                                                                                    gTGCGCGCGAGAACGTGGCTCATAACTATCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTCTTTTTTGTTTCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGttt    >  1:551456/1‑141 (MQ=255)
                                                                                                                                    gTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACTTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAGGttt    >  2:149345/1‑141 (MQ=255)
                                                                                                                                    gTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGttt    >  1:499704/1‑141 (MQ=255)
                                                                                                                                    gTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGttt    >  1:50399/1‑141 (MQ=255)
                                                                                                                                    gTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGttt    >  2:577921/1‑141 (MQ=255)
                                                                                                                                    gTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGttt    >  1:499430/1‑141 (MQ=255)
                                                                                                                                     tGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATAGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTc   <  2:500678/141‑1 (MQ=255)
                                                                                                                                      gcgcgcGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCg  >  1:222936/1‑141 (MQ=255)
                                                                                                                                       cgcgcgACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACaa                  >  2:680146/1‑124 (MQ=255)
                                                                                                                                       cgcgcgACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACaa                  <  1:680146/124‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                           
TTTATCTGAATCATTTACAACATCTGGGTCACTTAACAGTTCATTTAACTGTTCGTATCTTTCTTCTACAATATCTAATTGATCAAACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCG  >  CP000730/2235038‑2235312

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