Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,447,461 G→C 12.5% L269L (CTG→CTC nagZ → beta‑N‑acetylglucosaminidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,447,4610GC12.5% 112.9 / 8.9 40L269L (CTG→CTCnagZbeta‑N‑acetylglucosaminidase
Reads supporting (aligned to +/- strand):  ref base G (16/19);  new base C (3/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

ACCAGCCCGCTGGCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAACCGATTATCGC  >  NC_002947/2447316‑2447605
                                                                                                                                                 |                                                                                                                                                
aCCAGCCCGCTGGCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAgcg                                                                                                                                              >  2:637738/1‑150 (MQ=255)
   aGCCCGCTGGCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCg                                                                                                                                           <  2:188725/150‑1 (MQ=255)
   aGCCCGCTGGCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGGGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCg                                                                                                                                           >  2:714731/1‑150 (MQ=255)
     cccGCTGGCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCAGTATCGAGGCGGCGCTCAGCGCCgcc                                                                                                                                         <  1:204091/150‑3 (MQ=255)
        gCTGGCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTtgt                                                                                                                                      <  2:973845/150‑1 (MQ=255)
            gCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGata                                                                                                                                  <  2:892646/150‑1 (MQ=255)
              ttCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATg                                                                                                                                <  2:145904/150‑1 (MQ=255)
               tctcGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATgg                                                                                                                               <  2:801126/150‑1 (MQ=255)
                 tcGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGc                                                                                                                             <  1:223672/150‑1 (MQ=255)
                   gcggcgTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCt                                                                                                                           >  1:240555/1‑150 (MQ=255)
                            gttgCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCaa                                                                                                                  >  1:1067131/1‑150 (MQ=255)
                             ttgCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAAt                                                                                                                 >  1:1064877/1‑150 (MQ=255)
                             ttgCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGAGCTGAGCGCAGGGTGTTATATGGGCCTGGTGTGCAAt                                                                                                                 >  2:289187/1‑150 (MQ=255)
                                        ccTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGggcggc                                                                                                      <  1:937841/150‑1 (MQ=255)
                                                   aGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTgg                                                                                           <  2:1012470/150‑1 (MQ=255)
                                                    gCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGgc                                                                                          <  2:547042/150‑1 (MQ=255)
                                                         gCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGGCGCGGCCAGCCGTATCGAGGCGGCGCTGGGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGAACGGGCGGCGGCGGAGCTGGCGCTGa                                                                                     <  2:196173/150‑1 (MQ=255)
                                                           tttGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGt                                                                                   >  1:255813/1‑150 (MQ=255)
                                                            ttGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTg                                                                                  <  2:1067131/150‑1 (MQ=255)
                                                               aTGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTgcgg                                                                               >  2:494894/1‑150 (MQ=255)
                                                                                  tgaCCTGTCGATGGCCGGTGCGCATGTGGTGGGCGCCGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCaa                                                            <  1:289187/150‑1 (MQ=255)
                                                                                  tgaCCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCaa                                                            <  2:588645/150‑1 (MQ=255)
                                                                                    aCCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGACGGTCAAGc                                                          >  2:562990/1‑150 (MQ=255)
                                                                                     ccTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGcc                                                         <  2:255813/150‑1 (MQ=255)
                                                                                                ccGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGa                                              <  1:723752/150‑1 (MQ=255)
                                                                                                      cgcATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTgcgc                                        <  1:637738/150‑1 (MQ=255)
                                                                                                                    gACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCaa                          >  2:853730/1‑150 (MQ=255)
                                                                                                                            cAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTgc                  >  1:1022287/1‑150 (MQ=255)
                                                                                                                                  tATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAAc            >  2:144502/1‑150 (MQ=255)
                                                                                                                                      gAGGCGGCGCTGAGCGCCGGTAGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGGGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGCATTGCGCGGATGCGTTGGCAAGGGGTTGCGCGAACCGAt        >  2:398662/1‑150 (MQ=255)
                                                                                                                                       aaccggcgCTCAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAACCGAtt       >  1:802882/4‑150 (MQ=255)
                                                                                                                                        accggcgCTCAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAACCGATTa      >  2:829665/3‑150 (MQ=255)
                                                                                                                                         gcggcgCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCTGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGAtt                                            <  2:1065400/111‑1 (MQ=255)
                                                                                                                                         gcggcgCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCTGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGAtt                                            >  1:1065400/1‑111 (MQ=255)
                                                                                                                                         gcggcgCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAACCGATTAt     >  1:836344/1‑150 (MQ=255)
                                                                                                                                         ccggcgCTCAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAg                                                                 >  1:890240/2‑90 (MQ=255)
                                                                                                                                         ccggcgCTCAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAg                                                                 <  2:890240/89‑1 (MQ=255)
                                                                                                                                           ggcgCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGgcgcagc                                                                         >  1:875494/1‑80 (MQ=255)
                                                                                                                                           ggcgCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGgcgcagc                                                                         <  2:875494/80‑1 (MQ=255)
                                                                                                                                            gcgcTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAACCGATTATCGc  <  2:240555/150‑1 (MQ=255)
                                                                                                                                                 gAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAACCGa         <  1:581037/138‑1 (MQ=255)
                                                                                                                                                 gAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGAAAGGGTTTGCGCGAACCGa         >  2:581037/1‑138 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACCAGCCCGCTGGCTTCTCGCGGCGTTGGTTGCAGGACATCCTGCGTGGCGAGCTGGGCTTTGATGGGGTGATTTTCAGTGATGACCTGTCGATGGCCGGTGCGCATGTGGTGGGCGACGCGGCCAGCCGTATCGAGGCGGCGCTGAGCGCCGGTTGTGATATGGGCCTGGTGTGCAATGACCGGGCGGCGGCGGAGCTGGCGCTGAGTGCGGCGCAGCGCATGAAGGTCAAGCCATCGCCGCGGATTGCGCGGATGCGTGGGCAAGGGTTTGCGCGAACCGATTATCGC  >  NC_002947/2447316‑2447605

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