Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,553,315 A→T 13.6% L552Q (CTG→CAG)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,553,3150AT13.6% 118.0 / 10.9 44L552Q (CTG→CAG) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base A (14/24);  new base T (4/2);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

ATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCTGCAGCTGT  >  NC_002947/4553168‑4553462
                                                                                                                                                    |                                                                                                                                                   
aTTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAgg                                                                                                                                                   >  2:722567/1‑150 (MQ=255)
 ttCGCTACCCACAACTCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGc                                                                                                                                                  <  1:489120/150‑1 (MQ=255)
    gCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGa                                                                                                                                               <  2:596690/150‑1 (MQ=255)
     cTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAg                                                                                                                                              >  1:45869/1‑150 (MQ=255)
      tACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCTGCTGGCCGaac                                                                                                                                             <  2:105007/150‑3 (MQ=255)
           gcacagtgcCCGG‑CGCGCGACCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGtt                                                                                                                                        >  1:851134/8‑150 (MQ=255)
            caACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTg                                                                                                                                       <  1:527190/150‑1 (MQ=255)
                      ggcctcgct‑cCCGGCACTGGTCGCGGTTCTGACGCCTCGCGGTGTTGATGTTGCCTATGTGCCTGAAGGCTTTAGGTAATTTGCCGAGCATGGCCTGGGCGGCCGATGCTTACTGTTCTGCCAGCAGGCCGAGGTCGTTGCACAGGCCgg                                                                                                                             <  2:333929/145‑1 (MQ=255)
                              g‑cCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGgtgt                                                                                                                     <  2:45869/150‑1 (MQ=255)
                                 cGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTg                                                                                                                   <  2:693379/150‑1 (MQ=255)
                                            cgcgGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGAGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCtt                                                                                                        >  1:367792/1‑150 (MQ=255)
                                               ggTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGcc                                                                                                                                                 >  1:266891/1‑106 (MQ=255)
                                               ggTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGcc                                                                                                                                                 <  2:266891/106‑1 (MQ=255)
                                                gTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGc                                                                                                    >  1:294532/1‑150 (MQ=255)
                                                   cAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCttt                                                                                                 <  2:829966/150‑1 (MQ=255)
                                                         cAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTg                                                                                           <  2:294532/150‑1 (MQ=255)
                                                           gcgcGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGAACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGgc                                                                                         <  1:412493/150‑1 (MQ=255)
                                                              cgGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGgcgca                                                                                      <  2:669301/150‑2 (MQ=255)
                                                                   ttGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCtt                                                                                                  >  1:374731/1‑133 (MQ=255)
                                                                   ttGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCtt                                                                                                  <  2:374731/133‑1 (MQ=255)
                                                                     gatgatGCCTATGTGGCTGAAGGCGTGGGGGAAGTTGACGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGTCAGTAGGCCGAGGTCGTTGCACAGGCAGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTaa                                                                               <  2:139276/150‑1 (MQ=255)
                                                                     gatgatGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTaa                                                                               <  2:485695/150‑1 (MQ=255)
                                                                          tgCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACAt                                                                          >  2:580676/1‑150 (MQ=255)
                                                                                tCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGt                                                                    >  2:41617/1‑150 (MQ=255)
                                                                                  tGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca                                                                  <  2:204693/150‑1 (MQ=255)
                                                                                     cTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagcc                                                               >  2:48659/1‑150 (MQ=255)
                                                                                      tGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagcca                                                              <  1:722567/150‑1 (MQ=255)
                                                                                                        tGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCa                                            <  1:827619/150‑1 (MQ=255)
                                                                                                           cGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGa                                         >  1:2079/1‑150 (MQ=255)
                                                                                                             aGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAg                                       <  1:48659/150‑1 (MQ=255)
                                                                                                                 tGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGc                                   <  1:580676/150‑1 (MQ=255)
                                                                                                                    gcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagc                                                                >  2:588432/1‑118 (MQ=255)
                                                                                                                    gcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagc                                                                <  1:588432/118‑1 (MQ=255)
                                                                                                                            gccggGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGc                        <  1:250004/150‑1 (MQ=255)
                                                                                                                                   cGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGc                 <  2:367792/150‑1 (MQ=255)
                                                                                                                                   cGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGc                 <  1:676142/150‑1 (MQ=255)
                                                                                                                                      aCTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca                                                                  <  1:669306/98‑1 (MQ=255)
                                                                                                                                      aCTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca                                                                  >  2:669306/1‑98 (MQ=255)
                                                                                                                                      aCTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGt              >  1:314190/1‑150 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca                                                                  >  2:666487/2‑94 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca                                                                  <  1:666487/93‑1 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCt          >  2:91522/2‑150 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGCTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACCCCAGGACGGTGCCCTCTCCCGGCGTCACGccctccct          >  2:344529/2‑150 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCACGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCt          >  1:509129/2‑150 (MQ=255)
                                                                                                                                                  gcagGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCTGCAGctgt  <  2:851134/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCTGCAGCTGT  >  NC_002947/4553168‑4553462

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