Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,297,638 G→A 10.2% R154C (CGC→TGC)  PP_3768 ← shikimate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,297,6380GA10.2% 129.5 / 6.2 49R154C (CGC→TGC) PP_3768shikimate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base A (2/3);  total (22/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 = 7.43e-01

TCGGGAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGGTG  >  NC_002947/4297490‑4297786
                                                                                                                                                    |                                                                                                                                                    
tCGGGAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCACATGCCGCGCGTCTCAAGCACCCAGCGATAACCACACACATCCGCCCCCTGACGCCCGGCCTGCTCGTTGCGt                                                                                                                                                     >  1:216026/1‑150 (MQ=255)
tCGGGAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGt                                                                                                                                                     <  2:390458/150‑1 (MQ=255)
           cAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTAACGCCCGGCCTGCTCGTTGCGTGCAACgatgat                                                                                                                                          <  2:541011/150‑1 (MQ=255)
              ctctTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagcaggcc                                                                                                                                       <  2:62736/150‑8 (MQ=255)
               tcttctGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCCCCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgat                                                                                                                                             <  1:496072/143‑1 (MQ=255)
               tcttctGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgat                                                                                                                                             >  2:496072/1‑143 (MQ=255)
                       cTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgatgat                                                                                                                                          >  1:710187/1‑138 (MQ=255)
                       cTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgatgat                                                                                                                                          <  2:710187/138‑1 (MQ=255)
                          cAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAAccc                                                                                                                           >  1:492876/1‑150 (MQ=255)
                              ccGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgat                                                                                                                                             >  1:38076/1‑128 (MQ=255)
                              ccGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgat                                                                                                                                             <  2:38076/128‑1 (MQ=255)
                                    gCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCgcgc                                                                                                                 <  2:210932/150‑1 (MQ=255)
                                     ccATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCa                                                                                                                >  2:289560/1‑150 (MQ=255)
                                      cATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCaa                                                                                                               >  2:76470/1‑150 (MQ=255)
                                      cATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCaa                                                                                                               <  2:415850/150‑1 (MQ=255)
                                       aTACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAg                                                                                                              >  2:70928/1‑150 (MQ=255)
                                            gATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCg                                                                                                                                   <  1:660195/124‑1 (MQ=255)
                                            gATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCg                                                                                                                                   >  2:660195/1‑124 (MQ=255)
                                                 gCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagca                                                                                                                                           >  2:111338/1‑108 (MQ=255)
                                                 gCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagca                                                                                                                                           <  1:111338/111‑4 (MQ=255)
                                                        gTTCCCCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGCCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCtt                                                                                             <  1:564249/150‑1 (MQ=255)
                                                        gTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCtt                                                                                             >  1:364575/1‑150 (MQ=255)
                                                                gCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGc                                                                                     >  1:674473/1‑150 (MQ=255)
                                                                  aTGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAGCTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCa                                                                                   >  1:635267/1‑150 (MQ=255)
                                                                    gggcgggcAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCAcc                                                                                 >  1:148349/1‑150 (MQ=255)
                                                                      gcgggcAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGc                                                                               >  1:517561/1‑150 (MQ=255)
                                                                                gTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAg                                                                     <  2:110344/150‑1 (MQ=255)
                                                                                gTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAg                                                                     <  2:148349/150‑1 (MQ=255)
                                                                                    ccTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgag                                                                                                                                             <  1:631684/74‑2 (MQ=255)
                                                                                    ccTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgag                                                                                                                                             >  2:631684/1‑73 (MQ=255)
                                                                                         cTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGc                                                            <  1:289560/150‑1 (MQ=255)
                                                                                             ggCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGTAAAGGCCGTg                                                        <  2:517561/150‑1 (MQ=255)
                                                                                                   cAGCGATAACCACAGACATCCGCCAACTGGCGCCCGGCCTGCTCGTTGCGGGCAACGATGATGGCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGGCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGa                                                  <  1:76196/150‑1 (MQ=255)
                                                                                                   cAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCgg                                                                                                                         <  2:296043/79‑1 (MQ=255)
                                                                                                   cAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCgg                                                                                                                         >  1:296043/1‑79 (MQ=255)
                                                                                                        aTAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGAc                                             <  2:364575/150‑1 (MQ=255)
                                                                                                        aTAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGAc                                             <  1:748549/150‑1 (MQ=255)
                                                                                                            ccACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGAt                                                 >  2:333677/1‑142 (MQ=255)
                                                                                                            ccACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGAt                                                 <  1:333677/142‑1 (MQ=255)
                                                                                                            ccACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGa                                         >  2:315795/1‑150 (MQ=255)
                                                                                                                  aCATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGa                                   <  2:305564/150‑1 (MQ=255)
                                                                                                                        gCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGAc                                             >  1:618327/1‑134 (MQ=255)
                                                                                                                        gCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGAc                                             <  2:618327/134‑1 (MQ=255)
                                                                                                                            aCTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAAtt                         >  1:442715/1‑150 (MQ=255)
                                                                                                                              tGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAACGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTgg                       <  1:499093/150‑1 (MQ=255)
                                                                                                                                        ccTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCa             <  1:244926/150‑1 (MQ=255)
                                                                                                                                         atcatCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGa                                                  >  1:380799/5‑112 (MQ=255)
                                                                                                                                         atcatCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGa                                                  <  2:380799/108‑1 (MQ=255)
                                                                                                                                                   cGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGgtg  <  1:569761/150‑1 (MQ=255)
                                                                                                                                                    gTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGAt                                        >  1:46832/1‑111 (MQ=255)
                                                                                                                                                    gTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGAt                                        <  2:46832/111‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGGGAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGGTG  >  NC_002947/4297490‑4297786

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