Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,245,543 T→A 15.3% intergenic (‑59/‑36) nudC ← / → rsd NADH pyrophosphatase/stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,245,5430TA15.3% 215.3 / 27.6 104intergenic (‑59/‑36)nudC/rsdNADH pyrophosphatase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base T (35/53);  new base A (3/13);  total (38/66)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 104 total aligned reads displayed.

GCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGG  >  minE/2245486‑2245606
                                                         |                                                               
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTTAGCAGtttt                                                       >  1:910216/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:1128930/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:1440277/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:19528/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:820616/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:642375/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:1109361/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:1108156/1‑68 (MQ=255)
gCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:1080087/1‑68 (MQ=255)
  tctTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTg                                                     <  1:1092909/68‑1 (MQ=255)
  tctTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTg                                                     <  1:1857942/68‑1 (MQ=255)
  tctTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTg                                                     <  1:2827938/68‑1 (MQ=255)
  tctTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTg                                                     <  1:3087882/68‑1 (MQ=255)
  tctTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTg                                                     <  1:3127431/68‑1 (MQ=255)
  tctTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTg                                                     <  1:1583139/68‑1 (MQ=255)
      gCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGtttt                                                       >  1:2065381/1‑62 (MQ=255)
      gCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATa                                                 >  1:2725833/1‑68 (MQ=255)
      gCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATa                                                 >  1:77835/1‑68 (MQ=255)
      gCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATa                                                 >  1:2043370/1‑68 (MQ=255)
         cTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:1325424/68‑1 (MQ=255)
         cTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:982603/68‑1 (MQ=255)
         cTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:414266/68‑1 (MQ=255)
         cTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:44812/68‑1 (MQ=255)
         cTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:100790/68‑1 (MQ=255)
           aCCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACt                                           <  1:1542423/69‑1 (MQ=255)
           aCCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACt                                           <  1:1542327/69‑1 (MQ=255)
           aCCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACt                                           <  1:717619/69‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGa                                     >  1:518660/1‑68 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:206747/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:3111142/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:772938/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:2051754/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:1291451/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:128544/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:2588786/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:2676885/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:890068/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:2025842/68‑1 (MQ=255)
                  cATCACTGGCATGTTTAACATGGTTTTTACATTTCTAACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:944529/68‑1 (MQ=255)
                   aTCACTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGa                                     <  1:496832/67‑1 (MQ=255)
                       cTGGCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTTAATACAAACTTGCGGa                                     <  1:61593/63‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:2601830/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:269562/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:804200/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:1022667/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:2282546/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:1897111/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:1968239/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:1832928/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:1998254/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:1837027/68‑1 (MQ=255)
                        tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:1679693/68‑1 (MQ=255)
                          gCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCAt                            >  1:1502926/1‑69 (MQ=255)
                          gCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCAt                            >  1:825972/1‑69 (MQ=255)
                          gCATGTTTAACATGGTTTTTACATTTCTCACAGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCAt                            >  1:2954008/1‑69 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACa                                               >  1:2489201/1‑45 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACa                                               >  1:1738804/1‑45 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACa                                               >  1:2352486/1‑45 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACa                                               >  1:2705388/1‑45 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACa                                               >  1:382509/1‑45 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCAt                            >  1:1451141/1‑64 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCAt                            >  1:2013063/1‑64 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCAt                            >  1:863035/1‑64 (MQ=255)
                               tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                       >  1:2000167/1‑69 (MQ=255)
                                 tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCa                                 >  1:2722448/1‑57 (MQ=255)
                                 tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                       >  1:401235/1‑67 (MQ=255)
                                 tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                       >  1:2201160/1‑67 (MQ=255)
                                 tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                       >  1:1038481/1‑67 (MQ=255)
                                 tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                       >  1:2523421/1‑67 (MQ=255)
                                 tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                       >  1:2887560/1‑67 (MQ=255)
                                 tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                       >  1:2328161/1‑67 (MQ=255)
                                   aCATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAAcc                    <  1:1124976/68‑1 (MQ=255)
                                   aCATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAAcc                    <  1:3315141/68‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:2971245/37‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:3402938/37‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:1494704/37‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:1769202/37‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:1518669/37‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:463469/37‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACaa                                              <  1:1704828/37‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                               <  1:413065/52‑1 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTc               >  1:3106948/1‑68 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTc               >  1:776863/1‑68 (MQ=255)
                                        gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTc               >  1:2048131/1‑68 (MQ=255)
                                                aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt       <  1:1341441/68‑1 (MQ=255)
                                                aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt       <  1:2436984/68‑1 (MQ=255)
                                                aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt       <  1:1212632/68‑1 (MQ=255)
                                                aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt       <  1:1171312/68‑1 (MQ=255)
                                                aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt       <  1:2231849/68‑1 (MQ=255)
                                                aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt       <  1:3047724/68‑1 (MQ=255)
                                                aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt       <  1:428876/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:744703/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:1382957/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:2639243/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:1266334/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:2645747/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:1977462/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:3251802/68‑1 (MQ=255)
                                                     tcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:3225724/68‑1 (MQ=255)
                                                      cACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACgg  <  1:4300/67‑1 (MQ=255)
                                                         |                                                               
GCTCTTGCACTACCTTTGCATCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGG  >  minE/2245486‑2245606

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