Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,124 +G intergenic (+58/+26) 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,1231.G97.5% 181.1 / ‑1.3 81intergenic (+57/+27)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (35/44);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGG  >  NC_000913/3957051‑3957193
                                                                             |                                                                      
tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                       >  1:4414096/1‑76 (MQ=18)
 gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                      <  2:2305205/76‑2 (MQ=17)
 gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                      <  2:3756495/76‑2 (MQ=17)
 gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                      >  1:2133084/1‑75 (MQ=17)
  cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                                     <  1:4433576/76‑3 (MQ=17)
   tGTTGCGGGTTAAGGGCGCGCTGATGCCCTCACCCCGACCCGCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                                    >  1:703250/1‑73 (MQ=17)
      tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                  >  2:1941298/1‑76 (MQ=14)
      tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                  <  2:2925065/76‑1 (MQ=14)
       gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                 <  2:1612000/76‑1 (MQ=14)
       gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                 <  2:2368892/76‑1 (MQ=14)
       gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                 >  2:2722957/1‑76 (MQ=14)
       gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                 >  1:955500/1‑76 (MQ=14)
        cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                >  1:2074472/1‑76 (MQ=14)
        cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                <  1:2565759/76‑1 (MQ=14)
         gggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                               >  1:1461407/1‑76 (MQ=14)
         gggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                               >  1:2235355/1‑76 (MQ=14)
         gggTTAAGTGCGCGCTGATGACCTCACACCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                               >  1:901981/1‑76 (MQ=14)
          ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                              <  2:3257407/76‑1 (MQ=14)
          ggTTAAGTGCGCGCTGATGCCCGCACCCCGACCCTCTCCCACAGGGAGAGGGAGCACACACTCAAGGCCTTCTCCt                                                              <  2:703250/76‑1 (MQ=14)
             tAAGTGCGCGCTGATGCCCTCTCCCCGACCCTCTCCCACAGGTAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                           <  1:2589307/76‑1 (MQ=14)
               aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCCAGGCCTTCTCCTGgaga                                                         >  1:3318988/1‑76 (MQ=14)
                gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                        <  2:4635747/76‑1 (MQ=14)
                gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                        <  2:3931756/76‑1 (MQ=14)
                gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                        <  2:426167/76‑1 (MQ=14)
                 tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                       <  2:1289684/76‑1 (MQ=14)
                 tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                       <  1:1973819/76‑1 (MQ=17)
                 tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCCTCTCCTGGAGAAg                                                       >  2:990502/1‑76 (MQ=14)
                 ggcgcgcTGATGCCCTACCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                       <  2:2198243/75‑1 (MQ=14)
                   cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGCAGGc                                                     >  2:1746703/1‑76 (MQ=14)
                    gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAAGCCTGCTCCTGGAGAAGGcc                                                    >  1:883833/1‑76 (MQ=17)
                    gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                    >  1:3941187/1‑76 (MQ=17)
                    gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                    >  2:1309370/1‑76 (MQ=17)
                     cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                   >  1:1166489/1‑76 (MQ=17)
                      tcggATGCCCTGACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                  <  1:700137/73‑1 (MQ=17)
                      gcTGATGCCCTCACCGCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                  <  1:956227/76‑1 (MQ=17)
                         gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                               <  1:4145504/76‑1 (MQ=18)
                         gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                               <  2:832927/76‑1 (MQ=18)
                         gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                               >  1:1322609/1‑76 (MQ=18)
                          aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                              <  2:2774749/76‑1 (MQ=18)
                            gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                                            <  2:1682834/76‑1 (MQ=18)
                              ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                          >  1:3385060/1‑76 (MQ=21)
                               cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                         >  1:3049342/1‑76 (MQ=21)
                                tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                        >  2:4365379/1‑76 (MQ=21)
                                tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                        >  1:1918893/1‑76 (MQ=21)
                                ccACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                        <  2:684935/75‑1 (MQ=21)
                                 cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                       >  1:1400570/1‑76 (MQ=25)
                                        accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                <  1:1192640/73‑1 (MQ=38)
                                        accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                <  1:3965145/73‑1 (MQ=38)
                                        accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                >  1:2191565/4‑76 (MQ=38)
                                         ccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                               >  2:1669039/3‑76 (MQ=38)
                                         ccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                               <  1:2266195/74‑1 (MQ=38)
                                          cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                              >  2:3087169/2‑76 (MQ=255)
                                          cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                              <  2:1751870/75‑1 (MQ=255)
                                           ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGcc                             >  2:1483754/1‑75 (MQ=255)
                                            tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc                            <  1:2320381/76‑1 (MQ=255)
                                            gctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc                            <  2:2296007/75‑1 (MQ=255)
                                             ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                           >  2:3458080/1‑76 (MQ=255)
                                             ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                           <  1:1941299/76‑1 (MQ=255)
                                              tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt                          <  2:4224151/76‑1 (MQ=255)
                                                ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa                        >  2:1012290/1‑76 (MQ=255)
                                                  acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                      <  1:3849552/76‑1 (MQ=255)
                                                  acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                      >  2:1615497/1‑76 (MQ=255)
                                                    agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatga                    >  2:365944/1‑76 (MQ=255)
                                                    agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGGACAGCACCTTAatga                    >  2:726576/1‑76 (MQ=255)
                                                     gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                   <  2:1380563/76‑1 (MQ=255)
                                                      ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                  >  1:3444256/1‑76 (MQ=255)
                                                         gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa               <  2:4041260/76‑1 (MQ=255)
                                                          agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat              <  2:2319532/76‑1 (MQ=255)
                                                          agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat              <  2:2467262/76‑1 (MQ=255)
                                                          agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat              >  1:1962798/1‑76 (MQ=255)
                                                           gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata             >  1:1984821/1‑76 (MQ=255)
                                                           gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata             <  2:2877736/76‑1 (MQ=255)
                                                            ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat            <  2:3299150/76‑1 (MQ=255)
                                                             gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc           <  1:1509860/76‑1 (MQ=255)
                                                              agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc          <  2:2263968/76‑1 (MQ=255)
                                                              agaAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc          <  1:3748113/76‑1 (MQ=255)
                                                               gaAGACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg         <  1:2122368/76‑1 (MQ=255)
                                                               gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg         >  1:586364/1‑76 (MQ=255)
                                                                 aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa       <  1:1341982/76‑1 (MQ=255)
                                                                   acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACt     >  2:2441170/1‑76 (MQ=255)
                                                                    cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTg    <  2:1339943/76‑1 (MQ=255)
                                                                     acTCAAGGCCTTCTCCTGGCGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg   <  2:3280593/76‑1 (MQ=255)
                                                                     acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg   >  2:3154763/1‑76 (MQ=255)
                                                                     acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg   >  2:2475598/1‑76 (MQ=255)
                                                                     acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg   <  2:1918892/76‑1 (MQ=255)
                                                                      cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg  <  1:4814252/76‑1 (MQ=255)
                                                                      cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAAAGCACCTTAATGATGTGATATCCGAACTggg  >  2:182386/1‑76 (MQ=255)
                                                                             |                                                                      
TGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGG  >  NC_000913/3957051‑3957193

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