Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,887,744 C→T 7.2% A172A (GCG→GCA)  maa ← acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,887,7440CT 7.2% 163.8 / 2.1 55A172A (GCG→GCA) maaacetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base T (2/2);  total (30/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 = 8.70e-01

GTGTGCAACGCCAGCAACAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCTCTCCAGCC  >  NC_002947/5887598‑5887886
                                                                                                                                                  |                                                                                                                                              
gtgtGCAAAGCCAGCAACAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTcgc                                                                                                                                              >  3:94648/1‑149 (MQ=255)
   tgCAACGCCAGCAACAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGcc                                                                                                                                           >  7:217072/1‑149 (MQ=255)
        cGCCAGCAACAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCgg                                                                                                                                      <  4:94648/149‑1 (MQ=255)
        cGCCAGCAACAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCgg                                                                                                                                      <  3:263242/149‑1 (MQ=255)
               aaCAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATc                                                                                                                               <  4:157558/149‑1 (MQ=255)
               aaCAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCCCGCCCGCGGGCACCTc                                                                                                                               >  7:16505/1‑149 (MQ=255)
                    cTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCggg                                                                                                                          >  2:196031/1‑149 (MQ=255)
                        aTCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTcac                                                                                                                      <  8:69595/149‑1 (MQ=255)
                            aGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACcaca                                                                                                                  >  7:333780/1‑149 (MQ=255)
                              acAAACCTCAAGGCTTTGCTCGACGGCCGCCAGTTCAACGCCTCCAATACCCTGGCGTTAGGCTGCTACTGCCCTTCGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGTGCACATCCCGGGTCACCACACt                                                                                                                <  8:276630/148‑1 (MQ=255)
                                     tCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATccc                                                                                                                             <  4:178632/129‑1 (MQ=255)
                                     tCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATccc                                                                                                                             >  3:178632/1‑129 (MQ=255)
                                     tCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGgcg                                                                                                         >  6:182328/1‑149 (MQ=255)
                                     tCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGgcg                                                                                                         >  3:249866/1‑149 (MQ=255)
                                      cAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCggg                                                                                                                          <  2:213007/131‑1 (MQ=255)
                                      cAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCggg                                                                                                                          >  1:213007/1‑131 (MQ=255)
                                       aaGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGgcgc                                                                                                        <  4:114506/148‑1 (MQ=255)
                                       aaGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGcc                                                                                                       <  8:333780/149‑1 (MQ=255)
                                            tttGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCcgacga                                                                                                  >  2:39813/1‑149 (MQ=255)
                                              tGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGAt                                                                                                 <  7:212238/148‑1 (MQ=255)
                                              tGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGAt                                                                                                 <  7:348934/148‑1 (MQ=255)
                                               gCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATg                                                                                                <  8:217072/148‑1 (MQ=255)
                                                cTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTTGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATgg                                                                                               >  7:12636/1‑148 (MQ=255)
                                                cTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTTGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATgg                                                                                               >  7:7188/1‑148 (MQ=255)
                                                cTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGc                                                                                              <  1:39813/149‑1 (MQ=255)
                                                cTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCCGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATgg                                                                                               >  6:318933/1‑148 (MQ=255)
                                                 tCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCa                                                                                             <  4:249866/149‑1 (MQ=255)
                                                  cGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCAt                                                                                            >  4:293943/1‑149 (MQ=255)
                                                   gACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCAt                                                                                            >  2:298851/1‑148 (MQ=255)
                                                         cgccAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGATGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCAcc                                                                                     >  1:286770/1‑149 (MQ=255)
                                                          gccAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCAcc                                                                                     >  4:151871/1‑148 (MQ=255)
                                                             aGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTTGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATg                                                                                 <  8:7188/149‑1 (MQ=255)
                                                             aGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATg                                                                                 <  7:95474/149‑1 (MQ=255)
                                                             aGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATg                                                                                 <  3:151871/149‑1 (MQ=255)
                                                             aGGGCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTTGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATg                                                                                 <  8:12636/149‑1 (MQ=255)
                                                                        tcaaaTACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCACGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCTATGGTCACCCCCgg                                                                      >  6:356140/4‑149 (MQ=255)
                                                                        tCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCgg                                                                      >  1:132976/1‑149 (MQ=255)
                                                                           aaTACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAg                                                                   <  1:196031/149‑1 (MQ=255)
                                                                                               ctACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCaa                                               >  7:347744/1‑149 (MQ=255)
                                                                                                  cTGCCCTTCGTCGGGCTGACGGACCCTCGCCGGCTTACCCACCACCGTCGCGCCCGCGGGCCCAACCCGGGTCACCACACTGCGGGCGCCGACTATGGCATAATCACCGCTGGTCCCCCCCGGCAGGATGATCGCCGCACCGCCCATcc                                            >  5:267494/1‑149 (MQ=255)
                                                                                                      ccTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCacac                                        <  2:132976/149‑1 (MQ=255)
                                                                                                        ttGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGt                                      >  7:140665/1‑149 (MQ=255)
                                                                                                              gggCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGc                                >  7:76068/1‑149 (MQ=255)
                                                                                                                    aCGGACCCTCGCGGGGTTAACCACCTCGGTCGCGCCGGGGGGGAGATCCCGGGTCACCACACTGCTGGCGCCGACGATGGCATTATCACCGATGGTCACCCGCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGAt                            <  6:267494/147‑1 (MQ=255)
                                                                                                                    aCGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATg                           <  8:347744/148‑1 (MQ=255)
                                                                                                                    aCGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCCGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGAt                            <  5:318933/147‑1 (MQ=255)
                                                                                                                            tCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATcgccg                                                        <  8:81550/111‑1 (MQ=255)
                                                                                                                            tCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATcgccg                                                        >  7:81550/1‑111 (MQ=255)
                                                                                                                            tCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGgc                  >  5:101924/1‑149 (MQ=255)
                                                                                                                                      aCCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGctct        >  4:129200/1‑149 (MQ=255)
                                                                                                                                        ccaccaCCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCtctc       >  5:128093/1‑148 (MQ=255)
                                                                                                                                         caccacCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCtctc       >  5:201882/1‑147 (MQ=255)
                                                                                                                                          accaccGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCTCTCCAg    >  5:118140/1‑149 (MQ=255)
                                                                                                                                          accaccGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCTCTCCAg    >  1:68990/1‑149 (MQ=255)
                                                                                                                                           ccaccGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCTCTCCAGc   <  1:298851/149‑1 (MQ=255)
                                                                                                                                            caccGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCTCTCCAGcc  >  4:126899/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GTGTGCAACGCCAGCAACAGCTTCATCGAGCCAATCCTCAAGGCTTTGCTCGACGGCCGCCAGGTCAACGCCTCCAATACCCTGGCGGTAGGCTGCTACTGCCCTTGGTCGGGCTGACGGACCCTCGCCGGGTTACCCACCACCGTCGCGCCCGCGGGCACATCCCGGGTCACCACACTGCCGGCGCCGACGATGGCATTATCACCGATGGTCACCCCCGGCAGGATGATCGCCGCACCGCCAATCCACACGTTGTCGCCGATGGTTACCGGCCGCCCGCTCTCCAGCC  >  NC_002947/5887598‑5887886

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: ‑