Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,808,075 G→A 9.2% V42V (GTG→GTA kdsA‑I → 2‑dehydro‑3‑deoxyphosphooctonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,808,0750GA 9.2% 148.8 / 3.2 54V42V (GTG→GTAkdsA‑I2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (28/21);  new base A (2/3);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAGAACAAGGCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGG  >  NC_002947/1807930‑1808223
                                                                                                                                                 |                                                                                                                                                    
agaagaACAAGGCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGc                                                                                                                                                     <  3:114366/147‑1 (MQ=255)
 gaagaaCAAGGCCTGAACGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCggg                                                                                                                                                  >  7:84249/1‑149 (MQ=255)
    gaaCAAGGCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATa                                                                                                                                               >  1:210001/1‑149 (MQ=255)
     aaCAAGGCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTAc                                                                                                                                              >  7:106953/1‑149 (MQ=255)
         aGGCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCTGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCgag                                                                                                                                          <  7:106775/149‑1 (MQ=255)
         aGGCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCgag                                                                                                                                          <  8:22829/149‑1 (MQ=255)
           gCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGaa                                                                                                                                        <  8:170679/149‑1 (MQ=255)
                    aTGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTAt                                                                                                                               >  8:255212/1‑149 (MQ=255)
                        cccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCg                                                                                                                           >  2:117188/1‑149 (MQ=255)
                           agaagaTCATTCGCGTCGGTAACATCGGGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTAc                                                                                                                        <  5:278092/149‑1 (MQ=255)
                           agaagaTCATTCGCGTCGGTAACATCGAGATCGCCACCGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTAc                                                                                                                        >  1:243735/1‑149 (MQ=255)
                           agaagaTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCAc                                                                                                                        <  2:86237/149‑1 (MQ=39)
                              agaTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGGCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGTAACCCGCACgta                                                                                                                     <  1:40670/149‑2 (MQ=38)
                                   aTTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCa                                                                                                                 <  1:117188/148‑1 (MQ=255)
                                    ttCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGaa                                                                                                                                        <  6:87299/124‑1 (MQ=255)
                                    ttCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGaa                                                                                                                                        >  5:87299/1‑124 (MQ=255)
                                          tCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGAAGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGt                                                                                                         >  8:52975/1‑149 (MQ=255)
                                             gTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                      <  2:48431/149‑1 (MQ=255)
                                                 cATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGAGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACATGTTCAAGGCCAGTTTCGAc                                                                                                   >  7:61420/1‑148 (MQ=255)
                                                   tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGAc                                                                                                   >  2:306480/1‑146 (MQ=255)
                                                   tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACaa                                                                                                 <  1:306480/148‑1 (MQ=255)
                                                   tCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAg                                                                                                <  7:255212/149‑1 (MQ=255)
                                                     gagaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGt                                                                                                                                 <  2:263788/114‑1 (MQ=255)
                                                     gagaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTACGGGATACCGAGCTTCTCGGt                                                                                                                                 >  1:263788/1‑114 (MQ=255)
                                                          cGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTAACGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAAcc                                                                                         >  6:228802/1‑149 (MQ=255)
                                                              aaCGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtc                                                                                     >  2:13277/1‑149 (MQ=255)
                                                                    aaGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCgg                                                                                >  7:180056/1‑148 (MQ=255)
                                                                               ccTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTAcc                                                                    >  7:209130/1‑149 (MQ=255)
                                                                                cTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCg                                                                   >  5:207540/1‑149 (MQ=255)
                                                                                 tGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGt                                                                  >  5:222616/1‑149 (MQ=255)
                                                                                 tGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGCGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGt                                                                  >  2:84596/1‑149 (MQ=255)
                                                                                  gTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGANGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTg                                                                 >  1:96411/1‑149 (MQ=255)
                                                                                     cggcggCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGcc                                                              >  2:291778/1‑149 (MQ=255)
                                                                                       gcggcATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGACTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGATTCCACAAGGCCAACCGTTCGTCGGAAACGTCCTACCATGGCCCt                                                            >  3:44909/1‑149 (MQ=255)
                                                                                            aTGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCAt                                                       >  5:94116/1‑149 (MQ=255)
                                                                                              gAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATg                                                      >  1:127662/1‑148 (MQ=255)
                                                                                               aaCGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATgg                                                     <  8:209130/148‑1 (MQ=255)
                                                                                                aCGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATgg                                                     >  4:40341/1‑147 (MQ=255)
                                                                                                aCGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGaa                                                   <  1:291778/149‑1 (MQ=255)
                                                                                                      tGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAgggg                                             >  3:202484/1‑149 (MQ=255)
                                                                                                        gAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAgggg                                             <  6:222616/147‑1 (MQ=255)
                                                                                                          aTCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGa                                         <  6:207540/149‑1 (MQ=255)
                                                                                                             ccGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTg                                          >  6:170268/1‑145 (MQ=255)
                                                                                                                  aCCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTc                                 <  8:106953/149‑1 (MQ=255)
                                                                                                                    cTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGa                               <  5:170268/149‑1 (MQ=255)
                                                                                                                       gCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                              <  4:202484/147‑1 (MQ=255)
                                                                                                                            gAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATc                        >  4:283493/1‑148 (MQ=255)
                                                                                                                               ggTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAgc                    >  7:163970/1‑149 (MQ=255)
                                                                                                                               ggTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAgc                    <  2:127662/149‑1 (MQ=255)
                                                                                                                               ggTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCCTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAgc                    >  2:153356/1‑149 (MQ=255)
                                                                                                                                     cgacgaGTACGTGCGGGTTCCCGAGAAGCGCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCAcc               <  6:156398/148‑1 (MQ=255)
                                                                                                                                     cGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTACCGTCCTACCGTGGCCCGGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCAcc               >  5:156398/1‑148 (MQ=255)
                                                                                                                                     cGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCt              <  1:13277/149‑1 (MQ=255)
                                                                                                                                           gTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACg        <  3:40341/149‑1 (MQ=255)
                                                                                                                                                 gCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCgg  <  2:96411/149‑1 (MQ=255)
                                                                                                                                                 gCGGGTTACCGAGAAGCACGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCAAACCGTTCGTAGGTAACGTCCTACCGTGGCCCTGGCATCGAAGAGGGGCTGAAGAACTTCGAAGAGATCAAGCGCACCTTCAACGTGCCgg  >  3:255871/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGAAGAACAAGGCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGG  >  NC_002947/1807930‑1808223

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: