Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,902,856 G→A 18.7% intergenic (+30/+39) quiC → / ← PP_2555 3‑dehydroshikimate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,902,8560GA18.7% 65.9 / 8.3 32intergenic (+30/+39)quiC/PP_25553‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/14);  new base A (2/4);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

CGAGATCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTGCAGACTGACCAGCCTCAACCCCGCCATGTCCAGCGCA  >  NC_002947/2902712‑2902997
                                                                                                                                                |                                                                                                                                             
cGAGATCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTggg                                                                                                                                             <  4:251632/147‑1 (MQ=255)
cGAGATCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTggg                                                                                                                                             <  4:159012/147‑1 (MQ=255)
 gagaTCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGTCTGTCTTGTGGGAGTGGGGt                                                                                                                                           <  6:330008/148‑1 (MQ=255)
   gaTCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGTGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGt                                                                                                                                           <  1:4241/146‑1 (MQ=255)
   gaTCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGt                                                                                                                                           >  2:4241/1‑146 (MQ=255)
   gaTCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGAAGCAGGGttca                                                                                                                                        >  4:259929/1‑146 (MQ=255)
     tcatcCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAAcc                                                                                                                                      <  2:57497/149‑1 (MQ=255)
        tcCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTg                                                                                                                                   >  6:360378/1‑149 (MQ=255)
         ccAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGc                                                                                                                                  >  1:106872/1‑149 (MQ=255)
          cAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCt                                                                                                                                 <  4:208239/149‑1 (MQ=255)
               cAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGAAGCAGGGTTCACCCCACTCCcac                                                                                                                            <  7:250308/149‑1 (MQ=255)
                aaGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCcaca                                                                                                                           <  6:252904/149‑1 (MQ=255)
                 aGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAg                                                                                                                          <  5:360378/149‑1 (MQ=255)
                 aGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGAAGCAGGGTTCACCCCACTCCcacaa                                                                                                                          <  3:259929/149‑2 (MQ=255)
                   gCGGGTTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGTGTGAACCCTGCTTCCACAGGa                                                                                                                        <  6:355142/149‑1 (MQ=255)
                       ggTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGattat                                                                                                                    >  4:123158/1‑149 (MQ=255)
                         tACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTAGGGTGAACCCTGCTTCCACAGGAttatta                                                                                                                  <  7:302587/149‑1 (MQ=255)
                                 ttACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGa                                                                                                          <  2:106872/149‑1 (MQ=255)
                                                  ttGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGccc                                                                                                            >  7:292927/1‑130 (MQ=255)
                                                  ttGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGccc                                                                                                            <  8:292927/130‑1 (MQ=255)
                                                  ttGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGc                                                                                          >  5:347623/1‑148 (MQ=255)
                                                                    ccATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGAAGCAGGGTTCACCCCACTCCcacaa                                                                                                                          <  5:250336/98‑2 (MQ=25)
                                                                    ccATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGAAGCAGGGTTCACCCCACTCCcacaa                                                                                                                          >  6:250336/1‑97 (MQ=25)
                                                                                         gggCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGAtgc                                                  <  3:123158/149‑1 (MQ=255)
                                                                                                  gcAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTcc                                         >  8:31044/1‑149 (MQ=255)
                                                                                                  gcAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTcc                                         >  2:16891/1‑149 (MQ=255)
                                                                                                      gCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTGCa                                     <  6:347623/149‑1 (MQ=255)
                                                                                                       ccGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTg                                       >  6:112054/1‑146 (MQ=255)
                                                                                                        cGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTGCAGa                                   >  8:2365/1‑149 (MQ=255)
                                                                                                                tAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTGCAGACTGAcc                             <  7:31044/147‑1 (MQ=255)
                                                                                                                         gCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCAcc                                                                 >  2:263170/1‑102 (MQ=255)
                                                                                                                         gCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCAcc                                                                 <  1:263170/102‑1 (MQ=255)
                                                                                                                             cGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTGCAGACTGACCAGCCTCAACCCCGc               >  2:234398/1‑148 (MQ=255)
                                                                                                                                         tgGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTGCAGACTGACCAGCCTCAACCCCGCCATGTCCAGCGCa  <  1:16891/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CGAGATCATCCAGCGCAAGGCGGGGTACGCTGGTTACGGCGCTGCCAACGTTGCGGTGCGCCTGGCAGCCATGGCCAAGGCCCGTAGCGGGGCGGCGCGCAAGCCGGTGCTGTAGTCGCCTGCACCGGCCTGTCTTGTGGGAGTGGGGTGAACCCTGCTTCCACAGGATTATTACGGCCCGATTCAGTGACGCTCAGCCACCATCAACAACGACTGCGGCACCGCGCTCTGCGGATGCTGCGGTTCCTGCAGACTGACCAGCCTCAACCCCGCCATGTCCAGCGCA  >  NC_002947/2902712‑2902997

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: