Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,994 C→G 7.5% M235I (ATG→ATC)  PP_3718 ← aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,240,9940CG 7.5% 163.3 / 4.2 53M235I (ATG→ATC) PP_3718aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base G (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

TTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACGTGGT  >  NC_002947/4240853‑4241139
                                                                                                                                             |                                                                                                                                                  
ttCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCGA‑TCAtgg                                                                                                                                            >  7:262992/1‑147 (MQ=255)
    tCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGAtcggctc                                                                                                                                        >  1:280176/1‑149 (MQ=255)
    tCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGAtcggctc                                                                                                                                        >  4:195104/1‑149 (MQ=255)
       cTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCGA‑TCAtgggctcggg                                                                                                                                     <  8:262992/149‑2 (MQ=255)
             aagaagGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaagaa                                                                                                                               >  8:162208/1‑149 (MQ=255)
             aagaagGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaagaa                                                                                                                               >  7:20827/1‑149 (MQ=255)
                  ggTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATgccg                                                                                                                          >  3:78262/1‑149 (MQ=255)
                      cccAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCAc                                                                                                                      >  8:250547/1‑149 (MQ=255)
                          ggCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTg                                                                                                                  >  7:315498/1‑149 (MQ=255)
                          ggCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTg                                                                                                                  <  7:162208/149‑1 (MQ=255)
                             gcCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCg                                                                                                               >  2:83220/1‑149 (MQ=255)
                             gcCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCg                                                                                                               >  7:171075/1‑149 (MQ=255)
                                cGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGc                                                                                                             >  2:311318/1‑148 (MQ=255)
                                       ggTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCg                                                                                                     <  7:250547/149‑1 (MQ=255)
                                           accacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCACTGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                                  <  1:83220/149‑1 (MQ=255)
                                             cacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                                  <  3:256829/146‑1 (MQ=255)
                                             cacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                                  >  4:256829/1‑146 (MQ=255)
                                                ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaa                                                                                                                                  >  5:29523/1‑111 (MQ=255)
                                                ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaa                                                                                                                                  <  6:29523/111‑1 (MQ=255)
                                                ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATgccgcc                                                                                                                        >  5:310347/1‑121 (MQ=255)
                                                ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATgccgcc                                                                                                                        <  6:310347/121‑1 (MQ=255)
                                                  cgtcgGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCtt                                                                                           <  2:280176/148‑1 (MQ=255)
                                                    tcgGCGACGAACAGGATGTCTTAGGTCTGGCACAGCTGCCAGATGCGCTGGAAGTAGCCCTCGGGCGTAATGATCACCCCGCCCGAGCCCA‑TTATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTc                                                                                        <  2:306616/149‑1 (MQ=255)
                                                      ggCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTcc                                                                                       <  8:171075/148‑1 (MQ=255)
                                                       gCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTc                                                                                     <  3:136433/149‑1 (MQ=255)
                                                         gacgaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCg                                                                                    >  3:282280/1‑148 (MQ=255)
                                                              aCAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCg                                                                              >  5:138675/1‑149 (MQ=255)
                                                                aGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGc                                                                            <  1:311318/149‑1 (MQ=255)
                                                                        gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTc                                                                                     >  4:129370/1‑132 (MQ=255)
                                                                            gTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATc                                                                                             >  8:156501/1‑120 (MQ=255)
                                                                            gTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATc                                                                                             <  7:156501/120‑1 (MQ=255)
                                                                             tCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAATAATGCCGCCACGTTGTCGGCACCCAGCgagatg                                                                                                <  7:197567/116‑3 (MQ=255)
                                                                             tCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAgg                                                                                                >  8:197567/1‑116 (MQ=255)
                                                                                 gCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGa                                                              <  6:138675/146‑1 (MQ=255)
                                                                                 gCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAgg                                                            <  8:315498/148‑1 (MQ=255)
                                                                                 gCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAg                                                             <  4:282280/147‑1 (MQ=255)
                                                                                                   cTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATgt                                         <  6:74485/149‑1 (MQ=255)
                                                                                                         gTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTg                                   <  5:132598/149‑1 (MQ=255)
                                                                                                                  gggCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCgg                           >  4:251162/1‑148 (MQ=255)
                                                                                                                      ggAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAAc                                                                                 >  6:13431/1‑90 (MQ=255)
                                                                                                                      ggAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAAc                                                                                 <  5:13431/90‑1 (MQ=255)
                                                                                                                       gAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGCCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGcc                                                   >  6:153858/1‑119 (MQ=255)
                                                                                                                       gAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGCCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGcc                                                   <  5:153858/119‑1 (MQ=255)
                                                                                                                        aaTGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggt                    >  6:221605/1‑149 (MQ=255)
                                                                                                                          tGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttg                  >  6:270314/1‑149 (MQ=255)
                                                                                                                            aTCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttggg                >  4:130834/1‑149 (MQ=255)
                                                                                                                               aCCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTAc                        <  3:158588/138‑1 (MQ=255)
                                                                                                                                ccccGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttggggttg            >  5:38734/1‑149 (MQ=255)
                                                                                                                                  ttcgccgAGCCGA‑TCATGGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTcc                                                      <  1:173829/101‑1 (MQ=255)
                                                                                                                                  ttcgccgAGCCGA‑TCATGGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTcc                                                      >  2:173829/5‑105 (MQ=255)
                                                                                                                                         aGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACgtgg   <  5:270314/149‑1 (MQ=255)
                                                                                                                                          gCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACgtggt  <  6:38734/149‑1 (MQ=255)
                                                                                                                                          gCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACgtggt  <  5:221605/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACGTGGT  >  NC_002947/4240853‑4241139

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