Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 349,785 A→G intergenic (+213/‑227) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,7850AG97.6% 101.6 / ‑4.9 43intergenic (+213/‑227)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (21/20);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTAC  >  NC_000913/349652‑349919
                                                                                                                                     |                                                                                                                                      
aaCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATg                                                                                                                                  >  2:596550/1‑140 (MQ=38)
     tAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGa                                                                                                                             >  2:123437/1‑140 (MQ=255)
         tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATgcg                                                                                                                         <  2:786957/140‑1 (MQ=37)
          cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATgcgc                                                                                                                        >  1:142864/1‑140 (MQ=255)
          cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATgcgc                                                                                                                        >  1:476651/1‑140 (MQ=255)
            gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGcgccg                                                                                                                      >  1:186835/1‑140 (MQ=255)
               aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTg                                                                                                                   >  2:100126/1‑140 (MQ=255)
                     gATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGa                                                                                                             <  1:124321/140‑1 (MQ=25)
                          cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                                        <  1:385647/140‑1 (MQ=17)
                          cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                                        <  1:165625/140‑1 (MQ=17)
                          cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                                        <  2:186835/140‑1 (MQ=17)
                          cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                                        <  1:123437/140‑1 (MQ=17)
                                aTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc                                                                                                  >  1:428776/1‑140 (MQ=17)
                                 tCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCt                                                                                                 <  1:786592/140‑1 (MQ=17)
                                    gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACg                                                                                              <  2:751836/140‑1 (MQ=17)
                                           gATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCa                                                                                       <  2:476651/140‑1 (MQ=18)
                                           gATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCa                                                                                       <  2:142864/140‑1 (MQ=18)
                                            gtGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCAc                                                                                      <  1:696614/139‑1 (MQ=21)
                                                   aTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACg                                                                                              >  1:399175/1‑125 (MQ=18)
                                                   aTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACg                                                                                              <  2:399175/125‑1 (MQ=17)
                                                    tGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                              >  2:667418/1‑140 (MQ=25)
                                                        ccgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAggc                                                                          >  2:696614/1‑139 (MQ=21)
                                                             gACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                                     <  1:596550/140‑1 (MQ=25)
                                                             gACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                                     <  1:667418/140‑1 (MQ=25)
                                                                           ggCCTACGGGGTACACTGAACTGTAGGTCGGATAAGTCGCATAGCGTCGCATCCGACATCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                                       <  2:726223/140‑1 (MQ=12)
                                                                            gCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                      >  1:616440/1‑140 (MQ=17)
                                                                               tACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                                                   <  2:236230/140‑1 (MQ=32)
                                                                                 cGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc                                                 >  1:743107/1‑140 (MQ=17)
                                                                                   gggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                               <  2:772056/140‑1 (MQ=17)
                                                                                    ggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                              <  2:616440/140‑1 (MQ=17)
                                                                                      tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACa                                            >  2:206754/1‑140 (MQ=17)
                                                                                      tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACa                                            >  1:685142/1‑140 (MQ=17)
                                                                                           tGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGa                                       >  1:722393/1‑140 (MQ=21)
                                                                                            gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                      >  1:254200/1‑140 (MQ=25)
                                                                                                               aCGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTCTCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcg                   <  2:722393/140‑1 (MQ=255)
                                                                                                                cGCATAGCGTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgc                  >  2:579704/1‑140 (MQ=37)
                                                                                                                        gTCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc          >  1:189772/1‑140 (MQ=17)
                                                                                                                         tCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa         <  2:575334/140‑1 (MQ=17)
                                                                                                                         tCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa         >  2:517812/1‑140 (MQ=17)
                                                                                                                         tCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa         >  1:575334/1‑140 (MQ=17)
                                                                                                                          cGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg        >  1:236230/1‑140 (MQ=17)
                                                                                                                              tCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCt    <  2:749200/140‑1 (MQ=17)
                                                                                                                                cGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTAc  <  2:562672/140‑1 (MQ=17)
                                                                                                                                     |                                                                                                                                      
AACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTAC  >  NC_000913/349652‑349919

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