Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,092 +CC intergenic (+26/+58) ilvC → / ← ppiC ketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,0911.C98.4% 159.4 / ‑3.3 62intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (3/58);  total (3/59)
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.15e-01
*NC_0009133,957,0912.C98.4% 159.7 / ‑3.3 62intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (3/58);  total (3/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAA  >  NC_000913/3956999‑3957172
                                                                                               ||                                                                                  
ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑a                                                                                      <  2:2140753/93‑1 (MQ=255)
 tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑c                                                                                   <  1:2172536/93‑1 (MQ=255)
      aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctct                                                                              <  1:3638144/93‑2 (MQ=255)
            gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGCCC‑‑CTCTCCcaca                                                                          >  1:686502/1‑93 (MQ=37)
             taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACag                                                                         <  2:1989207/93‑1 (MQ=37)
              aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagg                                                                        <  2:2791016/93‑1 (MQ=37)
              aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagg                                                                        <  2:637054/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACaggga                                                                      <  1:929194/93‑1 (MQ=37)
                aaaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACaggga                                                                      <  1:3612524/92‑1 (MQ=37)
                  aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggaga                                                                    <  2:2463948/93‑1 (MQ=37)
                   aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagag                                                                   <  1:2980668/93‑1 (MQ=37)
                      gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggga                                                                <  1:1117799/93‑1 (MQ=25)
                      gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggga                                                                <  2:4466141/93‑1 (MQ=25)
                      gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggga                                                                <  2:4537271/93‑1 (MQ=25)
                          ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaa                                                            <  2:2479171/93‑1 (MQ=25)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaaa                                                           <  1:3318265/93‑1 (MQ=21)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaaa                                                           <  1:3318266/93‑1 (MQ=21)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaaa                                                           <  1:3726013/93‑1 (MQ=21)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaaa                                                           <  1:4106262/93‑1 (MQ=21)
                            cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAac                                                          <  2:342729/93‑1 (MQ=21)
                            cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAac                                                          <  1:1496496/93‑1 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCCCAGGGAGGGGGGGAAAaca                                                         >  2:4731513/1‑93 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAaca                                                         <  1:3315249/93‑1 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAAGGAGAGGGAGaaaaaa                                                         >  2:796281/1‑91 (MQ=21)
                              atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                        <  1:3560124/93‑1 (MQ=21)
                              atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                        <  2:542502/93‑1 (MQ=21)
                               tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCCCTCCCACAGGGAGAGGGGGaaaaacac                                                       >  2:4695305/1‑89 (MQ=21)
                                 tGACAGATATGAAACTTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                     <  1:320225/93‑1 (MQ=21)
                                  gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                    <  1:1814825/93‑1 (MQ=21)
                                   aaaGATATGAATGTTATTGCTGTTGCTGGTTAAGTGCGCGCTGAGGCCCTCCCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                  <  1:1324441/91‑1 (MQ=17)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                  <  2:2121640/93‑1 (MQ=21)
                                    cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                 <  2:2757683/93‑2 (MQ=21)
                                     aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACCCAAggc                                                <  2:3588179/93‑3 (MQ=21)
                                      gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                               <  1:2876473/93‑4 (MQ=21)
                                      gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                               <  1:3281217/93‑4 (MQ=21)
                                       atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                              <  1:98335/93‑5 (MQ=21)
                                         atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                             <  2:1026483/93‑1 (MQ=18)
                                           gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                           <  1:1789410/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                          <  1:2337875/93‑1 (MQ=18)
                                              aCGTATTGCGGTGGGGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                        <  2:1217104/93‑1 (MQ=17)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                       <  2:3637162/93‑1 (MQ=18)
                                                gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                      <  1:461609/93‑1 (MQ=21)
                                                gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCACCCCG‑‑ACCCTCTCCCACACGGGGAGGGGGAAAAAACTCAAGGACTTCTTCTGga                                      >  2:3231444/1‑93 (MQ=12)
                                                     gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                 <  2:1006014/93‑1 (MQ=18)
                                                     gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                 <  2:2113099/93‑1 (MQ=18)
                                                     gCTGTTGCGCGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                 <  2:645493/93‑1 (MQ=18)
                                                       tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                               <  2:3372444/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                              <  2:1675889/93‑1 (MQ=18)
                                                         ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGATAGGTAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                             <  1:1619048/93‑1 (MQ=18)
                                                          tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                            <  1:1652904/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                           <  2:2292213/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                           <  2:649743/93‑1 (MQ=18)
                                                             gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                         <  1:494222/93‑1 (MQ=18)
                                                               gTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                       <  1:2546902/93‑1 (MQ=18)
                                                                ttAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                      <  2:111654/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                    <  1:581871/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                    <  1:4531956/93‑1 (MQ=18)
                                                                      gcgcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                <  1:1301489/93‑1 (MQ=21)
                                                                      gcgcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                <  1:4549574/93‑1 (MQ=21)
                                                                       cgcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg               <  1:1987246/93‑1 (MQ=21)
                                                                       cgcgcTGATGCCCGCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg               <  1:1987245/93‑1 (MQ=21)
                                                                        gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt              <  2:4723522/93‑1 (MQ=21)
                                                                         cgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa             <  1:4013827/93‑1 (MQ=21)
                                                                               tGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc       <  2:813318/93‑1 (MQ=25)
                                                                                gCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc      <  1:980248/93‑1 (MQ=25)
                                                                                 cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt     <  2:3024969/93‑1 (MQ=25)
                                                                                  ccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt    >  1:1652386/1‑93 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGGGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa  >  2:3754238/1‑93 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa  <  1:444171/93‑1 (MQ=37)
                                                                                               ||                                                                                  
TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAA  >  NC_000913/3956999‑3957172

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