Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,484,163 C→A 8.7% T165N (ACC→AAC)  PP_3976 → isochorismatase family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,484,1630CA 8.7% 141.2 / 4.0 46T165N (ACC→AAC) PP_3976isochorismatase family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (17/25);  new base A (2/2);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

GGTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCA  >  NC_002947/4484017‑4484295
                                                                                                                                                   |                                                                                                                                    
ggTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCt                                                                                                                                    <  6:247885/149‑1 (MQ=255)
 gTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTg                                                                                                                                   <  5:293615/149‑1 (MQ=255)
      cAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGgcagc                                                                                                                              <  4:246245/149‑1 (MQ=255)
            caccTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAggag                                                                                                                         >  1:243606/1‑146 (MQ=255)
               cTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTAc                                                                                                                     >  5:58600/1‑149 (MQ=255)
                      cTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTg                                                                                                              >  5:163356/1‑149 (MQ=255)
                         tCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTgg                                                                                                             <  2:221547/147‑1 (MQ=255)
                           cGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAgc                                                                                                         >  3:179069/1‑149 (MQ=255)
                               aTGCACGATGCAGGGTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCTAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGa                                                                                                     <  2:159741/149‑1 (MQ=255)
                               aTGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGa                                                                                                     <  2:69418/149‑1 (MQ=255)
                                  cACGTTGCAGGATATGAGATCTATATGGTTGCCGACTCTTACTGGCGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACt                                                                                                   <  4:171383/148‑1 (MQ=255)
                                  cACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACt                                                                                                   <  6:131402/148‑1 (MQ=255)
                                   aCGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGAACAATGAGGCAA‑ACCACTACCCCATTCAGCGCAAGATCCACGCGGGCGGGGTTCCAATTTCCATGCAGCAGGTGCAACTTGACTTGCAGCCCGACTg                                                                                                  >  1:80092/1‑148 (MQ=255)
                                         cAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGaaaaa                                                                                            <  5:96168/148‑1 (MQ=255)
                                            gCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTAcc                                                                                                                                     >  8:176107/1‑104 (MQ=255)
                                            gCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACcgc                                                                                        >  7:235442/1‑149 (MQ=255)
                                            gCTATGAGATCTATATGGTCGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACcgc                                                                                        <  7:252380/149‑1 (MQ=255)
                                              tATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGt                                                                                                               <  3:261595/124‑1 (MQ=255)
                                              tATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGt                                                                                                               >  4:261595/1‑124 (MQ=255)
                                                tGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGcagcag                                                                                                                            >  6:291261/1‑109 (MQ=255)
                                                tGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGcagcag                                                                                                                            <  5:291261/109‑1 (MQ=255)
                                                        atatGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAAGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGa                                                                             <  6:283164/148‑1 (MQ=255)
                                                             gTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGt                                                                       >  2:92784/1‑149 (MQ=255)
                                                             gTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGt                                                                       >  2:125952/1‑149 (MQ=255)
                                                                 ccGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATg                                                                   <  6:163356/149‑1 (MQ=255)
                                                                      gCTTCCGGGGGCACGACCAAGGAGGCTCTACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGt                                                                       >  1:19188/1‑141 (MQ=255)
                                                                      gCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGt                                                                       <  2:19188/140‑1 (MQ=255)
                                                                      gCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCt                                                              >  8:139325/1‑149 (MQ=255)
                                                                             gtggCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCaaa                                                       <  6:58600/147‑1 (MQ=255)
                                                                               ggCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCTCATGATCCAGGCGTGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAaga                                                     <  6:258511/149‑1 (MQ=255)
                                                                                    gACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCAtt                                                <  7:139325/149‑1 (MQ=255)
                                                                                         aggaggCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGgc                                           <  4:4978/149‑1 (MQ=255)
                                                                                                tC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACg                                    >  3:131262/1‑149 (MQ=255)
                                                                                                   aCGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGt                                                                       <  6:129198/112‑1 (MQ=255)
                                                                                                   aCGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGt                                                                       >  5:129198/1‑112 (MQ=255)
                                                                                                   aCGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGc                                  <  2:80092/149‑1 (MQ=255)
                                                                                                     gACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTg                                                                                                                   <  3:21936/66‑1 (MQ=255)
                                                                                                     gACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTg                                                                                                                   >  4:21936/1‑66 (MQ=255)
                                                                                                       cTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATgg                              >  1:130672/1‑149 (MQ=255)
                                                                                                             cATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTgg                                                            >  6:229749/1‑113 (MQ=255)
                                                                                                             cATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTgg                                                            <  5:229749/113‑1 (MQ=255)
                                                                                                                      cATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGc                 <  1:202820/147‑1 (MQ=255)
                                                                                                                              aGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCTTGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATg       <  1:92784/149‑1 (MQ=255)
                                                                                                                                     cctgctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCAtt                                                <  6:4185/96‑1 (MQ=255)
                                                                                                                                     cctgctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCAtt                                                >  5:4185/6‑101 (MQ=255)
                                                                                                                                      ctgctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGca  >  7:16989/5‑146 (MQ=255)
                                                                                                                                      ctgctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGca  <  8:16989/142‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                    
GGTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTC‑ACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCA  >  NC_002947/4484017‑4484295

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