Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,150,529 C→G 7.8% P44P (CCC→CCG)  PP_3655 → cytosine/purine/uracil/thiamine/allantoin permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,150,5290CG 7.8% 136.6 / 2.3 48P44P (CCC→CCG) PP_3655cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base G (2/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTGGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATC  >  NC_002947/4150383‑4150676
                                                                                                                                                   |                                                                                                                                                   
gccttgGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCa                                                                                                                                                    <  2:19572/148‑1 (MQ=255)
        gtcagtcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCg                                                                                                                                           >  3:212550/1‑149 (MQ=255)
           agtcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgcg                                                                                                                                        >  2:105456/1‑149 (MQ=255)
             tcaTGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGgcgcg                                                                                                                                      >  7:238098/1‑149 (MQ=255)
                tGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTg                                                                                                                                   <  1:29825/149‑1 (MQ=255)
                     aGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCGATCGTCCCCGGCGCGCTGGgtgt                                                                                                                              >  3:195644/1‑149 (MQ=255)
                          aGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGc                                                                                                                         >  6:286822/1‑149 (MQ=255)
                                        cccTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCGATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAAc                                                                                                           <  4:195644/149‑1 (MQ=255)
                                        cccTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGgt                                                                                                                                >  6:37863/1‑128 (MQ=255)
                                         ccTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGt                                                                                                                      <  8:103457/137‑1 (MQ=255)
                                         ccTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGt                                                                                                                      >  7:103457/1‑137 (MQ=255)
                                          cTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCa                                                                                                                  >  2:154157/1‑140 (MQ=255)
                                          cTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCa                                                                                                                  <  1:154157/140‑1 (MQ=255)
                                             aTAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTg                                                                                                        <  4:212550/147‑1 (MQ=255)
                                             aTAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTg                                                                                                        <  6:190009/147‑1 (MQ=255)
                                             aTAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTg                                                                                                        <  5:286822/147‑1 (MQ=255)
                                               aGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCACCATGGAAAGGTAAGGGAGTCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTa                                                                                                     >  1:252462/1‑149 (MQ=255)
                                                gAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGCTAAGGGA‑TCTGTTCACACTCTGGTACGCCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGAGCTGGGTGTGCAGCTGTTCCACCTCAACCTCGTAtc                                                                                                   >  3:233424/1‑148 (MQ=255)
                                                  aaaaCACACGATCGGCTACGTTCACCCGCAAGACCGGCATGTTAAGGTATGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAAGCTGGTATggg                                                                                                 <  4:233424/149‑1 (MQ=255)
                                                   aaaCACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATgggg                                                                                                >  3:122051/1‑149 (MQ=255)
                                                          cGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCAtcgtc                                                                                         <  4:122051/149‑1 (MQ=255)
                                                           gATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTgg                                                                                                       <  6:240883/134‑1 (MQ=255)
                                                           gATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTgg                                                                                                       >  5:240883/1‑134 (MQ=255)
                                                           gATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGCGCAtcgtcg                                                                                        >  3:284331/1‑149 (MQ=255)
                                                                 cTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATcc                                                                                  <  2:115300/149‑1 (MQ=255)
                                                                             gCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCAcc                                                                      >  8:154716/1‑149 (MQ=255)
                                                                                   ccGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTcg                                                                >  7:82289/1‑149 (MQ=255)
                                                                                    cGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTcgg                                                               >  6:209429/1‑149 (MQ=255)
                                                                                        aTGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCggg                                                           <  8:82289/149‑1 (MQ=255)
                                                                                          ggAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCCCCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGt                                                         >  7:86355/1‑149 (MQ=255)
                                                                                                     ggA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGccc                                              <  8:238098/149‑1 (MQ=255)
                                                                                                      gA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCt                                             <  7:154716/149‑1 (MQ=255)
                                                                                                       a‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTg                                            >  4:213512/1‑149 (MQ=255)
                                                                                                            gTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACt                                        >  6:2680/1‑149 (MQ=255)
                                                                                                             ttCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTc                                       >  5:146947/1‑149 (MQ=255)
                                                                                                              tCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCg                                      >  2:95888/1‑149 (MQ=255)
                                                                                                                      tGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAggg                              >  2:299885/1‑149 (MQ=255)
                                                                                                                        gTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGcc                            >  7:127736/1‑149 (MQ=255)
                                                                                                                              gggcAACGGCGCGATGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCaaa                      >  5:20296/2‑149 (MQ=255)
                                                                                                                                 cAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTc                                                                 <  5:259198/103‑1 (MQ=255)
                                                                                                                                 cAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTc                                                                 >  6:259198/1‑103 (MQ=255)
                                                                                                                                 cAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATgg                   <  5:2680/149‑1 (MQ=255)
                                                                                                                                 cAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATgg                   <  1:95888/149‑1 (MQ=255)
                                                                                                                                   aCATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCAGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGc                 <  2:252462/149‑1 (MQ=255)
                                                                                                                                     aTCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCAt               >  3:74958/1‑149 (MQ=255)
                                                                                                                                      gcgcgCCGGTGACGATGGGCAACGGCGCGATGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGcc                                                                         >  2:262341/2‑90 (MQ=25)
                                                                                                                                      gcgcgCCGGTGACGATGGGCAACGGCGCGATGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGcc                                                                         <  1:262341/89‑1 (MQ=25)
                                                                                                                                             gTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGTGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAga       <  4:151251/149‑1 (MQ=255)
                                                                                                                                                  ccATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATc  >  3:254052/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCTTGGAGTCAGTCATGACCAGTGCAGCCAATTCGGCACCCCTGATAGAAAAACACACGATCGGCTACGTGCCCCCGCAAGACCGCCATGGAAAGGTAAGGGA‑TCTGTTCACACTCTGGTTCGGCGGCAACATCGCGCCGTTGCCCATCGTCACCGGCGCGCTGGGTGTGCAGCTGTTCCACCTGAACCTGGTATGGGGCATCGTCGCCATCCTTGTCGGCCACCTGGTCGGCGGGGTGCTGATGGCCCTGCACTCGGCCCAGGGCCCGCAAATGGGCATTCCGCAGATGATC  >  NC_002947/4150383‑4150676

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